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Chapter 1 



General Overview of the Thesis 



This thesis contributes to the field of Gravitational lensing (GL) and observational cosmology. 
We have investigated the possibility of detecting the existence of matter violating the weak 
energy conditions through its lensing effect on background sources. In a different context we 
have investigated GL in an alternative cosmology with a linearly evolving scale factor. We 
also considered gravitational lensing statistics in such a cosmology and its compatibility with 
existing observations. We have studied propagation of light in strong gravitational fields and 
the equations of motion for a vector particle with spin. Using geometrical optics approach 
we have shown that in Schwarzschild and Kerr geometries, massless particles deviate from 
geodesic motion. 

Einstein's General Theory of Relativity (GR) predicts that light rays are deviated from 
their straight path when they pass close to a massive body. This prediction was experimentally 
verified in 1919. Although the deflection is small, its effect can be enhanced by the passage of 
light over long distances. The deflected rays have enough time to intersect with one another to 
form caustics, and for objects which are compact and massive or are at cosmological distances 
there is a possibility of observing the effects of bending of light. Due to the bending of light, 
background objects appear distorted and, in extreme cases, form multiple images. This infor- 
mation can be used to obtain the distribution of mass in the lens in a completely novel way. 
The images of background objects are magnified by the action of lensing which makes them 
appear bigger (and therefore brighter). Thus, a gravitational lens acts as a natural telescope 
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providing us information about the distant objects which are otherwise too dim to be detected. 

GL is a powerful tool for exploring the universe. It can be used for the detection of exotic 
objects as well as for testing alternative theories of gravity. Proposals have been made to 
discover cosmic strings, boson stars, neutralino stars or wormholes through their gravitational 
lensing effects. There is no compelling evidence that any of the observed GL systems are due 
to these objects. However, it is essential to develop new lens models with objects which are 
not forbidden on theoretical grounds. 

The list of multiple-imaged gravitational lens systems has been growing steadily since the 
discovery of the first lens system in 1979 (the famous 'Old Faithful' QSR 0957+561 A&B). At 
present, more than 30 multiple-image systems are confirmed, or are very likely to be, as grav- 
itational lens systems. These lens systems can provide us with information about the universe 
as a whole. The global geometry of the universe, usually specified by its matter density and a 
cosmological constant, remains a significant source of uncertainty in modem cosmology. The 
possibility of using GL as a tool for the determination of cosmological parameters, either by a 
detailed study of specific lensing systems or through statistical analysis of samples of lenses, 
has been long and frequently discussed. One of the results of the previous works was that the 
mean image separations of lens systems have different dependence on source redshift in differ- 
ent cosmologies and that it may therefore be possible to measure the curvature of the universe 
directly. Besides, the expected frequency of multiple image lensing events for high redshift 
sources turned out to be quite sensitive to some cosmological parameters. All this makes the 
gravitational lensing statistics an interesting method to test different cosmological models. 

CHAPTER TWO is an introduction to the gravitational lensing theory. We give a brief 
introduction to the basic mathematical formalism for studying gravitational lensing and the 
background cosmology that we use to describe spacetime in which light propagation takes 
place. We introduce the terminology and concepts used in GL calculations, and review the 
basic equations of gravitational lensing. The basic types of lenses are presented and their 
properties are discussed. The second part of the chapter reviews different astrophysical and 
cosmological applications of gravitational lensing. 

hi CHAPTER THREE, the substance of which has appeared as Refs. [1,2], we discuss 
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the use of gravitational lensing as a tool in search of exotic objects in the Universe. In the first 
section we give an introduction to the energy conditions (EC) of classical GR and examine 
the consequences of their possible violations. Typically, observed violations are produced by 
small quantum systems and are of the order of h. One recent experimental study, however 
reported a violation, which could, in principle, arise because of the existence of classically 
forbidden regions carrying negative energy. It is currently far from clear whether there could 
be macroscopic quantities of such an exotic, EC-violating matter. 

Of all the systems which would require violations of the EC, wormholes are the most 
intriguing. The salient feature of these objects is that an embedding of one of their space- 
like sections in Euclidean space displays two asymptotically flat regions joined by a throat. 
Since wormholes have to violate the null EC in order to exist, the hypothesis underlying the 
positive mass theorem no longer applies, and there is nothing, in principle, that can prevent the 
occurrence of a negative total mass. In other words, we need to have some negative mass near 
the throat to keep the wormholes open. If wormholes exist, they could have formed naturally 
in the Big Bang, "inflated" from the "quantum foam" that is thought to underlie spacetime. 
Alternatively, they could have been constructed by an advanced extraterrestrial civilization as 
terminuses for, say, pan-galactic subway system. 

Discovery of any object with negative mass will not prove the existence of wormholes 
for sure, although it will certainly enhance the possibilities for wormholes to exist. In this 
chapter we analyze the gravitational effects light experiences while traversing the regions with 
negative mass. 

In the second section we provide an in-depth study of the theoretical peculiarities that arise 
in effective negative mass lensing, both for the case of a point mass lens and source, and for 
extended source situations. We describe novel observational signatures arising in the case of a 
source lensed by a negative mass. We show that a negative mass lens produces total or partial 
eclipse of the source in the umbra region and also show that the usual Shapiro time delay is 
replaced with an equivalent time gain. We describe these features both theoretically, as well 
as through numerical simulations. In the third section we provide negative mass microlensing 
simulations for various intensity profiles and discuss the differences between them. The light 
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curves for microlensing events are presented and contrasted with those due to lensing produced 
by normal matter. Presence or absence of these features in the observed microlensing events 
can shed light on the existence of natural wormholes in the Universe. 

In the last section we present a set of simulations of the macrolensing effects produced by 
large-scale cosmological violations of the energy conditions. These simulations show how the 
appearance of a background field of galaxies is affected when lensed by a region with an en- 
ergy density equivalent to a negative mass ranging from 10^^ to 10^^ solar masses. We compare 
with the macrolensing results of equal amounts of positive mass, and show that, contrary to the 
usual case where tangential arc-like structures are expected, there appear radial arcs — runaway 
filaments — and a central void. These results make the cosmological macrolensing, produced 
by space-time domains where the weak energy conditions is violated, observationally distin- 
guishable from standard regions. Whether large domains with negative energy density indeed 
exist in the universe can now be decided by future observations of deep fields. 

In the FOURTH CHAPTER we explore GL in an alternative cosmological model and the 
concordance of this theory with the current gravitational lensing observations. This chapter 
is divided into three sections. The first section introduces the problems standard model ex- 
periences and ways of their resolution. In the second section, part of which has appeared in 
ref. [3], we investigate the concordance of the gravitational lensing statistics with the cosmol- 
ogy in which the scale factor is linearly evolving. The use of gravitational lensing statistics as a 
tool for the determination of cosmological parameters either by a detailed study of specific lens 
systems or through a statistical analysis of a sample of lenses has been frequently discussed. It 
has been pointed out that the expected frequency of multiple imaging lensing events is sensitive 
to a cosmology. We use this test to constrain the power index a of the scale factor. We cal- 
culate the expected number of multiple image gravitational lens systems in a particular quasar 
sample with a known distribution of redshifts. This is compared with the observed frequency 
of lens systems found. Expected number of lens systems depends upon the index a through 
the angular diameter distances. We derive the expressions for angular diameter distances for 
this cosmology and use them in the lensing probability functions. By varying a, the number 
of lenses changes, which on comparison with the observations gives us a constraint on a. We 
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find that the value a=l corresponding to the coasting cosmology is in concordance with the 
number of observed lenses in the considered sample. In the last section, substance of which 
has appeared as ref [4], we introduce the non-minimally coupled effective gravity theory in 
which one can have non-topological soliton solutions. A typical solution is a spherical region 
having Gen = outside, and the canonical Newtonian value inside. Such a spherical domain 
(gravity -ball) is characterised by an effective index of refraction which causes bending of light 
incident on it. The gravity ball thus acts as a gravitational lens. We consider the gravity ball to 
be of a size of a typical cluster of galaxies and show that even empty (without matter) gravity 
ball can produce arc-like images of the background source galaxy. In the case of background 
random field the ball produces distortions ('shear') of that field. We also obtain constraints 
on the size of the large gravity ball which can be inferred form the existing observations of 
clusters with arcs. 

The FIFTH CHAPTER is dedicated to studies of the propagation of light in strong grav- 
itational fields. Most cosmological studies of lensing are performed in the weak field ap- 
proximation. However, there are interesting astrophysical situations where light propagates 
in a strong gravitational field. Weak field approximation becomes invalid in the vicinity of 
compact objects like black holes and pulsars. Thus, gravitational lensing can have additional 
effects. This chapter is divided into 5 sections. In the first two sections we present introduction 
to the question and motivations for the study. Third section describes the formulations of the 
theory, where the equation of wave propagation coupled to the curvature of the spacetime is 
derived. We also derive the modified geodesic equation and present a way to solve the equa- 
tion. Fourth section discusses the application of the results of the third section to the velocity of 
photons propagation in the field of the Kerr black hole. We conclude that the velocity remains 
subluminal contrary to recent claims. Fifth section, an edited version of Ref. [5], is dedicated 
to another approach to the problem and presents the derivation of the Papapetrou equation for 
massless particles from a simple Lagrangian. 

CHAPTER SIX is the concluding chapter and presents a summary and remarks with di- 
rections for the future work. 

Some technical details are included in the appendices. 
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Chapter 2 

Gravitational Lensing and Its 
Applications 

2.1 Gravitational lensing as a cosmic telescope 

"Do not Bodies act upon Light 
at a distance, and by their action 
bend its Rays; and is not this action 
strongest at the least distance?" 
I. Newton, Opticks, 1704 

Our Universe is controlled by gravity which doesn't limit its effects on matter — light rays 
can be deflected and bent. In 1704 Newton proposed that a light ray passing close to a massive 
body would be attracted and its path bent. In 1911 Einstein obtained the so-called "Newto- 
nian" value for the deflection angle from the Principle of Equivalence and the undisturbed 
Euclidean metric. However with the full equations of General Theory of Relativity (GR) Ein- 
stein obtained an angle twice the "Newtonian". The 1921 solar eclipse expedition confirmed 
the "Einstein" value, making the Einstein theory of General Relativity a new paradigm and 
Einstein famous. 

In fact, during those early years many scientists contributed to the subject of light deflec- 
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tion. It is worth mentioning such names as John Michell and H. Cavendish, P. S. Laplace, A. 
Eddington, who appear to be the first to point out that multiple (double) images can occur if 
two stars are sufficiently well aligned, O. Chowlson, who is actually responsible for what we 
now call the "Einstein ring". The first to mention that action of gravity of a massive body on 
light is similar to the refraction of light in an optical lens and called it "gravitational lens", was 
a Czech engineer R. Mandl in his letter to Einstein. In 1936 Einstein published a paper about 
this effect, where he remarked that "there is no great chance of observing this phenomenon". 
However, in 1937 the famous "prophet" of astrophysics, Fritz Zwicky, published a very opti- 
mistic paper about real possibility of discovery of gravitational lens (GL) in case of a galaxy 
lensed by foreground galaxy. He was the first to point out the usefulness of GL as an astro- 
physical tool which will allow a deeper look into the universe. He predicted many important 
applications of GL, pointed out that magnification leads to a selection bias and estimated the 
probability of detecting lensing to be very high. 

In the mid 60's the discovery of quasars (QSOs) renewed interest in GL. The subject was 
revived by S. Refsdal, who was later called by R. D. Blanford "the most reliable prophet in 
gravitational lensing". In his first paper [6] Refsdal gave a full account of the properties of 
the point-mass GL and calculated the time delay for the two images and mentioned a compact 
object as a candidate for a lens. In a subsequent paper he considered the application of GL 
for estimating the mass of the bending galaxy and the Hubble constant through observable pa- 
rameters of the source-lens system [7]. In 1965 it was suggested by J. Bamothy that QSOs are 
in fact Seyfert 1 galaxies, made to appear extra bright through other foreground galaxy acting 
as a gravity lens. It is now believed that QSOs are indeed the nuclei of galaxies, though not 
especially magnified by GL (with the possible exception of BL Lac objects, that can, in fact, be 
magnified QSOs). However, only after 1979 when the first GL system (QSO 0957+561 A & B, 
"The Old Faithful") was discovered did a systematic search for lenses begin. By now many 
GL systems have been identified which we can roughly divide into three classes: (i) more than 
30 proposed multiple images of QSOs; (ii) several tens of arcs and arclets; (iii) several radio 
and optical rings (and nearly rings) with one possible candidate for an X-ray ring [8]. One can 
now say that gravitational lensing significantly affects our view and physical understanding of 
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the distant Universe and of its major constituents. 

In spite of great theoretical and observational inspiration in GL field, there exist certain 
problems, as in any other new field. One of the basic problems is the amplification bias. The 
possible magnification of the source, associated with the deflection, can fool the observer: 
some QSOs seen through a foreground galaxy are much fainter than they appear to be. If 
one conducts a flux-limited QSO search, some of these sources get boosted in flux above 
the threshold of the sample. The net effect is that frequency of multiple imaging appears 
much greater in flux-limited samples than in volume-limited ones. Another direct consequence 
of the amplification bias is that one would naturally expect to find an excess of amplifying 
galaxies near distant and bright QSOs selected from a flux-limited sample. This has indeed 
been reported [9]. 

Another problem is the "verification", or "when is a lens a lens". This became important 
after the discovery of several binary QSOs. There exist several basic characteristics of images 
that are signatures of gravitational lensing, though this list is by no means exclusive: (i) mul- 
tiple images of the same object; (ii) background image seen as a nearly complete ring or as an 
extended arc; (iii) spectroscopic similarities of the images; (iv) detection of the lensing galaxy 
or cluster in the right location and of sufficient mass to create the image splitting; (iiv) images 
exhibiting brightness variations characteristic of compact objects crossing the line of sight. 

However, in spite of the problems and the uncertainties, the theoretical and observational 
achievements of last few decades have made the GL one of the most active and exciting fields 
of research. All possible applications of gravitational lensing are surely impossible to list. 
From determination of the Hubble constant and mass of the bender to using microlensing 
for dark matter search to testing alternative gravitation theories and exotic objects search — 
there is hardly an area of cosmology and astrophysics where GL has not been applied. In 
this Introduction we will concentrate on the basic theory underlying the gravitational lensing 
subject and the new progress of GL as well as describe some of its interesting applications in 
astrophysics. 
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2.2 Elements of GR and cosmology and propagation of light 



2.2.1 Basic notions of GR 

The basic principles of Einstein's GR are considered to be the Principle of Relativity and 
Principle of Equivalence (EP). However, the decisive step for the construction of the general 
relativity formalism was the suggestion by Einstein and Grossman in 191 2- 1913 that the grav- 
itational field must be identified with the non-Minkowskian metric of the spacetime. Usually 
this suggestion is deduced from the EP. However, some authors consider the metricity of grav- 
itation as the independent, if not, the main, principle of GR [10, 1 1]. If we accept the metricity 
of the gravitational field, then a pseudo-Riemannian spacetime must be chosen as the model 
for our spacetime. (It is characterised by the structure of the manifold, pseudo-Riemannian 
metric, connection and curvature). 

Differential manifold We will assume that the spacetime has the properties of a continuum, 
i.e. it is a four-dimensional differential affine manifold X4. Any point of it can be labeled by 
real coordinates x'^ with = 0, 1,2, 3, where refers to the time coordinate, and 1,2,3 to the 
space coordinates. Any affine space has a defined notion of the interval between its points. At 
any point there exists an independent metric tensor field gfj,^, = gnu(x) in order to allow for local 
measurements of distances and angles. The square of the infinitesimal interval ds between 
and x^ + dx^ is then determined by 

ds^ = g^^,dx^dx'' . (2.1) 
Thus, we can adopt the definition [12]: 

Definition 1 A spacetime is a four-dimensional manifold equipped with a Lorentzian (pseudo- 
Riemannian ) metric with signature (-,+,+,+). (The manifold should be Hausdorjfand para- 
compact.) This definition is important in the studies of wormholes (see Chapter 3.1). To admit 
the construction of a wormhole, a spacetime shall be either "almost-everywhere Lorentzian" 
or be non-Hausdorff. 
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Connection In order to do physics in such a spacetime, we should have additional structures 
onXi. In an Xi it does not make sense to say "a (nonzero) vector field is constant." To give such 
a statement a meaning, one must introduce the notion of parallel transfer of vectors. Parallely 
displaced from and x^+dx^, a vector C changes according to the prescription 

rfC^ = -r^^(x)C"J;c^ . (2.2) 

Here dC^ is assumed to be bilinear in C" and dxi^, the set of the 64 coefficients V^p is the 
affine connection. An Xt equipped with a F is called a linearly connected space or L4. The 
parallel transport law (2.2) can be extended to higher rank tensor fields and densities, and it 
is possible to define their covariant differentiation with respect to F; for a vector, V^C^ = 
duC^ + T'^^C°' . Now it makes sense to state that "a field is constant over spacetime" — its 
covariant derivative has to vanish. In general, metric and connection are two independent 
geometrical constructions on a manifold. But in GR, since the gravitational field is identified 
with the metric, a connection should satisfy two additional conditions: 

• Symmetry Kv = '^ta^ 

• Metricity ^vgap = '^- 

The last postulate guarantees that lengths, in particular the unit length, and angles are 
preserved under parallel displacement. This enables a choice of a coordinate system with 
connection coefficients vanishing at a point x and metric g^uix) = r]^^ (the Minkowskian 
metric) at this point. Such a coordinate system is called the locally inertial system at x. 

The above conditions allow us to express the components of such a connection through the 
metric: 

Components (2.3) are called Christoffel symbols. 



Curvature Parallel transfer is a path-dependent concept. If we parallely transfer a vector 
around an infinitesimal area back to its starting point, we find that its components change. 
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This change is proportional to the Riemann curvature tensor 



Physical meanings Loosely speaking, the connection governs the "acceleration" of a freely 
falling particle in a gravitational field. A small test particle, free to move under the unfluence 
of gravity alone, will follow a geodesic of the metric.^ According to (2.1), the length between 
two given points depends only on the metric field. Therefore, the differential equation for the 

1 /2 

extremals can be derived from 5 Jds = S J (^g/^udx^dx^) = and results in 

dhi" ^„ dx"" dx^ ^ 

the equation of geodesic. The Riemann tensor governs the difference in acceleration of two 
freely falling particles that are near to each other. 

Einstein equations By contracting the Riemann tensor we obtain 
• Ricci tensor 

D — E) a/3 



Ricci scalar 



a/3 



and Einstein tensor 



Finally, the Einstein field equations relate the curvature of spacetime (as measured by the 
Einstein tensor G^^) to the distribution of matter and energy (as measured by the stress-energy 



'a correction to this definition is in order (the meaning will become more obvious in Chapter 5): geodesic 
equation is an equation for motion of a moving free spinless test particle in an external gravitational field (for ex. 
[13]. 
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tensor T^^). Explicitly, 

SttG 

The components T'^, T'j are, correspondingly, energy density, the energy flux and the 
stress. For most astrophysical and cosmological purposes, one idealises bulk matter as a "per- 
fect fluid", for which 

= {pc^ + p) u'^u'^-pg'^l^ , (2.7) 

where p denotes mass density and p the pressure, both measured by a comoving observer, and 
is the 4-velocity, normilized to unity 

^„;3mV = 1. (2.8) 
The Einstein equations can be derived from the Einstein-Hilbert action principle 

where Lm is the matter Lanrangian. 



2.2.2 Weak fields 

Metric If the gravitational field is "weak" — the metric of the spacetime differs little from the 
flat, Minkowskian metric, one can write in approximately Cartesian coordinates 

g/ii/ = ?7m<^ + V' I VI << 1 • (2-10) 

GR then can be reduced to a "linearized theory", where the metric for any static distribution of 
matter can be given as [14] 




(2.11) 
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where is the Newtonian gravitational potential and \dx\ = dl denotes the Einstein spa- 
tial line element. This "post-Minkowskian " metric satisfies the weak-field condition 2.10, if 
|$| «c^ and matter moves slowly |v| «c. 

Effective refraction index From 2. 1 1 , by putting ds = and solving for dl /dt, we can obtain 
an effective speed of light, v^n, to the first order 



The speed of light as measured in a locally inertial frame is c, but since the coordinate system 
(?,x) is not an inertial frame, the apparent coordinate speed of light is different from its value 
in vacuum. We can characterize the effect of light propagation in the presence of gravitational 
potential by the effective refractive index n given by 



The gravitational potential for a massive object is a negative quantity, therefore, the apparent 
speed of light is slower in the presence of a gravitational field. We assume that the lens is 
stationary, therefore its gravitational field is a function of only space and is independent of 
time. The light rays effectively move through a region of space with spatially varying refractive 
index. This causes bending of light in analogy with the usual optical phenomenon. We note 
that this effective refractive index is independent of the wavelength of light, and, thus, to a 
very good approximation, gravitational lensing is achromatic. 



The full Einstein equations are nonlinear. This is the principal difficulty in extracting the exact 
solutions. Nevertheless, many solutions have been found; among them are the Schwarzschild 




n(x) = — = 1 - 



2$(x) 
~c^ 



(2.12) 



2.2.3 Strong fields 



, the Friedmann-Robertson- Walker (FRW) and the Kerr solutions. 
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ADM split In any well-behaved coordinate patch one can use the "time" coordinate to de- 
compose the (3 -I- l)-dimensional Lorentzian metric via the Arnowitt-Deser-Misner (ADM) 
split [15]. The ADM split yields: 



(2.13) 



Sij 



The function N(t,x) is known as the lapse function, while /3(t,x) is known as the shift func- 
tion. The three-metric gij(t,x) describes the geometry of "space", while the lapse and the shift 
functions describe how the space slices are assembled to form a spacetime . This ADM split 
allows one to adopt quick and dirty definition of a horizon. In every asymptotically flat region 
^ 1, /3 — > and gij Sjj asymptotically as one approaches spatial infinity. Associated with 
each asymptotically flat region one may define a putative horizon by vanishing of the lapse 
function. Roughly speaking, when = time has slowed to a stop. ADM split is essential to 
define the ADM mass (see App. A). 

2.2.4 Standard cosmological model 

FRW Universe Our universe is believed to be homogeneous and isotropic on large scales. 
This is borne out by the observations of the distribution of galaxies and the remarkable isotropy 
of the Cosmic Microwave Background Radiation. These assumptions, together with the field 
equations of GR, give solutions for the geometry of the universe. The space-time metric for 
such a universe is given by the FRW line element (for example, [16]) 



where (r, 9, cf)) are the spatial comoving spherical polar coordinates of a space-time point and t 
is the cosmic time, a{t) denotes the scale factor of the universe at time t and k is the curvature 
index of the spatial hypersurfaces t = constant. In this form of the metric we can rescale the 
coordinates in such a way that constant is -l-l, -1 or 0, corresponding to spatial sections of 



ds^ = c^df--a^(t) 



+ r\sin^ede^ + d(f)^) 



(2.14) 
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constant positive, negative or vanishing curvature, respectively. With such a rescaling, the 
coordinate r in the metric is dimensionless and a{t) has dimensions of length. The dynamics 
of the universe is obtained from the field equations of GR. These equations are presented later 
in this chapter. But many properties of the universe, which are kinematic in nature, can be 
obtained solely from the FRW metric. 

Light from distant objects appears redshifted due to the expansion of the universe. The 
observed redshift is defined as z = (Ao-Ae)/Ae, where Ae and Aq are the emitted and the 
observed wavelengths, respectively. It is related to the expansion parameter by 

^ = l+z, (2.15) 
a 

where uq is the present value of the scale factor and a is the value of the scale factor when the 
light ray was emitted from the source. We define the Hubble parameter H{t), which measures 
the rate of change of the scale factor at any time t as H{t) = d{t)/a{t). The present value 
of the Hubble parameter is denoted as Hq. Its latest numerical value is ascertained to be 
H(, = 70±7 kms~^Mpc~^ [17]. The dynamics of the universe depends on the matter content of 
the universe. This can be specified by the energy density p(t) and the pressure p(t), which are 
often related by an equation of state of the form p = wp; the classic examples are 



• Non-relativistic matter: (w = 0) Galaxies are the tracers of the expansion of the uni- 
verse in the sense that they follow the general expansion of the universe. Treated as a 

fluid, they exert negligible pressure, therefore, to an excellent approximation, they can 
be treated as pressureless dust. 

• Radiation: (w = 1 /3) A major component of the early universe was in the form of 
the radiation. It is believed that after the radiation era, the universe has undergone a 
long period of matter domination, with the radiation providing only 10^ of the closure 
density today. 



Cosmological Constant (A): (w = -1) Though originally introduced as an arbitrary 
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constant by Einstein ^, it has made a comeback in the recent times [18]. The current 
observations suggest that about 2/3 of the present day density is in the form of cosmo- 
logical constant A, and only 1 /3 is in the form of matter. 



In general though there need not be a simple equation of state. There maybe more than one 
type of material, such as combination of radiation and non-relativistic matter. Certain types of 
matter, such as a scalar field, cannot be described by an equation of state at all. 

Dynamical Equations The crucial equations governing the evolution of the scale factor of 
the universe and the matter-energy content of the universe, are the Friedmann equations: 

— + 5 Ac2 = -p (2.16) 

dV kc^ SttG Ac^ 



a J 3 3 

where we have also included a constant A term. The spatial geoometry is flat if ^ = 0. For a 
given H, this requires that the density equals the critical density 

= 3//oV87rG , 

Densities are often measured as fractions of p^: 

0(0 = P- , 

Pc 

the dimensionless density parameter as Om = Po/Pc and the Lambda parameter as ilLambda = 

Distance Measures In an expanding, curved space-time, the measure of distance is not 
uniquely defined. The distance between two points can be defined by: 



^The original motivation was that the universe was beheved to be static and therefore a cosmic repulsive 
force was needed to balance the attractive force of gravity. This was later on abandoned by Einstein when it was 
discovered that the universe is, in fact, expanding. 
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• The light travel distance. 

• The flux received from a standard candle. 



• The angle subtended by a standard ruler. 

In a Euclidean space all these three measures would coincide. However, in a non-Euclidean 
space-time the three measures are all different. We need to define all the distances separately 
in complete analogy with their Euclidean counterparts. 

(a) Comoving Coordinate Distance: A quantity of interest is the coordinate distance r(z) 
up to a redshift z. Since light rays travel along the null geodesies of the space-time, 
ds = 0. For a FRW metric we obtain 



We can convert the integrals over t into integrals over z by differentiating \+z = ao/ato 
obtain dz/dt- -{\ + z)H{z) and substituting in the previous equation. This gives 



where we defined the dimensionless Hubble parameter h{z) = H{z)/Hq. For a flat space- 
time, where ^ = 0, we obtain 



(b) Luminosity Distance: The Luminosity distance Dl is defined in such a way as to pre- 
serve the Euclidean inverse-square law of diminishing of light with distance from a point 
source. This gives 




(2.18) 




(2.19) 




(2.20) 



Di^^aQr{\ + z) . 



(2.21) 



(c) 



Angular Diameter Distance: The Angular Diameter Distance Z)a is defined in such a 
way as to preserve a geometrical property of Euclidean space, namely, that the angular 
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size subtended by an object should fall off inversely with dp^. This gives 



aor 



(2.22) 



l+z 



In later sections we will use the angular diameter distance measured by an observer at zi 
up to zi- This is given by 



which gives 



2.3 Basic concepts of gravitational lensing 

2.3.1 Approximations 

In the following section we describe results on gravitational lensing based on the Refs. [14, 19- 
23]. The formal description of GL is based on several approximations. In the previous section 
we have already made an assumption of weak gravitational fields, thus, justified the use of 
linearized field equations of GR. Indeed, even in clusters of galaxies, the deflection angles are 
well below V and the maximum image separation in multiple-imaged systems are not more 
than 10". Thus, we can express the approximations used as 

• We assume that the gravitational field can be described by the linearized metric: 



where $ is the Newtonian potential due to the gravitational field. 

• Geometrical optics approximation — the scale over which the gravitational field changes 
is much larger than the wavelength of the light being deflected. 




(2.24) 




• Small-angle approximation — the total deflection angle is small. The typical bending 
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angles involved in gravitational lensing of cosmological interest are < 1'; therefore we 
describe the lens optics in the paraxial approximation. 

• Geometrically-thin lens approximation — the maximum deviation of the ray is small 
compared to the length scale on which the gravitational field changes. Although the 
scattering takes place continuously over the trajectory of the photon, appreciable bend- 
ing occurs only within a distance of the order of the impact parameter. 



Figure 2.1: Schematic diagram of gravitational lensing. The two hypothetical planes are de- 
fined to erect a convenient coordinate system. 

We begin the discussion of gravitational lensing by defining two planes, the source and the 
lens plane. Fig. 2. 1 describes a typical lensing situation. A convenient origin, passing through 
the lens is chosen on the sky. The planes, described by Cartesian coordinate systems, pass 
through the source and deflecting mass and are perpendicular to the optical axis (the straight 
line extended from the source plane through the deflecting mass to the observer). These planes 
are hypothetical and are solely for the purpose of visualization. The coordinates of the image 
with respect to the origin are (^1,^2) and that of the source are (7^1,772), respectively. Since 





Source ^PQme 
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the components of the image and the source positions are much smaller in comparison to the 
distances to lens and source planes, we can write the coordinates in terms of the observed 
angles. Therefore, the image coordinates can be written as (61,62) and those of the source as 

(A,/?2). 

2.3.2 Deflection angle 

To calculate the bending angle produced by a gravitational field, we use the ray equation which 
describes the path of a light ray through a spatially varying refractive index n(x). To derive the 
ray equation we start from Fermat's principle, which states that the light travel time from the 
source to the observer is an extremum. 



[•to 

> / dt = Q, 

Jts 



(2.25) 



where the subscript "5" stands for the source and "o" stands for the observer. The integral is 
evaluated along the trajectory of the light ray. To obtain the ray equation we parameterize the 
ray path by s, therefore 



dx dx ds 
dt ds dt 



(2.26) 



Substituting \dx/dt\ = c/n{x) in Eq. 2.25, we obtain 

b / n(x) Y^{dxi/dsf ds = Q. (2.27) 

J Sem \ ; I 



1/2 



Using the Euler-Lagrange equation and choosing the parameter s to be the path length / we 
obtain 

^[n(x)£]=V/i, (2.28) 

where x = dx/dl is the unit tangent vector along the path of the ray x(/). In astrophysical 
applications of GL the bending angles are small, therefore, to obtain the deflection angle using 
this formula we can integrate equation (2.28) along the unperturbed path of the ray 7 to obtain 

ocd = 0^em -%bs)± =- V±ndl , (2.29) 
J J 
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where the component perpendicular to the unperturbed ray is used in the calculation. Using 
this formula we can obtain the deflection angle angle due to a point mass. We set our coordinate 
system such that the source and the lens lie along the z-axis, and the origin is chosen at the 
position of the lens, the component of the position vector perpendicular to the z-axis being 
denoted as ^. The impact parameter is the distance of the unperturbed ray from the centre 
orthogonal to the direction of propagation. From (2.29) we obtain the deflection angle as 



ai = -2GM I Vg| J__, | Jz=^ / |l + (f)l (2.30) 



/ 




where, since the main contribution to the integral comes from the range -,^o <z < we have 
put the limits of the integral as minus and plus infinity. Thus, the Einstein deflection angle of 
a light ray passing near a compact mass M at a distance $, is 

4GM 

«d(0 = — ^- (2.31) 

This bending angle is twice the value of what would be expected from the Newtonian theory. 
Einstein's General Theory of Relativity was vindicated when this angle was measured for the 
case of the Sun, where the predicted value (~ 1.75") was confirmed by observations during a 
total solar eclipse in 1919. 
For an extended mass we can obtain this angle by integrating individual deflections due to 
all mass elements constituting the lens. In thin-lens approximation the deflection angle 
can be obtained by projecting the volume mass density of the deflector onto the lens plane 
€ - (6) '^2), which results in a surface mass density S(^) = J p{^,z)dz. The deflection angle 
is a superposition of Einstein angles for mass elements dm = E(^) J^^. Considering all the 
deflecting mass to be concentrated only in the lens plane and the deflection taking place only 
in the lens plane, a deflection angle can be expressed as a two dimensional vector 

.2 I I |^_^/|2 



ccM) = -r 7. l,n d'i' . (2.32) 
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From here we obtain 



where we have used the fact that ^ = Di^d and defined displacement vector a as 

a = -^ad. (2.34) 

There is no unique definition of distances in a curved spacetime. Distances which should 
be used in this equation are the angular diameter distances to ensure that the equation 
remains valid for a more general spacetime (see 2.2.4). Defining the critical density as 
Scr = (c^/47rG)(Ds/Z)L^Ls) and the dimensionless quantity K(d) - Il(0)/Scr, we can write 
equation (2.33) as 

«=^//-(^')t^^'^'- (2.35) 



2.3.3 Lens equation and the lensing potential 

Considering the projection of the light ray on the two planes, we can derive a relation between 
the source coordinates and the image coordinates in terms of the bending angle 

/3 = 0-^ad, (2.36) 

Equation (2.35) can be written as 

a = Vip, ip = ^ JJ K(d')ln\d-e'\d^e' . (2.37) 

Using the identity V^ln |x-x'| = 2-k6^(x-x') we obtain the equation which the dimensionless 
relativistic lens potential i/j satisfies: 



V^(d) = 2K(d) 
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(2.38) 



In terms of this potential the lens equation (2.36) can be written as 



(3 = d-Vip{d). 



(2.39) 



In general, angles /3, ^, a and may not be coplanar and so, Eq. 2.39 is a vector equation. 
Given the matter distribution of the lens and the position of the source the lens equation may 

have more than one solution, which means that the same source can be seen at several positions 
in the sky. The lens equation describes a mapping 6 /3 from the lens plane to the source 
plane. 

2.3.4 Magnification 

Besides multiple imaging, the differential deflection across a light bundle affects the properties 
of the images. In particular, the cross-sectional area of the bundle gets distorted and the flux of 
the images is influenced. The source subtends a solid angle Aco^ at the observer in the absence 
of lensing. In the case of lensing AuJi is a solid angle subtended by the image. Gravitational 
lensing preserves the surface brightness of the source (we assume that during deflection no 
absorption or emission of light is taking place and that deflection by a nearly static deflector 
introduces no additional frequency shift between the source and observer, except a cosmolog- 
ical redshift). The flux is 5 = 7 • Ao;, I being specific intensity. For an infinitesimally small 
source the ratio between the solid angles gives the flux amplification due to lensing 



A solid-angle element Sp'^ of the source is mapped to the solid-angle element 59^ of the image, 
and the magnification is given by 




5, dui 



Local properties of the lens mapping are described by its Jacobian matrix A 




(2.40) 
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The Jacobian matrix A is thus the inverse of the magnification factor 



det^ 

de 



^ (2.41) 



det^(6>) 




Eq. 2.40 shows that the matrix describes deviation of the lens mapping from the identity 
mapping. From (2.38) we have 

Tri^ij = 2k . (2.42) 
Two additional combinations of V'o are important: 

71 = ^(^11-^22); 

72 = "012 = -021 • 

With these definitions we can write the Jacobian matrix A as 

. l-K-71 

^=1 I ; (2.43) 

-72 

and the magnification factor 

The eigenvalues of A are ai 2 = 1 - k ^ 7, where 7 = V7f+72 determinant is det^ = 

(l-K^)- 7^. When the line of sight competely misses the deflector, E(0) = and the k term 
vanishes in Eq. 2.38. So, k represents the amplitude of the convergence due to the matter 
within the light-ray (also called Ricci focusing), while the 7 term is the amplitude of the shear 
due to the matter outside the beam (also called Weyl focusing). Eigenvalues of A describe 
the image distortion in the radial and tangential directions for a circular source, resulting in an 
ellipse. 

The zeroes of the Jacobian of the lens mapping are called the singular points of the lens 
mapping. For the isolated lenses the lens mapping would go to identity at large distances from 
the lens mapping. For such lenses the zeroes of the Jacobian are either points or closed curves 
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in the image plane. These curves are called the critical curves. Their images in the source 
plane are calles the caustics. Caustics separate the regions of different image multiplicities. 
When a source crosses a caustic the number of images changes by two. 



2.3.5 Lens models 

The mass distribution inside galaxies and clusters of galaxies is in general quite complicated 
and may not have any symmetry. However, since the circular mass distributions are easier to 
handle analytically they are very convenient to use in gravitational lensing. Symmetry allows 
the lens equation to be separated in the polar coordinates making the equations analytically 
tractable. Besides, for many celestial bodies, like "Jupiters", stars, black holes, and even 
galaxies when the light rays pass outside the deflector, a point mass approximation is valid. 
Galaxies and even clusters of galaxies are also well approximated by the singular isothermal 
sphere. Below we give a description of these two models with corresponding lensing equa- 
tions. 

Point Mass (Schwarzschild) Lens^ 
(/) Lens equation 

Due to axial symmetry, the propagation of light reduces to one dimensional problem. Let us 
consider the situation described in Figure 2.2. 
From the Figure follows the geometric relation 

eOs^pDs + aaDi^s- (2.45) 
Substituting for a from (2.31), we rewrite it as 

Dls 4gm_ 



^In lens theory the term "point mass" is used whenever one is concerned with light rays deflected with the 
impact parameters greater than the Schwarzschild radius of a static spherical object; the exterior of such an object 
is always described by the Schwarzschild metric, hence the term Schwarzschild lens. 
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Figure 2.2: Basic geometry of gravitational lensing 



where 




4GMDls 



(2.47) 



Angle 6'e denotes the angular radius, called the Einstein radius. It provides a natural angular 
scale to describe the lensing geometry. Sources which are closer than about 6e to the optic 
axis are significantly magnified, whereas sources which are located well outside the Einstein 
ring are magnified very little. Besides, it is the radius of a tangential critical curve. In the 
given case of, the caustic is a point on tye optical axis. If a source is displaced slightly off the 
axis, two bright images are created on opposite sides of the lens centre, one just inside, and the 
other just outside the critical radius. The equation (2.46) has two real roots: 



which correspond to two physical images of the source S. The angular separation between the 
images is 



The separation between the source and the deflector is related to the image position by 




(2.48) 




(2.49) 
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(2.50) 



Thus, the lens equation has two solutions of the opposite sign. The source has an image on 
each side of the lens, one inside Einstein radius, one outside. If the source is a disk of radius (^s, 
the images will represent ellipses, squeezed along the axis connecting them and stretched in 
perpendicular direction. For example, if a, b — semi-major axes, the area of ellipse is Q = i^ab. 
The relations with the radius of the source we can write as 



a _9i b _ dOi 

ips /3' (ps dp ' 



(2.51) 



In the case of perfect alignment between source, lens and observer, (/? = 0), an observer will 
see a ring with radius 9e and thickness 



1 1 

+ 



2 2^/;5^T4^ 



= ^2(/?s = (/9s , 



(2.52) 



equal to the source radius. The solid angle which it subtends on the sky is then IttOe^Ps 
(ii) Magnifications 

For a circular symmetric lens, the magnification factor /x (Eq. 2.41) is reduced to 



11 = 



e_de_ 

pdp 



(2.53) 



For a point mass lens, which is a special case of a circular symmetric lens, we substitute for /? 
using the lens equation (2.46) to obtain the magnifications of the two images. 
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+ 2 ^1 
2uVu^ + 4 2 



(2.54) 



where u is the angular separation of the source in units of Eistein angle. The total magnification 
of the two images is 



u^ + 2 



(2.55) 



M = 1, or = 6'e, is often taken to be a typical case that characterizes the efficiency of the 
lens. This corresponds to // = 1.34 or A/j, = 0.32 in apparent magnitude. Point-like masses 
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play an important role in the study of microlensing, which arises when the separation of the 
images is too small to be resolved and the lensing effect can only be observed through the 
lensing-induced time variability of the source. 

Singular Isothermal Sphere 

When we consider galaxies as lenses we need to allow for the distributed nature of the matter. 

The studies of the flat rotation curves of galaxies and the galaxy /gas distributions in clusters of 
galaxies suggest that the total matter profiles in these systems follow the singular isothermal 
sphere (SIS) model very well 

2 

p(r) = , (2.56) 

where measures the line-of-sight velocity dispersion. In this model it is assumed that the 
mass components behave like particles of an ideal gas, confined by their combined spherically 
symmetric gravitational potential. It is assumed also that the gas is isothermal, so that is 

constant across the galaxy mcr^ = kT . Upon projecting along the line-of-sight, we obtain the 
surface mass density 

E(0=^. (2.57) 

where ^ is the distance from the centre of the two-dimensional profile. Referring to Eq. 2.31, 
we find 



a^ = An^, (2.58) 
which is independent of the impact parameter. The Einstein angle in this case is 

6'e = 47r^— = ttd-p— • (2.59) 
Ds Ds 

The solution to the lens equation (3 = 9-9^is 

e± = eE±f3 for/3<^E- (2.60) 



Thus, the lens has two images on the oposite sides of the lens centre. For (3 > 6^ only one 
image appears at 9 = 6+ = /3+9e. The image separation is just the diameter of the Einstein ring: 
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A6 = 6+ + 9- = 26e. Although the surface mass density is infinite at ^ = 0, the behaviour of 
the model for larger values of ^ seems to approximate the matter distribution of galaxies fairly 
well. Real galaxies, however, cannot follow the density law (2.56) due to an infinite density at 
the centre and an infinite mass. Other models exist and are frequently employed. For example, 
if the singularity is removed from the centre, the model is called a softened SIS — a SIS with 
a finite core radius (ISC). In this case a lens is capable of producing either one, or three 
images. However, the usual absence of the third image implies that even if there is a core, it 
must be small, <, 200 parsecs [24]. 

The magnifications of the two images follow from Eq. 2.53 and are 





1^1 




1^1 




Oe 



the circle |^| = 6'e is a tangential critical curve. Images are stretched in the tangential direction 
by a factor |//|, whereas the distortion factor in the radial direction is unity (see Section 2.3.4) 

2.4 Astrophysical applications of gravitational lensing 

Listing already existing and possible future astrophysical and cosmological applications of 
gravitational lensing is nearly an impossible task, so vast has become this field in the last 
years. Zwicky's idea of the gravitational lensing as a cosmic telescope is proving itself with 
each observational discovery. We can see the magnified distorted images of galaxies which 
otherwise are far too dim to be observed, not to say, studied. Gravitational lensing effect 
allows us to test the General Relativity Theory, to probe the nature of the lensing object, the 
source and the intermediate space, and to test the large-scale structure of the universe. We will 
describe several interesting applications of gravitational lensing. 

2.4.1 Determination of Hubble parameter and mass of the deflector 

Hubble parameter One of the first applications of gravitational lensing, suggested by Refs- 
dal [7], is the determination of the Hubble constant via the direct measurement of the time 
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delay A? between the observed light curves of multiply imaged quasars. For axially symmet- 
ric lens this method can be described using the wavefront picture. Wavefront characterizes 
the locus of all points with equal light-time-travel from the source. Light rays in vacuum are 
perpendicular to the wavefronts, which are spherical close to the source. However, they be- 
come deformed by the gravitational field of deflector, may intersect themselves and cross the 
observer several times, producing multiple imaging. Every passage past the observer corre- 
sponds to the image of the source in the direction normal to the wavefront. The time delay for 
pair of images is the time between two crossings of the wavefront. 

The wavefronts from a distant, doubly imaged quasar cross each other at the symmetry 
point. They represent the same light propagation time and for an observer, located at a distance 
X from the symmetry axis, the time delay must be equal to the distance between the wavefronts 
at the observer divided by the velocity of light. The deflection law can be written as [7] 

«dOc|erS (2.62) 

with e = for a point mass lens. Using Hubble relation for small redshifts 

cz = DHo, (2.63) 

one obtains the expression for the Hubble parameter Hq in terms of observable quantities 

HoAt = ^-^^AeH2-e) . (2.64) 



Determination of mass One more direct application of gravitational lensing is the determi- 
nation of the mass of the deflector. The simplest situation here is when the lens is a spherically 
symmetric object and the source lies exactly behind the lens centre. The lens can then form an 
Einstein ring. In this case, the bending angle is 

4GM« 9) 
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and lensing equation with the source at the origin becomes 



e = am 



(2.66) 



Combining these two equations we obtain 



M« 9) = 7r(DL^)^S, 



(2.67) 



The mass inside the Einstein ring can be determined, once its angular diameter and redshifts 
of the lens and the image are known. Even if the alignment of the source, deflector and the 
observer is not perfect, and the ring is not observed, this mass estimate may be very useful 
and rather accurate. For example, the mass inside the inner 0".9 of the lensing galaxy in the 
quadruple quasar QSO 2237 + 0305 ("Einstein cross") has been determined with an accuracy 
of a few percent [14], with the largest uncertainty being due to the estimate of the Hubble 
constant. This method doesn't depend on the nature or state of matter, but it measures only the 
projected mass and only in the inner part of the lensing galaxy. 

Another method of determining the mass of the galaxy-deflector is from making use of the 
Eq. 2.64. From the observed image separation ~ 26^ (Eq. 2.49) we find 



Thus, the observed image separation oc MHq, if the redshifts of the lens and the source are 
known. Using the relation between At and Hq (2.64) we obtain M oc At. Thus, from the direct 
measurement of the At one can determine the mass M of the galaxy-lens located within an 
angular radius A9/2. 

Determination of the mass and mass distribution of the cluster of galaxies has become 



2^E ^ 01 + 02 



So that. 




(2.68) 
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posible since the discovery of arcs and arclets. Arcs are the result of very strong distortion of 
background galaxies (when a part of an extended source covers different parts of the diamond- 
shaped caustic, associated with cluster of galaxies as a deflector). Assuming the cluster mass 
distribution to be axially-symmetric, we can have a rough mass estimate from (2.67), where 9 
now is the distance of the arc from the cluster centre, and Ecr can be determined, if the redshift 
of the arc can be determined. Also, since the arc roughly traces the Einstein radius, we can use 
the Eq. 2.47. This method loses accuracy if the cluster is highly asymmetric or has significant 
substructure (clumps of dark matter, for example). However, it is believed that the assumption, 
describing clusters as isothermal spheres with finite cores, works well [21]. 

Discovery of arclets (less elongated images of background galaxies than arcs) and weakly 
distorted images of background galaxies opened up the possibility of studying the mass dis- 
tribution in the outer parts of the clusters. Shape of a galaxy image is affected by the tidal 
gravitational field along its corresponding light bundle. This distortion is small and since 
galaxies have intrinsically different shapes, the effect cannot be determined in the individ- 
ual galaxy image. However, with the sky densely covered by randomly oriented faint galaxy 
images, a statistical study of the distortions of these far-away sources is possible. From the 
coherent alignment of images of an ensemble of galaxies one can draw the distortion pattern 
which traces the gravitational field of the foreground cluster. By reconstruction techniques one 
can measure the tidal field related to the gravitational potential of the cluster and obtain the 
surface mass density. 

2.4.2 A candidate string-lensing field 

One interesting aspect of the lensing phenomena is lensing by a straight cosmic string. Grav- 
itational interaction of strings is characterized by the parameter Gji, where ji is mass per unit 
length of the string n ~ rf. Here r] is the energy scale of symmetry breaking ~ 10^^ GeV for 
grand unification scale. Thus, 

G//~(77/mpi)2, (2.69) 
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Figure 2.3: The conical space around a straight string can be obtained from Eucledean space 
by cutting out the wedge of angular size SnGfi and identifying the exposed surfaces. Light rays 
emitted by quasar intersect behind the string and observer O sees a double image of quasar Q. 
If I and d are the distances from the string to the quasar and to the observer, respectively, then 
the angular separation between the images is 5 = S7TGnl(d+l)~^. 

nipi — Planck mass ~ 10^^ GeV and G = m~f. For grand unification strings Gfi ~ 10"^ [25] . The 
metric around the straight string is flat, so it cannot be detected by gravitational interaction. 
However, the space around the string is actually a conical space, that can be made out of a 
Euclidean space by cutting out a wedge of angular size 5 and by identifying the opposite sides 
of the wedge (Fig. 2.3). 

The deficit angle 5 is 5 = SnGfi. As seen from the Figure 2.3, the conical nature of space 
around the string can give rise to lensing phenomena. If there is an intervening string between 
us and background galaxy, such a string should produce an identical twin pair of images of the 
galaxy over a strip of space owing to an angle deficit 5 around the string. The discovery of the 
field with a peculiar group of 4 "twin" galaxies was actually reported in a field near the quasar 
UM 679 [26]. 

These 4 pairs are remarkably twinlike with characteristic separations ~ 2.4". The sepa- 
ration of the twins (and the width of the strip over which the splitting occurs) is determined 
by /i and a value ~ 2.4" roughly corresponds to a string mass of 10"^ in dimensionless units. 
Further investigation of the candidate string-lensing field [27] revealed that in an area of ±1' 
around the original twins there are 7 twinlike galaxies, which satisfy the magnitude and color 
difference criteria for being lensing pairs. These twins vary wildly in magnitudes and colors 
but the distribution of separations is strongly peaked at 2" -2.5". 

One more interesting phenomena can occur if the background source lies partially out of 
the wake since the galaxies are comparable in sizes with the lensing strip. Whereas a fully 
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lensed galaxy should have the same color and spectrum in both images, in a partially lensed 
one the presence of strong color gradients, such as in a UV excess nucleus, can produce color 
differentiation between the images. Such an event was observed in the reported case. Here 
one member is bluer in continuum light, while the pair is identical in images of emission lines. 
Thefore, this object can be a partially lensed galaxy. 

Strings can exist also in loops, though the exact metric of a long-lived loop string is un- 
known today. Its lensing properties in the linearised gravitational approximation were studied 
by Gott [28] and Wu [29]. It was shown that the loop can produce three images if the source is 
inside the loop; one is the original source seen through the loop, and two images on opposite 
sides from the light rays which passed outside the loop. 

2.4.3 Detection of gravitational waves by lensing 

Gravitational lenses can be used to detect gravitational waves, as a gravitational wave affects 
the travel-time of a light ray. In a gravitational lens, this effect produces time delays between 
the different images. Such "detectors" are 22 orders of magnitude larger than any of the ex- 
isting or contemplated detectors, and they are sensitive to much lower frequences [30]. Time 
delay, produced by the wave, can be measured if the source (e.g., quasar) has variable bright- 
ness, thus producing images whith brightness variations correlated to a time-shift. For this 
purpose, the most useful systems would be those which are highly symmetric, so that the time 
delay due to the difference in the path length is small. 

A gravitational wave affects the time delay because it perturbes the metric tensor, and 
therefore, modifies the path lenght of the two light rays. The metric is given by 

gab = r]ab + hab (2.70) 

where rjab is the Minkowski metric, and hab is a small perturbation. One can calculate the 
time delay by examining the influence of the metric perturbation on the equation of the null 
geodesic. The measured time delay in the lens system can then be used to put an upper limit 
on the amplitudes of stochastic background of the gravitational waves at low frequencies. 
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Sensitivity of such a "detector" is greatest at wavelengths comparable to the overall size of 
the lens system. Allen [30] made such an estimation for the lens 0957 + 561 for the frequency 
range ijj>1ttx 10~^^ Hz. The amplitude of the gravitational waves must be less than 



/j < 2 X 10" 



2 X 10-17 



sec 

— T— (2.71) 
radian 



or the expected time delay A?„„ would exceed the measured value of 420 days. 

Though it may be difficult to separate the "intrinsic", geometrical time delay, from the 
delay caused by the gravitational wave, this idea may be still useful if the gravitational wave 
amplitude h is larger than 77, the angle between the images. The spatial motion of the geodesies 
that form the two images, induced by the gravitational wave, becomes significant. One effect of 
the gravitational perturbation is to change the angle between the images, usually by increasing 
it. The relative intensities of the two images also change. 



2.4.4 Determination of the lens parameters from gravitationally lensed 
gamma-ray bursts 

Recent results [31] from the Burst and Transient Source Experiment (BATSE) on the 
Compton Gamma Ray Observatory show an isotropic distrubution of gamma-ray bursts across 
the sky and rule out a population of sources within the Galaxy. The most natural explanation 
is that the bursts have a cosmological origin. If 7-ray bursts occur at high redshifts, then some 
bursts are likely to undergo gravitational lensing by foreground matter. This would lead to the 
detection of multiple bursts with identical profiles but with different time delays and magni- 
fications from a single event. Given a set of a multiple bursts (two or four), produced by a 
gravitational lens, what one one deduce about the nature of the lensing mass? Narayan and 
Wallington [32] showed that, if the lens is compact and pointlike, the quantity (1 +z\)M can be 
determined directly from the observations without any information about the angular diameter 
distance to the lens or the source, and without knowledge of the source redshift. Here M is 
mass of the lens and z\ is its redshift. What does a determination of {\+z\)M mean? First, it 
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gives an upper bound on the lens mass 




(2.72) 



where 



Pi = 



(l+Zi)4GM 

C3 



(2.73) 



This bound is independent of the size or geometry of the universe and of the redshift of the 
source or the lens. Second, if one can obtain an upper bound on zu then one will also have 
a lower bound on M. Using the V/Vmax data from BATSE, most cosmological models of 
gamma-ray bursts currently estimate the redshifts of the faintest observed sources to be Zs < 
1.5. Accepting this estimate, we have Zi < 1.5, and obtain the lower bound 



The ability to bound the mass of the lens from both sides to within a factor of ~ 2.5 is an im- 
pressive accomplishment. For a point mass lens the measurable parameter pi has a magnitude 
given by 



which means that we can hope to detect point lenses with masses > IO^Mq. If a significant 
fraction of the mass of the universe is in the form of compact objects with masses of up to ~ 
10^ Mq in the dark halos of galaxies or in the intergalactic medium, then 7-ray bursts will reveal 
their presence through lensing and will provide accurate mass determinations. Moreover, once 
a sufficient number of lensed bursts has been detected, statistical techniques may be used to 
determine the fractional mass density ili^ns of the lenses [33]. 

If the lens is not pointlike but has an extended mass distribution, the observations can be 
used to obtain the velocity dispersion of the lens. If galaxy lenses are not singular, but have 
finite cores, then the image configuration will consist of three or five images. The extra burst 
will have the longest time delay and will generally be significantly weaker than the rest. The 
relative magnification of this burst will give useful information on the core radius of the lens. 



Mmin = Mjnax/(1 +Zl,max) = . 



(2.74) 




(2.75) 
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2.4.5 Discovering planetary systems through gravitational microlensing 

Traditional methods for the search of planetary systems involve either indirect observations 
(e.g. astrometry measurements) or direct infrared observations (search for dust lanes around 
the stars). A new method, based on the gravitational lensing of the bulge star by an intervening 
disk star was suggested in 1992 [34]. Planetary systems of galactic disk stars can be detected 
through microlensing of stars in the Galactic bulge. Planets in a solar-like system located 
half-way to the Galactic Centre should leave a noticeable signature on the light curve of a 
gravitationally lensed bulge star. 

The gravitational lensing of distant sources by intervening individual compact masses, typ- 
ically in the mass range 10^ Mq >M> XQT^Mq, is called "microlensing". This term originates 
from the fact that the undetectable separation between images is of the order of microarcsec- 
onds for a solar mass located at a cosmological distance. For Galactic stars, however, the 
separation is of the order of a milliarcsecond. This angular separation is too small to observe. 
However, the resulting magnification can change the integrated light from the images for the 
time duration of the microlensing event. The brightness of the lensed star increases, peaks, and 
then decreases. The resulting light curve is smooth and completely described by three parame- 
ters: the temporal width, the maximum magnification and the time of maximum magnification. 
The duration of such an event is from several weeks to several months. It is symmetric in time 
about its maximum magnification and is achromatic, which allows one to distinguish it from 
variable stars. 

If there is a planet around the lensing star, the light curve may be significantly altered. The 
planet of mass m will typically affect the image magnification only for a fraction of {m/My/^ 
of the duration of the entire event. With the typical stellar velocities (~ 200 km/sec), this is a 
day or so for a Jupiter-mass planet. The observed light curve will look almost exactly like a 
light curve of an isolated star, except for a sharp spike during the fraction of time when a source 
moves inside a planetary Einstein radius. The planet affects appreciably the microlensed image 
only if the planet and the unperturbed image are separated by a distance of the order of the 
planet's own Einstein radius, (4GMZ)eff/c^)^/^, where Z)eff = DlDs/^ls- 

37 



The probability of detecting such events in the total number of microlensing events can 
be estimated. If a solar-like planetary system lay at a random position along the line of sight 
to the Galactic bulge, and if the bulge source came within one Einstein radius of the central 
star of this system, then the system could be detected ~ 20% of the time (assuming minimum 
detectable perturbation 6,„in = 5%). The largest contributor will be a Jupiter-like planet, ~ 17%, 
Saturn-like will give ~ 3%, and all the other planets <^ 1%. 

From the light curve the ratio of planetary to stellar mass, e = m/M, can be determined. 
If the lensing star is a G dwarf or earlier, its spectrum can be taken. From the spectral type 
and the luminosity one may determine the mass and distance and, thereby, infer the mass of 
the planet and its projected distance from the star. The typical planetary signal lasts for a day 
or less. Thus, to detect a Jupiter-mass planet, observations should be taken every 4 hours, 
detection of a Neptune-mass planet will require hourly observations. 

2.4.6 Light deflection in strong gravitational fields 

Up to now we were considering weak gravitational fields, where the essential assumption for 
lensing formalism is that the bending angles of light rays are very small. All tests of GTR 
within the Solar System, including the bending and delay of light rays passing the Sun, have 
examined the gravitational interaction only in connection with weakly self -gravitating objects 
(for example, the Sun has a surface gravitational potential GM/c^R fa 2 x 10"^). The mea- 
sured relativistic effects are but small perturbations to Newtonian mechanics and these tests 
say nothing about the strong field situation. There are, however, astrophysical systems where 
gravitational fields are strong and light bending leads to new interesting effects. Observing 
such systems and measurement of their parameters can yield tests of GTR to a greater preci- 
sion. 

GL Effects in Accreting Systems 

Relativistic "looks" of a neutron star General relativistic effects are quite substantial for 
neutron stars (NS) of radii smaller than about 2rs, where = 2GMns/c^ is the Schwarzschild 
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radius of NS. One expects these effects to play a major role in the interpretation of the spectrum 
and light curves of such stars. The characteristic quantity here is the ratio p = r^/r^, where 
is NS's radius. For stars with p<2 radiation emitted deep inside the strong gravitational field 
of such a star will be significantly modified as seen by the distant observer. Light ray, emitted 
on or near the surface will be redshifted and may be deflected by more than 45" [35]. 

It was shown that rays reach the observer with larger impact parameters than they would 
in flat space and, in addition, photons from the "back" half of the star can reach the observer. 
As a result, the part of the star which would be visible in flat space now appears larger and we 
may see parts of the star which would otherwise be hidden from view. While in flat spacetime 
exactly half of the surface of the star is visible, for p > 1.7 the whole star becomes visible, the 
point at 180° appears as the circular boundary of the disk. 

For the NS with emission from polar regions the visibility of the hot spots is markedly 
different from the flat space case. Ref. [35] shows the simulated picture of a pair of antipodal 
hot spots at 12 different phases of rotation with the angle between rotation axis and a hot spot 
to be 45°. 

Light curves from relativistic neutron stars Pechenick et al [36] have investigated the 
influence of GL effects on the beaming of radiation from a hot spot on the surface of a slowly 
rotating star in the limit where Schwarzschild metric is applicable. It was demonstrated that 
the deflection of light in the vicinity of a NS produces large deviations of model light curves 
from those expected in the absence of gravitational effects. For thermal emission, gravity 
tends to flatten the light curve for the NS with p>2, giving very little observed variation for 
nominal pulsar values. For NS with p < 2 a new feature appears in the form of a spike at 6* = 0°, 
whose width is about the angular diameter of a polar cap and whose relative height increases 
with decreasing of a polar cap radius. For the most relativistic case this feature represents a 
jump of l.lmag over the essentially flat continuum. The sharp peak is actually due to the 
gravitational lensing effect of the star on the cap at ^ = 180° to the line of sight. It arises from 
photons emitted near the tangential plane and bent through large angles (up to 227°). The 
region near 9 = 180° appears as a ring whose apparent brightness exceeds that of a comparable 
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region a.t6 = 0°. Thus, gravity alone is capable of producing strongly beamed radiation even 
for isotropically radiating polar caps. 

These results are applicable to the structure of X-ray pulsars in binary systems where ac- 
cretion is assumed to give rise to hot polar regions. That the X-ray light curve arises simply 
from rotational eclipsing of these caps seems unlikely. One interesting possibility is that emis- 
sion originates above the NS surface. It was demonstrated that for the NS with p = 2.5 any 
emitting region higher than 0.19i?NS, regardless of shape and size, will always be in view. 
For a typical NS no point that is more than 2 or 3 km above the surface is ever out of sight. 
For accreting pulsars it would be relevant for the case when the infalling matter decelerates at 
some distance above the surface due to shocks or radiation pressure. In this case, one has an 
accretion "column" rather than a polar cap or hot spot. In the general relativistic treatment of 
the emission from the column, there appear two new qualitatively different effects [37]. First, 
the frequency redshift is different for radiation arising from different heights. Second, the star 
and the accretion column will both produce some shadowing of the light rays. This effect 
depends on the emission height and on the direction of observation. In general, one sees that 
the column beam is strongly backwardly bent for the most relativistic cases. 

Emission from accretion disk The spectrum of X-rays produced by an accretion disk 
around a black hole is influenced markedly by GL effects [38]. Due to the forward "peak- 
ing" of the emitted radiation by the rapidly moving gas in the inner disk and to the gravita- 
tional focusing, the radiation is concentrated toward the equatorial plane. The concentration 
toward the plane from the inner disk is more severe for large values ofa/M, where a is angular 
momentum per unit mass and M is mass of the black hole. As the disk thickens in its outer 
regions, a distant observer directly in equatorial plane is in a shadow and in Newtonian case 
(disk accretion model neglecting the relativistic effects) receives no radiation. In actuality, 
he sees some blueshifted radiation from every radius outside the radius of marginal stability 
rms (radius at which gas begins to plunge into the hole). Radiation seen by other observers is 
less blueshifted, and that seen by the axial observer is always redshifted. For a Schwarzschild 
black hole, the effects of redshift and focusing are minor, since a disk around such black hole 
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has a large inner radius (rms = 6M). Only the equatorial observer sees a spectrum different 
from Newtonian spectrum (since this observer sees only focused radiation he always sees a 
spectrum dominated by high-energy radiation from the inner disk). For Kerr black hole with 
a < 0.9M redshift and focusing effects on the observed spectrum are striking. Even though 
all observers receive practically the same integrated flux, the average photon energy differs 
by approximately an order of magnitude between the equatorial and axial observer. The axial 
observer sees radiation from the cold, outer regions of the disk primarily; radiation from the 
inner regions is redshifted and defocused. Consequently, he sees a spectrum, attenuated at 
high energies in comparison to the Newtonian spectrum. The equatorial observer sees radia- 
tion from the hot, inner region to be blueshifted and strongly focused, consequently, spectrum 
is enhanced at high energies, compared with the Newtonian. If emission is not isotropic, for 
the axial observer «e ~ (rie — the angle of the emitted radiation with the surface normal) for 
radiation from large radii. This angle increases as decreases. For other observers, the emis- 
sion angle decreases as increases, reaching a minimum for = 3-4M. The radius for min 
is smallest for the equatorial observer. Thus, radiation reaching the observer from very small 
radii has large emission angles. Near the horizon radiation must be emitted parallel to the disk 
surface to escape being trapped by the disk or the hole. 

GL Effects in Compact Binaries 

Compact binaries present a unique laboratory for testing GR in the strong-field limit. This has 
become possible due to increasing discoveries of binaries, where one compact object is a neu- 
tron star (pulsar, beaming at us) and the other is a white dwarf (WD) (for example, 1855+09), 
a neutron star (NS) (for example, 1913+16 system) or a black hole (BH) (for example, a binary 
PSR B0042-73 in SMC is argued to have a massive 10 — 30 Mq BH companion [39]). 

To date, there exists one test of relativistic gravity [40], the a;- 7 -P^, test. It is obtained 
by combining the five timing parameters of the binary: eccentricity — e, orbital period — P and 
three "post-keplerian" (PK) ones: advance of periastron a), time dilation and gravitational shift 
parameter 7 and the orbital period decay P. All these are linked by the theory-dependent 
constraint, which can be defined as P/f^^°'y(e, P, 7, a)) = 1. For PSR 1913+16 Damour [40] 
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finds 

pobs 

^ = 0.995 ±0.021 (2.76) 

By contrast, in other theories of gravity the influence of strong-field effects on the function 
ytheory^^^ P, 7, cu) is enormous and it changes drastically both in sign and in magnitude [40]. 
Therefore, this result constitutes good (better than 0.5% accuracy [41]) confirmation of GR 
and, above all, a very selective confirmation of the ability of GR to describe the strong (and/or 
rapidly varying) gravitational fields. The uj-^-Pb test is a mixed test, which combines strong- 
field effects (related to a) and 7) with radiative effects (related to P). Thus, one cannot logically 
conclude, when it is satisfied, that both the specific strong-field and radiative predictions of 
GR have been independently confirmed [42]. In the case of a binary's orbit inclination to 
the sky plane approaching 90° the pulsar's two more PK parameters, characterizing the range 
r = Gnic/r' and the shape s = sini of the Shapiro time delay, may be measured. This has 
already allowed the determination of the neutron star mass with a precision better than 18% 
(B1855-h09 [43]). 

Self-lensing by binaries Eclipsing binaries present a remarkable situation where, if the or- 
bital plane inclination to the plane of the sky is close to 90°, the gravitational amplification of 
the companion by the compact object is possible. The observable effect here is periodic bright- 
ness enhancement which would depend on the geometry of the system. Schneider [44] gives 
for the binary pulsar PSR 1957-1-20 the periodic brightness amplification between 10^ and 
10"^, with a probability that amplification really exists at about 10%, if an adopted radius of 
optical companion is Rc = OARq. For an impact parameter b = 0, the duration of amplification 
is of Ta ~ 400 sec. Though the variability in brightness of an optical companion was observed 
in the case of PSR 1957-1-20, Schneider warns that it is still difficult to make the presently pre- 
dicted amplification observable. The star can have intrinsic variability due to the geometry of 
the "evaporation" mechanism driven by the pulsar. Gould [45] suggested monitoring millisec- 
ond binaries to search for lensing events. He showed that self-lensing in binary pulsars could 
be used to probe the structure of the emitting region of pulsars, in particular, to check whether 
the emission originates from the light cylinder or from a smaller locus. For pulsars of period P, 
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the light cylinder has a radius Pc/li: ^ 10^ [P/ (20 m*)] cm. If the lensing magnification curve 
follows the classic form for a point source, this would show that the source is small on scales 
of the Einstein ring, i?E ~ 3 x 10^ {a/W^f/'^ cm. 

Deflection of the light beam by companion If the inclination of the sky plane of the orbit 
of binary pulsar / is close to 90°, the pulsar beam during a fraction of the orbital revolution will 
be gravitationally deviated by the companion and this effect will manifest itself in the timing 
formulae. Bending of the pulsar beam by a companion can lead to the deviation of the apparent 
position, to brightness amplification of the source and to time delays. As in the Shapiro effect, 
the deflection of pulsar's beam shows in the timing as a rapid, sharp growth of the magnitude 
of the post-fit residuals of times-of-arrival (TOAs) on a short time interval during the superior 
conjunction of the pulsar and its companion. This effect is superimposed on the Shapiro effect 
[46]. 

Pulsar's TOAs are determined by measuring the phase offset between each observed pro- 
file and a long-term average one. In addition, a set of post-fit residuals are determined — the 
differences between measured TOAs and those calculated using the classical spin-down model 
of pulsar rotation. When many TOAs are available, spaced over months or years, it generally 
follows that, at least the pulsar's celestial coordinates, spin parameters and Keplerian orbital 
elements will be measurable with high precision. The binary systems most likely to yield 
measurable PK parameters are those with large masses, high eccentricities and astrophysically 
"clean", so that orbits are mostly dominated by gravitational interactions between two com- 
pact masses. Apart from mass determination, measurement of the light deviation effect in near 
edge-on compact binaries would help to get further constraints on the position of the pulsar's 
spin. In this effect the shape of TOAs residuals strongly depends on the spatial orientation of 
the pulsar's rotational axis [46]. 

Epstein [47] was, perhaps, the first who paid attention to the phenomenon of pulsar beam 
deviation in the gravitational field of it's companion in the case of cos / < 10~^. To be detected, 
the PSR beam must, after each revolution, point in the direction defined by the null geodesic 
connecting the orbital location of the pulsar to the observer. Due to the light bending produced 
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by the gravitational field, this direction depends on the orbital phase. The result of this grav- 
itational shift is a delay due to the excess angle swept out by the pulsar in each spin period. 
The influence of two relativistic effects (Shapiro and gravitational shift) should be compared 
with the special relativistic effect — abberation, which also causes shifts of the pulsar beam 
and delays. The discussion on the importance of the different delays appears in [48]. Author 
discussed the importance of the different delays which must be included in the timing formu- 
lae of binary system with orbit inclination close to 90° to the line of sight. He studies the 
gravitational shift delay for compact binaries at orbital phases close to superior conjunction. 
For PSR 0655-1-64, Schneider derives delays of 1 - lO/is (due to the gravitational shift effect) 
during a small fraction of the orbital period. Observation of this delay in a binary system 
must be considered as evidence that pulsars are indeed rotating beacons, instead of stars with 
periodic isotropic flares [49]. All derivations are done in standard Schwarzschild formalism, 
neglecting the changes imposed by the mass of the pulsar in the system, which bears some jus- 
tification. When the pulsar mass is assumed to be negligible, there exists only one geodesic, 
radial, joining the centres of the two stars. More than one such path would imply that a photon 
emitted radially from the centre of the companion star can aquire a transverse motion, which is 
not possible in the Schwarzschild solution. However, once the mass of the pulsar is "switched 
on", there exists a family of geodesies between the two stars, all of which have geometric path 
lengths greater or equal to the initial radial one. The addition of the pulsar mass must increase 
all geodesic path lengths and those passing near to both stars the most. Thus, it appears that 
the corrections due to the pulsar mass can only increase the magnitude of the predicted effects 
[50]. 

Doroshenko and Kopeikin [46] reduced the favourable inclination angle to cos/ < 0.003 
and showed that this estimate becomes less restrictive as the pulsar's spin axis approaches 
the line of sight. For the first time they showed the time behaviour of the residuals of the 
TOA in the vicinity to the moment of superior conjunction of the pulsar and its companion. 
Their numerical estimates showed that beam deviation effect is too small to be detected in the 
presently the best available example of PSR B 1855-1-09 with cos/ = 0.04. 
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Chapter 3 

Gravitational Lensing as a Tool in Search 
of Natural Wormholes and Negative 
Matter in the Universe. 

3.1 Introduction 

Energy conditions and their violations The energy conditions (EC) of classical GR for a 

case of FRW spacetime are a set of simple constraints on various linear combinations of the 
energy density and pressure. They are the null (NEC), the weak (WEC), the strong (SEC), and 
the dominant (DEC) energy conditions. For a FRW spacetime and a diagonal stress-energy 
tensor T^i, = (p,-p,-p,-p) with p the energy density and p the pressure of the fluid, they read: 

NEC <^ (p+P>0), 

WEC <^=^ (p > 0) and (p+ p > 0), 

SEC <^ (p+3/7>0)and(p+/7>0), 

DEC <S=^ (p>0)and(p±p>0). (3.1) 



The EC are, then, just simple constraints on various linear combinations of the energy density 
and the pressure of the matter generating the spacetime curvature. Since normal matter has 
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both positive energy density and positive pressure, it automatically satisfies the null, weak and 
strong energy conditions. If we impose the condition that speed of sound in normal matter is 
less than the speed of light, we would have in addition dp/ dp < 1. Assuming there is no cos- 
mological constant, integrating this gives p < p,so that dominant energy condition is satisfied 
as well. Roughly speaking, the violation of the DEC is usually associated with either a large 
and negative cosmological constant, or with superluminal acoustic modes. Quite recently, it 
was shown [51] that the presently favoured value of Hubble parameter [Hq = (65, 85)] implies 
that SEC must be violated somehere between the epoch of galaxy formation and the present. 
Since normal matter satisfies SEC, this would require the introduction of "abnormal matter" 
(we would need large quantities of abnormal matter, sufficient to overwhelm the gravitational 
effects of normal matter). Violating the SEC is usually associated with a positive cosmological 
constant or cosmological inflationary epoch. 

In the 60 's and 70 's these energy coonditions were widely used as key foundations for a 
number of important theorems. For example, for "the positive mass theorem" (which states that 
objects made of the matter that satisfies the DEC can never antigravitate), for a variety of theo- 
rems that predict the creation of singularities in stellar collapse and in cosmological scenarios. 
The discovery of Hawking that nonrotating black holes can evaporate and, correspondingly, 
their surface area can shrink (in direct violation of the area increase theorem — "second law 
of black hole mechanics", based on the SEC) forced physicists to face the fact that though 
in classical physics the EC are perfectly reasonable assumptions, quantum fields can violate 
them. Many physical systems, both theoretical and experimental, are known to violate one or 
more EC. Perhaps, the most quoted is the Casimir effect (for more information on this and 
other examples, see [12]). Typically, observed violations are of quantum nature and are of the 
order of h. It is currently far from clear whether there could be macroscopic violations of EC. 
If they do exist, macroscopic negative masses could be part of the ontology of the universe. 

Negative Masses The violations of the EC, in particular the weak one, would admit the 

existence of negative mass. The possible existence of negative gravitational masses has been 
investigated since the end of the nineteenth century [52]. From a Newtonian point of view, 
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we can differentiate four possible situations [53]: (i) both inertial and gravitational masses are 
positive, (ii) inertial mass is positive and gravitational mass is negative, (in) inertial mass is 
negative and the gravitational mass is positive, and (iv) both inertial and gravitational masses 
are negative. Most of the Nineteenth Century literature on negative masses is devoted to case 
(ii). From a relativistic point of view, however, the situation is quite different: the Equivalence 
Principle (EP)^ requires that gravitational and inertial masses cannot be considered distinct 
any longer. Test particles move along geodesic lines in accordance with the initial conditions, 
independently of the fact that their energy density is positive or negative. Negative energy 
densities or negative masses exert a repulsive force not only for ordinary matter, but for exotic 
matter as well. This might conceivably pose a stability problem for large amounts of exotic 
material, requiring large values of tension in order to keep an equilibrium configuration [54]. In 
more complicated astrophysical systems stability could be achieved through electromagnetic 
forces. 

Wormholes Of all the systems which would require violations of the EC in order to exist, 
wormholes are the most intriguing. The salient feature of these objects is that an embedding 
of one of their spacelike sections in Euclidean space displays two asymptotically flat regions 
joined by a throat. Wormholes basically represent bridges between otherwise separated regions 
of the spacetime (Fig. 3.1) and need a special kind of matter in order to exist. This matter, 
known as exotic, violates the EC. 

Wormhole solutions to the Einstein field equations have been extensively studied in the 
last decade (for the detailed description of the wormhole solution see Appendix A). A simple 
search in PRD online resources yields 57 papers in the last three years including the word 
"wormhole" in its full record. Msot of these are works deal with different trials to find analyt- 
ical solutions representing Lorentzian wormholes in space-time. There are known wormhole 
solutions representing rotating structures, charged, uncharged, static, evolving wormholes, 
wormholes within alternative gravity theories, etc. But, clearly, this is only one aspect of the 

'Definitions: Strong equivalence principle — Spacetime is everywhere a Lorentzian manifold. Freely falling 
particles follow geoesics of the metric. Weak equivalence principle — ^All freely falling particles follow the same 
trajectories independent of their internal composition. 
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Figure 3.1: Left: Embedding diagramme for a wormhole. Two mouths, joined by a tunnel, 
can connect regions otherwise separated or disconnected (here the normal space should fold 
as a sheet of paper, whereas the wormhole would be a tunnel from one side of the sheet to the 
other). Right: Embedding diagram for a black hole. The singularity here is represented as a 
pinch off of the wormhole tunnel. 

problem. Once one has an analytical solution, it is not at all clear that the physical constructs 
it represents can, or do, indeed, exist. This second aspect of the problem is not commonly 
treated in the literature: very little work has been done to develope astrophysical, or other kind 
of tests, to see whether wormholes really are part of the ontology of the universe (see [12, 55] 
and references cited therein. 

Can a wormhole have total negative mass? Since wormholes have to violate the null en- 
ergy condition in order to exist, the hypothesis underlying the positive mass theorem no longer 
applies. There is nothing in principle that can then prevent the occurrence of a negative to- 
tal mass. In other words, we need to have some negative mass near the throat to keep the 
wormholes open [56]. Examples of wormhole solutions, both with positive and negative mass, 
followed the pioneering work by Morris and Thome [54] (see [12] for a review). Visser, for 
instance, suggested a particular class of solutions, lacking spherical symmetry, configurations 
in which there is a flat-space wormhole, framed by struts of an exotic material [57]. These 
kinds of wormholes and many others, including the simplest one known as absurdly benign 
wormhole [54], would have a negative mass density. Typically near the throat the following 
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relationship would hold: 

radius • 

mass at throat = — . (3.2) 

For a radius equal to 1 meter the mass of the mouth will be equal to -1 Jupiter mass. The 
total mass seen at infinity^ will depend on the details of the model, such as the neighboring 
matter, and can be positive, negative or zero according to the specific case. Whether or not 
mass separation is possible is not clear yet. 

Some speculations have been made about how inflation can be responsible to enlarge a 
microscopic wormhole — believed to exist in the Planck foam — out to macroscopic dimensions 
[58]. This mechanism could result in a population of natural wormholes. 



Large-scale violations and observational strategy The result of the already mentioned 
work [51] on the violation of SEC in the recent epoch could imply the existence of a massive 
scalar field or a positive cosmological constant. These are favoured by current observations. 
Even if the global energy density of the universe is WEC-respecting, it is still not clear whether 
there exist spacetime domains where large-scale violations of the EC occur, allowing the for- 
mation of physical systems with an energy density equivalent to a total negative mass of the 
size of a galaxy or even a cluster of galaxies. 

Fifteen years after the seminal paper by Morris and Thorne [54], we face the following 
situation: there is no observational evidence supporting the existence of natural wormholes 
or serious theoretical reasons for their impossibility. Blackholes shared such a status during 
years until the discovery of galactic x-ray sources and quasars in the 1960s. Clearly, we have 
no better way than devising observational tests for deciding the existence of negative masses. 
For instance, if natural wormholes with negative masses or spacetime domains having large- 
scale violations of the EC actually exist in the universe (e.g. if the original topology after 
the Big-Bang was multiply connected), then there could be some observable electromagnetic 
signatures that might lead to their identification. The approach of the present work points out 

^In any asymptotically flat spacetime it is possible to define the so-called "imprint at infinity". Traversable 
wormholes, by definition, reside in asymptotically flat spacetimes. To deduce wormholes conservation laws, the 
ADM mass is invoked. It can be defined in terms of a suitable limit of surface integrals at spatial infinity. The 
most important feature of ADM mass is that it is conserved (see App. A). 
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in this direction: we do not know whether there exist an astrophysical population of objects 
with negative mass, but if they indeed exit, then we provide specific observational signatures 
that we can expect. 

Previous Studies on Lensing by Negative Masses The idea that wormholes can act as gravi- 
tational lenses and induce a microlensing signature on a background source was first suggested 
by Kim and Sung [59]. Unfortunately, their geometry was that of a perfect alignment of a 
source, both wormhole's mouths and an observer, which is quite unlikely. They also consid- 
ered both mouths to be of positive mass. Cramer et al. [60] carried out more detailed analysis 
of negative mass wormholes and considered the effect they can produce on background point 
sources at non-cosmological distances. The generalization to a cosmological scenario was car- 
ried out by Torres et al. [61], although lensing of point sources was still used. As far as we are 
aware, the first and only bound on the possible existence of negative masses, imposed using 
astrophysical databases, was given by Torres et al. [61]. Recently, Anchordoqui et al. [62] 
searched in existent gamma-ray bursts databases for signatures of wormhole microlensing. 
Although they detected some interesting candidates, no conclusive results could be obtained. 
Peculiarly asymmetric gamma-ray bursts [63], although highly uncommon, might be probably 
explained by more conventional hypothesis, like precessing fireballs (see, for instance, Ref. 
[64]). 

In the following sections we describe microlensing by natural wormholes of stellar and 
sub-stellar masses. We provide an in-depth study of the theoretical peculiarities that arise in 
negative mass microlensing, both for a point mass lens and source, and for extended source 
situations. We present negative mass microlensing simulations, showing the resulting shapes 
of the images, the intensity profiles, the time gain function, the radial and tangential magni- 
fications, and other features. Our work extends and deepens previous papers in several ways, 
and gives a method of analyzing observational predictions quantitatively. 

In the last section of this chapter we present the results of a set of simulations showing 
macrolensing effects we could observe if such a large amount of negative energy density ex- 
isted in our universe. 

50 



3.2 Negative mass lensing formalism and basic equations for 
a point lens 



In this chapter we consider lensing only by a point negative mass lens, and thus we can use all 
the assumptions concurrent with the treatment of the Schwarzschild lens (see Section 2.3.5). 
The image coordinates can be written as (6i , 62) and those of the source as (/9i , ^2)- 



3.2.1 Effective refractive index of the gravitational field of a negative 
mass and the deflection angle 

The 'Newtonian' potential of a negative point mass lens is given by 



where b is the impact parameter of the unperturbed light ray and z is the distance along the 
unperturbed light ray from the point of closest approach. We have used the term Newtonian 
in quotation marks since it is, in principle, different from the usual one. Here the potential 
is positively defined and approaches zero at infinity [12]. In view of the assumptions stated 
above, we can describe light propagation close to the lens in a locally Minkowskian spacetime 
perturbed by the positive gravitational potential of the lens to first post-Newtonian order. In this 
weak field limit, we describe the metric of the negative mass body in orthonormal coordinates 
x° = ct,x = (x') by 



where dl = \x\ denotes the Euclidean arc length. The effect of the space-time curvature on the 
propagation of light can be expressed in terms of an effective index of refraction rien given by 



*(^,z) = 



G\M\ 



(3.3) 



(^2 + ^2)1/2 ' 




(3.4) 




(3.5) 
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Thus, the effective speed of light in the field of a negative mass is 

2 

Vefr = c//leff ~ (3.6) 

c 

Because of the increase in the effective speed of light in the gravitational field of a negative 
mass, light rays would arrive faster than those following a similar path in vacuum. This leads 
to a very interesting effect when compared with the propagation of a light signal in the grav- 
itational field of a positive mass. In that case, light rays are delayed relative to propagation 
in vacuum — the well known Shapiro time delay. In the case of a negative mass lensing, this 
effect is replaced by a new one, which we shall call time gain. We will describe this effect in 
more detail in the following subsections. Defining the deflection angle as the differenceof the 
initial and final ray direction 

•— Cin^Cout) (3.7) 

where e := dx/dl is the unit tangent vector of a ray x(/), we obtain the deflection angle as the 
integral along the light path of the gradient of the gravitational potential 

cx = ^Jv±^dl, (3.8) 

where V_l$ denotes the projection of V$ onto the plane orthogonal to the direction e of the 
ray. We find 

Vxmz) = -^^^|^. (3.9) 
Eq. 3.8 then yields the deflection angle 

a=-^. (3.10) 

It is interesting to point out that in the case of the negative mass lensing, the term 'deflection' 
has its rightful meaning — the light is deflected away from the mass, unlike in the positive mass 
lensing, where it is bent towards the mass. 
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3.2.2 Lensing geometry and lens equation 




s 



o 



Figure 3.2: Lensing geometry of a negative mass. O is the observer, S is the source, W is 
the negative mass lens, 7i is one of the images. P is the angle between the source and the 
lens — position of the source, 6 is the angle between the source and the image — position of the 
image, and a is the deflection angle, b is the impact parameter and Z)/, and Di^ are angular 
diameter distances. Other quantities are auxiliary. 



In Fig. 3.2 we present the lensing geometry for a point-like negative mass. From this figure 
and the definition of the deflection angle (Eq. 3.10), we can obtain the relation between the 
positions of the source and the image. We only need to relate the radial distance of the source 
and the image from the center, since due to circular symmetry, the azimuthal angle is not 
affected by lensing. This gives 



(/3-6l)A=-aAs 



(3.11) 



or 




(3.12) 
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With the deflection angle (Eq. 3.10), we can write the lens equation as 



^ , 4G\M\Dis , 4G\M\ Di, 1 



3.2.3 Einstein radius and the formation of images 

A natural angular scale in this problem is given by the quantity 

4CM»^_ (3.14) 

which is called the Einstein angle. In the case of a positive point mass lens, this corresponds to 
the angle at which the Einstein ring is formed when the source, lens and observer are perfectly 
aligned. As we will see later in this section, this does not happen if the mass of the lens is 
negative. There are other differences as well. A typical angular separation of images is of 
order 29e for a positive mass lens. Sources which are closer than about 9e to the optical axis 
are significantly magnified, whereas sources which are located well outside the Einstein ring 
are magnified very little. All this is different with a negative mass lens, but nonetheless, the 
Einstein angle remains a useful scale for the description of the various regimes in the present 
case and, therefore, we shall use the same nomenclature for its definition. 

The Einstein angle corresponds to the Einstein radius in the linear scale (in the lens plane): 



= (3.15) 



In terms of Einstein angle the lens equation takes the form 



P = + ^, (3.16) 



which can be solved to obtain two solutions for the image position 9: 



9,,2 = \(^P±^Jp'-49ly (3.17) 
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Unlike in the lensing due to positive masses, we find that there are three distinct regimes here 
and, thus, can classify the lensing phenomenon as follows: 

1: /3 < 29^ There is no real solution for the lens equation. It means that there are no images 
when the source is inside twice the Einstein angle. 

II: P > 29e There are two solutions, corresponding to two images both on the same side of 
the lens and between the source and the lens. One is always inside the Einstein angle, 
the other is always outside it. 

ni: P = 29e This is a degenerate case, 9i^2 = 9e', two images merge at the Einstein angular 
radius, forming the radial arc (see Section 3.2.5). 
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Figure 3.3: Formation of images by a point negative mass lens. S is the source, Ii is the outer 
image and I2 is the inner image. Lens is in the centre, 9-^ is the angular scale of the lensing. 
This Figure shows schamaticall the radial distortion of images in the case of a negative mass 
lens. 

We also obtain two important scales, one is the Einstein angle (^e) — the angular radius of 
the radial critical curve, the other is twice the Einstein angle (26'e) — the angular radius of the 
caustic. Thus, we have two images, one is always inside the 9^, one is always outside; and as 
a source approaches the caustic (2 ^e) from the positive side, the two images come closer and 
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closer together, and nearer to critical curve, thereby brightening. When the source crosses the 
caustic, the two images merge on the critical curve (6'e) and disappear. In Fig. 3.3 we show 
schematically the formation and positions of images for the negative mass lensing. 

3.2.4 Time gain and time-offset function 

Following [20], we define a scalar potential, ip(6), which is the appropriately scaled projected 
Newtonian potential of the lens. 




(3.18) 



For a negative point mass lens it is 



As 4G|M 



In 1^1 . 



(3.19) 



The derivative of tp with respect to is the deflection angle 




(3.20) 



Thus, the deflection angle is the gradient of ip — the deflection potential. 



a(6>) = Veijj . 



(3.21) 



From this fact and from the lens equation (3.12) we obtain 



(6>-/3) + VeV^(6>) = 0. 



(3.22) 



This equation can be written as a gradient. 




(3.23) 
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If we compare this equation with that for the Fermat's principle [20] 



Vem = 0, (3.24) 

we see that we can define the time-offset function (opposite to time-delay function in the case 
of positive mass lens) as 



K^) = iL+^i)AA 



~ ^geom "I" ^pot • (3.25) 



Here tg^om is the geometrical time delay due to the extra path length of the deflected light ray 
relative to the unperturbed one. It remains the same as in the positive case — increase of light- 
travel-time relative to an unbent ray. The coefficient in front of the square brackets ensures 
that the quantity corresponds to the time offset as measured by the observer. The second term 

fpot is the time gain a ray experiences as it traverses the deflection potential ■'tp{6), with an extra 
factor (l+zi) for the cosmological 'redshifting'. Thus, cosmological geometrical time delay is 



W=^^i(«-«^ (3.26) 
and cosmological potential time gain is 



(l+zi)AA 

C As 

In Fig. 3.4 we show the time delay and time gain functions. The top panel shows tgeom for a 
slightly offset source. The curve is a parabola centered on the position of the source, /3. The 
central panel displays tpot for a point negative mass lens. This curve is centered on the lens. 
The bottom panel shows the total time-offset. Images are located at the stationary points of 
ttotai- Here we see two extrema — maximum and minimum — on the same side (right) from the 
optical axis (marked by dots). We can find the time difference between the two images, Q\ and 
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"total 



angular position 



Figure 3.4: Geometric time delay, gravitational time gain and total time offset produced by a 
point negative mass lens for a source that is slightly off the optical axis. 
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02. If a source has intrinsic variability, it will appear in the two images at an interval 



(3.28) 



where by u we denoted the ratio of absolute values of magnifications of images. 



1^2 



\/u^-A + u 



(3.29) 



and u is the scaled angular position of the source u = (3/9^ and is the Schwarzschild radius 
of the lens. 



3.2.5 Magnifications 



For a point mass lens the magnification is given by (see Section 2.3.4) 



H = 



From the lens equation (3. 16), we find 



(5 _e^+el _ d(5 _e^-el 

1~ 9^ ' 19' 9^ 



Thus, 



^1,2 



(3.30) 



(3.31) 



(3.32) 



Using definition for u in (3.29), we find the total magnification (Fig. 3.5, right panel, continu- 
ous curve) as 



A^tot = \Hl\ + \H2\ = 



u'-2 



The tangential and radial critical curves follow from the singularities in tangential 



(3.33) 



A'tan 



9^ + 91 



(3.34) 
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8 [arcsec] 



Figure 3.5: The magnifications: tangential ^tan (dotted lines), radial /irad (dash-dotted lines), 
and total /x (continuous curves), are plotted as functions of the image position 6 for two cases; 
in the left panel for the positive mass, in the right panel for the negative mass. The singularities 
of /itan and /Xrad givc the positions of the tangential and radial critical curves, respectively. In 
the left panel the singularity is in the tangential critical curve. In the right panel, instead, in 
the radial critical curve. Here \M\ - IMq, Ds - 0.05 Mpc and A = 0.01 Mpc. Angles are in 
arcseconds. 



and radial magnifications 



dl3 



d0 



(3.35) 



^rad diverges when Q = 6^ — angular radius of the radial critical curve, //tan always remains 
finite, which means that there are no tangential critical curves — no tangential arcs can be 
formed by the negative point mass lens. In Fig. 3.5 we show the magnification curves (radial, 
tangential and total) for both positive (left panel) and negative mass lenses (right panel). The 
difference can be seen as follows — in the left panel there is no singularity in the radial curve 
(no radial arcs are formed by the positive mass lens), whereas in the right panel we see that the 
curve for the radial magnification experiences a singularity. 
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3.3 Microlensing 



3.3.1 Light curves of the point source 



When the angular separation between the images 



(3.36) 



is of the order of milliarcseconds, we cannot resolve the two images with existing telescopes 
and we can only observe the lensing effect through their combined light intensity. This effect 
is called microlensing. Both the lens and the source are moving with respect to each other (as 
well as the observer). Thus, images change their position and brightness. Of particular interest 
are sudden changes in luminosity, which occur when a compact source crosses a critical curve. 
For a positive mass the situation is quite simple (Fig. 3.6). Positive mass microlensing is 
extensively reviewed in [23]). 

For a negative mass lens the situation is different. We define a dimensionless minimum 
impact parameter Bq in terms f the Einstein radius as the shortest distance between the path 
line of the source and the lens (all necessary definitions are illustrated in the Fig 3.7). 

For three different values of Bq we have three different lensing configurations shown in 
Figs. 3.8, 3.9 and 3.10. Note the large difference in the shapes of the images for these three 
regimes. Fig. 3.11 shows the case of a minimum impact parameter equal to zero, Bq = Q (the 
path of the source goes through the lens). 

It can be assumed without any loss of generality that the observer and the lens are motion- 
less and the source moves in the plane perpendicular to the line of sight (thereby, changing 
its position in the source plane). We adopt the treatment given in [65] for the velocity V, and 
consider effective transverse velocity of the source relative to the critical curve. We define the 
time scale of the microlensing event as the time it takes the source to move across the Einstein 
radius, projected onto the source plane, .^o = G^D^, 




(3.37) 
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Figure 3.6: Schematic representation of the geometry of the positive mass lensing due to the 
motion of the source, lens and the observer (in this case we can consider only the motion of 
the source in the plane perpendicular to the optical axis). The lens is indicated with a dot at 
the center of the Einstein ring, which is marked with a dashed line. The positions of the source 
center are shown with a series of small open circles. The locations and the shapes of the two 
images are shown with a series of dark ellipses. At any instant, the two images, the source and 
the lens are all on a single line, as shown in the figure for one particular instant. 

Ae = 2^4^\ ■ Bg 




1 regime 



2 regime 



Observer is 
above this plane. 



Time scale for the positive 
mass is justRg/V 



Figure 3.7: Schematic representation of the geometry of the microlensing. All teh quantities 
are defined in the text. 
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Figure 3.8: True motion of the source and apparent motion of the images for Bq>2. The inner 
dashed circle is the Einstein ring, the outer dashed circle is twice the Einstein ring. The rest is 
as in Fig. 3.6. 




Figure 3.9: True motion of the source and apparent motion of the images for Bq = 2. The inner 
dashed circle is the Einstein ring, the outer dashed circle is twice the Einstein ring. The rest is 
as in Fig. 3.6. 
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Figure 3.10: True motion of the source and apparent motion of the images for Bq < 2. The 
inner dashed circle is the Einstein ring, the outer dashed circle is twice the Einstein ring. The 
rest is as in Fig. 3.6. 




Figure 3.11: True motion of the source and apparent motion of the images for Bq = 0. The inner 
dashed circle is the Einstein ring, the outer dashed circle is twice the Einstein ring. Images 
here are shown with the shaded ellipses. The rest is as in Fig. 3.6. 
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Angle l3 changes with time according to 



P(0=\l[^J +Po- (3.38) 

Here the moment t = corresponds to the smallest angular distance /3o between the lens and 
the source. Normalizing (3.38) to 6e, 



"(0=t/(^V+f^V, (3.39) 



where the dimensionless impact parameter u is defined in (3.29). Including the time scale ty 
and defining 

Bo-^, (3.40) 

f^E 

we obtain 



u(t)=xlBl+(^pj . (3.41) 
Finally, the total amplification as a function of time is given by 

u(tf-2 

A(t)= ^ (3.42) 

Comparing this analysis with that of Cramer et al. [60], we note that they wrote the equation 
for the time dependent dimensionless impact parameter as (cf. our Eq. 3.41) 



B(t)^Bo\ll+(^] , 



and defined the time scale for the microlensing event as the time it takes to cross the minimum 
impact parameter (cf. our Eq. 3.37) 

where bo is the minimum impact parameter and other variables carry the same meaning as in 
our case. While there is no mistake in using such definitions, there is a definite disadvantage 
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in doing so. Using Eq. 10 of [60] for B(t) we cannot build the light curve for the case of the 
minimum impact parameter Bq - 0. In this case their Eq. 8 diverges, although there is nothing 
wrong with this value of Eq (see our Figs. 3.11 and 3.12). In the same way, their definition of a 
time scale does not give much information on the light curves. With our definition (Eq. 3.37) 
we can see in Fig. 3.12 that in the extreme case of Bq = the separation between the half-events 
is exactly 26'e; it is always less than that with any other value of Bq. 

In Fig. 3.12 we show the light curves for the point source for four source trajectories with 
different minimum impact parameters Bq. As can be seen from the light curves, when the 
distance from the point mass to the source trajectory is larger than 26^, the light curve is 
identical to that of a positive mass lens light curve. However, when the distance is less than 
26'e (or in other terms, Bq < 2.0), the light curve shows significant differences. Such events are 
characterized by the asymmetrical light curves, which occur when a compact source crosses 
a critical curve. A very interesting, eclipse-like, phenomenon occurs here. We have a zero 
intensity region (disappearance of images) with an angular radius 9o 



In the next subsection we shall see how these features get affected by the presence of an 
extended source. 

3.3.2 Extended sources 

In the previous subsections we considered magnifications and light curves for point sources. 
However, sources are extended, and although their size may be small compared to the relevant 
length scales of a lensing event, this extension definitely has an impact on the light curves 
as will be demonstrated below. From variability arguments, the optical and X-ray continuum 
emitting regions of quasars are assumed to be much less than 1 pc [66], whereas the broad-line 
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Figure 3.12: Light curves for the negative mass lensing of a point source. From the center of 
the graph towards the corners the curves correspond to Bq = 2.5, 2.0, 1.5, 0.0. The time scale 
here is divided by the effective transverse velocity of the source. 



emission probably has a radius as small as 0. 1 pc [67]. The high energy gamma-spheres have a 
typical radius of 10^^ cm [68]. Hence, one has to consider a fairly broad range of source sizes. 

We define the dimensionless source radius, R, as 

R = i = r^ (3.45) 

where p and R are the angular and the linear physical size of the source, respectively, and is 
the length unit in the source plane (see Eq. 3.37). 



3.3.3 Comments on numerical method and simulations 

It is convenient to write the lens equation in the scaled scalar form 

y^x+-, (3.46) 

X 
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where we nomialize the coordinates to the Einstein angle :^ 

x=f ; y=^. (3.47) 

The lens equation can be solved analytically for any source position. The amplification factor, 
and thus the total amplification, can be readily calculated for point sources. However, as 
we are interested in extended sources, this amplification has to be integrated over the source 
(Eq. 3.48). Furthermore, as we want to build the light curves, the total amplification for an 
extended source has to be calculated for many source positions. The amplification A of an 
extended source with surface brightness profile /(y) is given by 

where Ao{y) is the amplification of a point source at position y. We have used the numerical 
method first described in [69]. We cover the lens plane with a uniform grid. Each pixel on this 
grid is mapped, using Eq. 3.46, into the source plane. The step width (5000 x 5000) is chosen 
according to the desired accuracy (i.e. the observable brightness). For a given source position 
Cyio^yio) we calculate the squared deviation function (SDF) 

= (yw-yi(xuX2)f + (y2o-y2(xuX2)f . (3.49) 

The solutions of the lens equation (Eq. 3.47) are given by the zeroes of the SDF. Besides, 
Eq. 3.49 describes circles with radii S around CvicJio) in the source plane. Thus, the lines 
S = constant are just the image shapes of a source with radius S, which we plot using standard 
plotting software. Therefore, image points where SDF has value correspond to those points 
of the circular source which are at a distance S from the center. The surface brightness is 
preserved along the ray and if I(Ro) - h for the source, then the same intensity is given to 
those pixels where SDF= R^. In this way an intensity profile is created in the image plane and 



^Note, that for the case in which x and y are expressed in length units, we obtain a different normalization in 
each plane, which is not always convenient. 
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integrating over it gives the total intensity of an image. Thus, we obtain the approximate value 
of the total magnification by estimating the total intensity of all the images and dividing it by 
that of the unlensed source, according to the corresponding brightness profile of the source. 
For a source with constant surface brightness the luminosity of the images is proportional to 
the area enclosed by the line S = constant. The total magnification is obtained by estimating 
the total area of all images and dividing it by that of the unlensed source. For calculations of 
light curves we used a circular source which is displaced along a straight line in the source 
plane with steps equal to 0.01 of the Einstein angular radius. 

3.3.4 Shapes of images and light curves for the uniform brightness 
source 

For a circular source of radius R and uniform brightness, equation (49) transforms into 

Figs. 3.13 and 3.14 show four projected source and image positions, critical curves/caustics 
in the lens/source plane and representative light curves for different normalized source sizes. 
The sources are taken to be circular disks with constant surface brightness. In order to get 
absolute source radii and real light curves we need the value of Oe, the normalized angular 
unit, the distance to the source, as well as the velocity V of the source relative to the critical 
curves in the source plane. We have used M = Mq, //q = 100 km s"^ Mpc"^ and a standard flat 
cosmological model with zero cosmological constant. Here and in all subsequent simulations 
the redshift of the source is Zs = 0.5 and the redshift of the lens is zi = 0. 1. 

We display two cases for two different impact parameters. It must be noted here that the 
minimum impact parameter Bq defines now the shortest distance between the line of path of 
the center of the source and the lens. For each Bq, the dimensionless radius of the source R in- 
creases from 0.01 to 2.0 in normalized units of Oe- The shape of produced images changes no- 
tably with the increase of the source size, as can be seen in the bottom right panel of Figs. 3.13 
and 3.14. At the same time the smaller the source the greater the magnification, since when 
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Figure 3.13: Four sets of lens-source configurations (upper panels) and corresponding ampli- 
fication as a function of source's center position (bottom panels) are shown for four different 
values of the dimensionless source radius R (0.01, 0.5, 1.0, 2.0, in normalized units, 9e). Each 
of the four upper panels display the time dependent position of the source's center, the shapes 
of images (shaded ellipses) and critical curves (dashed circles). The series of open small circles 
show the path of the source center. The lens is marked by the central cross. Minimum impact 
parameter Bq = 0. By replacing ^ with ^ObD^V'^ - ^ty we get corresponding time depending 
light curve. 
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Figure 3.14: Same as in Fig. 3.13, but for minimum impact parameter Bq = 2.0. 
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the source radius is greater than the Einstein radius of the lens, the exterior parts, which are 
amplified, compete with the interior ones, which are demagnified. 

It can be noted that despite the noise in some of the simulated light curves the sharp peaks 
which occur when the source is crossing the critical line are well defined even for the smallest 
source. Note that all infinities are replaced now by finite amplifications, and that the curves are 
softened; all these effects being generated by the finite size of the source. Indeed, while the 
impressive drop to zero in the light curve is maintained, the divergence to infinity, that happens 
for a point source, is very much reduced. Note, that in cases of a large size of the source, the 
magnification is very small. If we would like to see bigger enhancement than that, we should 
consider sources of smaller sizes, approaching the point source situation (cf. Fig. 3.13, upper 
left plot). It is also interesting to note here that for the impact parameter 5o = 2.0 the light 
curve of a small extended source, though approaching the point source pattern (Fig. 3.12), still 
differs considerably from it (Fig. 3.14, upper left plot). 

3.3.5 Shapes of images and light curves for the sources with non-uniform 
brightness profiles 

In order to compare a constant surface brightness source with more realisitc distributions, we 
simulate image configurations and calculate light curves for two different assumed profiles 
with radial symmetry. 

Source with Gaussian brightness distribution 

2/2 

For a Gaussian source we have I(r) = I^e''' '''o, where we normalized the profile such that the 
maximum value of / equals unity. We define the radius containing 90% of all the luminosity 
as the effective radius of the source, R^. To find the relation between and vq, we write the 
total luminosity as 




(3.51) 




R 



(3.52) 
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then 



= 0.9= l-g-^iAo 



(3.53) 



L(oo) 



Thus, effective radius relates to the parameter tq as 



= VinTo. 



(3.54) 



In Figs. 3.15 and 3.16 we show the images of an extended source with a Gaussian bright- 
ness distribution for two effective dimensionless source radii ^s, 3.0 and 1.0, in units of the 
Einstein angle (Fig. 3.15, frames a-e and Fig. 3.16, frames a-e, respectively), together with 
the corresponding light curves (Fig. 3.15, frame f and Fig. 3.16, frame f, respectively). Here 
the source path passes through the lens (Bq = 0), which lies exactly in the center of each frame. 
In Fig. 3.15 the source's extent in the lens plane is greater than the Einstein radius of the lens. 
We notice there that there is an eclipse-like phenomenon, occurring most notably when most 
of the source is near or exactly behind the lens. This is consistent with the light curve (frame f), 
where there is a de-magnification. For the source with radius smaller than the double Einstein 
radius of the lens (Fig. 3.16), the low intensity region is replaced by the zero intensity region; 
the source completely disappears from the view (frame c). 

Source with exponential brightness distribution 

We have for this brightness distribution 7(r) = I^e'''/^'^. In the same manner as above, Rs is de- 
fined as radius, containing 90% of total luminosity. In the same way as above, total luminosity 




(3.55) 



then 




(3.56) 



and 



= 0.9. 



(3.57) 



L(oo) 
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Figure 3.15: Image configurations (frames a to e) and a corresponding light curve (f) for a 
Gaussian source with effective radius Rs = 3.0, in units of the Einstein angle. The source 
is moving from the lower left corner (frame a) to the right upper corner (frame e), passing 
through the lens (Bo = 0). The lens is in the center of each frame. Size of each frame is 5 x 5, 
in the normalized units. Wedges to each frame provide the brightness scale for the images. 
Note the eclipse-like phenomenon, consistent with the incomplete demagnification showed in 
the light curve (frame f). 

74 



^ 1 1 T' 1 1 1 1 1 ^ 1 1 1 1 1 1 1 

^ > I I . I . I ^ , I . I : I . I 

Q 42. 04 09 44 1 n Hri- n.n n« 1 






Figure 3.16: Image configurations (frames a to e) and a corresponding light curve (f) for a 
Gaussian source with effective radius Rs = 1.0, in units of the Einstein angle. The source 
is moving from the lower left comer (frame a) to the right upper comer (frame e), passing 
through the lens (Bq = 0). The lens is in the center of each frame. Size of each frame is 3 x 3, 
in the normalized units. Wedges to each frame provide the brightness scale for the images. 
Note the complete disappearence of the source when it is inside the double Einstein radius of 
the lens (frame c), corresponding to the drop of magnification to zero in the light curve (frame 
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Figure 3.17: Image configurations (frames a to e) and a corresponding light curve (f) for a 
source with exponential brightness distribution and effective radius Rs = 1.0, in units of the 
Einstein angle. The source is moving from the upper left comer (frame a) to the upper right 
corner (frame e) with the impact parameter Bq = 2.0. The lens is in the center of each frame. 
Size of each frame is 2.5 x 2.5, in the normalized units. Wedges to each frame provide the 
brightness scale for the images. 
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From where we find that effective radius relates to the parameter ro as 

e-«sAo ^1 + :^^ =0.1 . (3.58) 

Solution to this equation gives Rs/ro ~ 3.89. This profile is also normalized such that the 
maximum value of I equals unity. 

In Fig. 3.17 we display the images of the source with the exponential brightness profile and 
the corresponding light curve (frame f). The effective radius of the source is 1.5. The impact 
parameter here is 5o = 2.0; the lensing regime corresponds to the one schematically depicted 
in Fig. 3.9. We see how shapes of the images change, becoming elongated and forming the 
radial arc (frames c and d). 

Comparing light curves for sources with different brightness profiles 

In Fig. 3.18 we compare light curves for three different radially symmetric source profiles: 
uniform, Gaussian and exponential, for two dimensionless source radii R = = iJg"^"" = 0.1 
and R - Rf"^^ = iJ^^^™ = 1 .0. As a reference curve we show the light curve of the point source. 
All curves are made for the impact parameter Bq = 0. We can see larger noise in the uniform 
source curve, since the source with uniform brightness has an extremely sharp edge, whereas 
Gaussian and exponential sources are extremely smooth. Though we considered the sources 
with the same effective radius, we can see from the plot that for a small source size, the 
maximum magnification is reached by the source with exponential profile (upper panel). This 
is explained by the fact that this profile has a narrower central peak than the Gaussian. 

For the larger source, this behaviour smoothens, though we still can see large differences in 
the light curves (bottom panel). Here the uniform source experiences darkening, while sources 
with other profiles only undergo demagnification. 

3.3.6 Time scales of microlensing 

Let the source have a transverse velocity Vs measured in the source plane, the lens a transverse 
velocity Vi measured in the lens plane, and the observer a transverse velocity Vobs measured in 
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Figure 3.18: Light curves for the point mass source (dashed line), source with constant surface 
brightness (solid), source with Gaussian brightness distribution (dash-dotted) and exponen- 
tial brightness distribution (dotted) for two different effective dimensionless source radii, 0. 1 
(upper panel) and 1.0 (bottom panel). 

the observer plane. The effective transverse velocity of the source relative to the critical curves 
with time measured by the observer is 

1 1 Ds 1 Du 

V= -— vs--— -^vi + -— -i^vobs . (3.59) 
1+Zs 1+ziA 1+ziA 

This effective velocity is such that for a stationary observer and lens, the position of the source 
will change in time according to = \At. 

We basically have two time scales of interest here. The first one is the typical rise time to 
a peak in the amplification. We can estimate that it corresponds to a displacement of Ay R 

78 



Table 1 : Time scales for several source radii. R is the dimensionless source radius, in units 

of Einstein angle, \M\ = I.OMq, redshift of the lens is zi = 0.1, redshift of the source is Zs = 
0.5, = 5.42 X 10^^ km is the normalized length unit in the source plane. The time scales 
correspond to apparent source velocity (see Eq. 3.59) V = 5000 km s"^ 



R 


R (pc) 


A/|2(sec)" 


Tl (yr) 


T-2(yr)" 


0.0 


point source 


2.0 X 10-4 




6.78 


0.01 


1.07 X 10-4 




0.03 


6.81 


0.1 


1.07 X 10-3 




0.34 


7.11 


1.0 


1.07 X 10-2 




3.38 


10.16 


2.0 


0.3 X 10-1 




6.75 


13.6 



" u=4.0 (definition in Eq. 3.29) 
* Bo = 0.0 

of the source across a critical line. The corresponding time scale is ri = t^R, with ty given by 
(3.37). In terms of the physical source size R = ^o^> 

n = ^ , (3.60) 

where effective transverse velocity of the source V is given by Eq. 3.59. The second time scale 
of interest is the time between two peaks T2. For a point source we can estimate it as T2 = tyA, 
where A is given in Eq. (3.44). For a source with radius R and impact parameter Bq = 0.0 it 
can shown to be 

r2 = ?v(A+^) . (3.61) 

In Table 1 we list the time scales ri and T2 for different values of a source radius, and the 
time delay between two images Atn for the point source. The value of V is estimated to be 
y = 5000kms-i. 



3.3.7 Concluding remarks 

In these two sections we have explored the consequences of gravitational microlensing follow- 
ing from the existence of matter violating the energy conditions. We have also quantitatively 
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analyzed, using numerical simulations, the influence of a finite size of the source on the grav- 
itational lensing negative-mass event. We have thus enhanced and completed previous works 
where the focus was put on the point source light curves and no discussion was given concern- 
ing the shapes of images, actual simulations of microlensing events, time gain function, and 
other features presented here. Figs. 3.5, 3.8-3.11, 3.13-3.18 and Table 1 comprise our new 
results: a useful comparison arena where possible existence of wormholes or any other kind 
of negative mass compact objects can be observationally tested. 

The next step would be to test these predictions using archival, current, and forthcom- 
ing observational microlensing experiments. The only search done up to now included the 
BATSE database of 7-ray bursts. There is still much unexplored territory in the gravitational 
microlensing archives. We suggest adaption of alert systems of these experiments in order to 
include the possible effects of negative masses as well. This would lead to a whole new world 
of discoveries. 

3.4 Macrolensing 

In this section we present a set of simulations showing macrolensing effects we could observe 
if a large amount of negative energy density exists in our universe. The physical systems could 
have an energy density equivalent to a total negative mass of the size of a galaxy or even a 
cluster of galaxies. In what follows, we used a background cosmology described by a FRW 
flat universe with fim = 1 and a zero cosmological constant. In all numerical computations 
a dimensionless Hubble parameter h is put equal to 1 (Hq = lOOh km sec"^ Mpc"^). The 
relationship Dgff = Z)lZ)ls^s' is a measure of the lensing efficiency of a given mass distribution. 
Z)eff peaks, quite independently of the cosmological model assumed, at a lens redshift of ~ 
0.2-0.4 for sources at a typical redshift Zg ~ 1-1.3 [70]. To be more realistic, we placed the 
lens at a redshift of zi = 0.3 and generated a random sample of galaxies in the redshift range 
0.3 < Zsource < 2.0. The redshift distribution conserves their comoving number density. This 
number density and the projected sizes for these background sources were taken to be close 
to the Tyson population of faint blue galaxies [71]. The luminous area of each galaxy was 
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taken to be a circular disk of radius R with a uniform brightness profile and orientations of 
disk galaxies randomly placed in space. This was done by defining in the code the ellipticity e 
(e = (l-r)/(l + r), r being the ratio of the minor axis to the major axis) and the position angle 
(f, and randomly chosing the values of e from the range < e < 0.7 and the values of (p from 
the range < ip <2n.^ 

We write the gravitational lens equation, which governs the mapping from the lens to the 
source plane, in dimensionless form 



(3.62) 



The lens equation (3.62) describes a mapping 9^/3, from the lens to the source plane. For 
convenience, we redefine the lens plane as x and the source plane as y. Then, Eq. (3.62) can 
be written as 



y = x|l + ^), (3.63) 



where x = |x| = y^cf+^. 

We now consider a source, whose shape — either circular or elUptical — can be described by 
a function x(y). Curves of constant x are the contours of the source. One can as well consider 
X as a function of x, x(y(x)), where y(x) is found using the lens equation. Thus, all points 
X of constant x mapped onto points y, which have a distance from the centre of the 
source. If the latter contour can be considered an isophote of a source, one has thus found the 
corresponding isophotes of the images (see Section 3.3.3 for details on the algorithm). 



3.4.1 Simulations results 

In the first set of figures (Figs. 3. 19-3.23), we show the results of our simulations. Some special 
precautions must be taken for the largest masses. The problem is that for a very massive lens 
the Einstein ring becomes very large. Since for the negative mass lensing all sources inside the 
double Einstein radius are shadowed (i.e. we can see the images of only those sources which 

"^The random number generator needed in the code was taken from the book by Press et al. [72], and we use 
the algorithm described in Ref. [73], p.298. PGPLOT routine PGGRAY was employed in the code. 

81 



are outside the double Einstein radius), if we were to use for lensing only the sources shown 
on the window, a fold-four symmetry pattern would appear (Fig. 3.24). Only the sources at 
the comers of the current window (and outside the double Einstein radius) are lensed. In order 
to solve this problem we have to consider also the sources from outside the current window; 

then the lensing picture is restored and the scale of the simulation is consistently increased. 
For this reason we increase the number of background galaxies in Figs. 3.21-3.23. We show 
this in detail in Figs. 3.24 and 3.25. 




Figure 3.19: Left: Appearance of a background field of sources in a range of redshifts (see 
text) (200 galaxies, intrinsic radius 7 Kpc), when lensed by a negative mass of |M|iens = 1 x 
IO^^Mq. Center. Unlensed background field. Right: Appearance of the same background 
field of galaxies when it is lensed by an equal amount of positive mass; redhsifts are the same 
as for the negative mass case. 
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Figure 3.20: Left: Appearance of a background field of sources (200 galaxies, intrinsic ra- 
dius 10 Kpc), when it is lensed by a negative mass of |M|iens = 1 x 10^"^ M©. Center: As in 
Figure 3.19. Right: As in Figure 3.19. 
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Figure 3.21: Le//: Appearance of a background field of sources (300 galaxies, intrinsic radius 
15 Kpc), when it is lensed by a negative mass of |M|ie„s = 1 x 10^^ Mq. The simulation was 
made taking into account sources located within 1.2 of the size of the shown window. Center. 
As in Fig. 3.19. Right: As in Fig 3.19. 



As a general feature of our simulations we can remark that, contrary to the standard positive 
mass case, where ring-like structures appear, the negative mass lensing produces finger-like, 
apparently "runaway" structures, which seem to escape from a central void. This is in agree- 
ment with the appearance of a central umbra in the case of a point-like negative mass lensing 
(see previous sections of this chapter). In the case of macrolensing, the umbra (central void) 
is maintained on a larger scale which, depending on the negative mass of the lens, can reach 
hundreds of acrsec in linear size. This umbra is always larger than the corresponding one gen- 
erated in positive macrolensing (see Figures 3.19-3.23) and totally different in nature [60]. At 
least qualitatively, the existence of a macroscopic amount of negative mass lens mimics the 
appearance of galaxy voids. 

In order to explore the influence of the adopted redshift values, we turn now to the case 
where Zsources = 0.08 and ziens = 0.05. The Bootes void is the closest void to us, and lies between 
the supercluster Corona Borealis (z ~ 0.08) and Hercules (z ~ 0.03) [74]. This serves as 
motivation for the selection of these redshift values. As an example of the results for different 
lens masses, we show in Figs. 3.26 and 3.27 the cases with |M|iens = 1 x 10^"^ M© and |M|iens = 
1 X IO^^Mq. 

It is interesting to note that since the background population of galaxies is very dense, one 
would expect a lot of lensing in the standard model of cosmology. However, there is still a 
surprising dearth of candidates for (positive mass) lensed sources [75]. Some of the richest 
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Figure 3.22: Left: Appearance of a background field of sources (300 galaxies, intrinsic radius 
50 Kpc), when it is lensed by a negative mass of |M|iens = 1 x IO^^Mq. The simulation was 
made taking into account sources located within 1.2 of the size of the shown window. Center. 
As in Fig. 3.19. Right: As in Fig. 3.19. 
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Figure 3.23: Left: Appearance of a background field of sources (500 galaxies), intrinsic radius 
60 Kpc, when it is lensed by a negative mass of |M|iens = 1 x 1O^^M0. The simulation was 
actually made taking into account sources located within 1.5 of the size of the shown window. 
Center: As in Fig. 3.19. Right: As in Fig. 3.19. 
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Figure 3.24: The need for the increase of the background number of sources is shown by 
the appearance of a four-fold symmetric pattern, which occurs because only the galaxies at the 
corners of the left window are being affected by lensing effects. Axis are marked in arcseconds. 
See next figure. 
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Figure 3.25: Left: Appearance of a background field of circular sources (300 galaxies), each 
of them of 15 Kpc radius, when it is lensed by a negative mass of |M|iens = 1 x lO^^M© with 
an Einstein angular radius equal to 6e - 84 arcsec. The dashed circle is the Einstein ring, 
the radial arcs are centered on it while the images inside it are demagnified. Here ziens = 0.4, 
^source = 1-4. The simulation was made taking into account sources located within 1.2 of the 
size of the shown window. Center. Background field in the absence of the lens, dashed circle 
is the double Einstein radius, all sources inside this radius are shadowed. Right: Macrolensing 
effects produced by a positive mass lens of 10^^ Mq, a dashed circle is the Einstein radius 
shown here for comparison, sources inside it are strongly lensed. 
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Figure 3.26: Left: Appearance of a background field of sources (200 galaxies), each of them 
5 Kpc radius, when it is lensed by a negative mass of |M|iens = 1 x 10'"* Mq with an Einstein 
angular radius equal to 6-^ = 47 arcsec. Center. Unlensed background field. Right: Appearance 
of the same background field of galaxies when lensed by an equal amount of positive mass, 
located at the same redshift (see text). 
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Figure 3.27: Left: Appearance of a background field of sources (300 galaxies), each of them of 
25 Kpc radius, when it is lensed by a negative mass of |M|iens = 1 x lO^^M© with an Einstein 
angular radius equal to 9^ = 467 arcsec. Center: As in Figure 3.26. Right: As in Figure 3.27. 
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clusters do not display arcs in the deepest images. CL0016+16 (z = 0.56), for instance, is one 
of the richest and strongest X-ray emitting clusters. It is rather extended on the sky and the 
light from many background sources should cross this cluster. Neither arcs nor arclets have 
been found, though weak lensing has been reported. This may be pointing towards a cautionary 
note: if the kind of finger-like structures displayed in our figures is not directly seen in its full 
pattern, that does not necessarily mean that they are absent. Even the presence of one radial 
arc (without tangential counter arc and/or tangential arcs) may be significant. 

3.4.2 Concluding remarks 

The null EC (NEC) is the weakest of the EC. Usually, it was considered that all reasonable 
forms of matter should at least satisfy the NEC. However, even the NEC and its averaged ver- 
sion (ANEC) are violated by quantum effects and semi-classical quantum gravity (quantized 
matter fields in a classical gravitational background). Moreover, it has recently been shown 
that there are also large classical violations of the energy conditions [76]. Here we have shown 
that disregarding the fundamental mechanism by which the EC are violated (e.g. fundamental 
scalar fields, modified gravitational theories, etc.), if large localized violations of NEC exist 
in our universe, we shall be able to detect them through cosmological macrolensing. Contrary 
to the usual case, where arc-like structures are expected, finger-like "runaway" filaments and 
a central void appear. In Figures 3.19-3.23 and 3.26-3.27 we compare the effects of negative 
masses with the case of macrolensing effects on background fields produced by equal amounts 
of positive mass located at the same redshift. The differences are obvious. These results make 
the cosmological macrolensing produced by matter violating the null energy condition obser- 
vationally distinguishable from the standard situation. Whether large-scale violations of the 
EC, resulting in space-time regions with average negative energy density indeed exist in the 
universe can therefore be decided through observations. 
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Chapter 4 

Gravitational Lensing as a Tool to Test 
Alternative Cosmologies 

4.1 Introduction 

Standard cold dark matter FRW cosmology presents serious theoretical and observational dif- 
ficulties as an acceptable description of the Universe. For example, a study conducted on a 
sample of 256 ultra-compact sources [77] showed that the canonical CDM model (Qq = 

= 0) is ruled out at the 98.5% confidence level. An overview in the literature shows the 
existence of a growing body of work discussing alternative cosmologies [78]. The first moti- 
vation comes from the conflict between the age of the Universe and the age of the oldest stars 
in Galactic globular clusters. The ages of the globular clusters typically fall in the interval 
12 - 14 Gyr ([79], and references within), while measurements of the Hubble parameter were 
recently updated toh = 0.7 ±0.1 [80] with the value for the age in the interval 8.1 Gyr <to< 
10.8 Gyr. "The age problem" is even more acute if we consider its variant based on the age 
constraints from old galaxies at high redshifts [81]. 

Another important motivation is the cosmological constant problem. Historically, it was 
Einstein who introduced the A term in his field equations in order to make them compatible 
with a static universe. The A terms was introduced several times in cosmology but was dis- 
carded when improved data became available. Recently, observations by Perlmutter et al [82] 
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and Riess et al [83] of more than 50 type la supemovae suggested Friedmann models with 
negative pressure matter, such as a cosmological constant, domain walls or cosmic strings 
[25, 84]. The main conclusion of these works was that the expansion of the universe is accel- 
erating. Though a cosmic acceleration can also be accounted for by invoking inhomogeneity 
(though at a cost of the cosmological principle [85,86]), a A-dominated model was revived 
again. The cosmological constant solves a lot of problems at once (for ex., [87]). It supplies 
the 'missing matter' required to make fitot = 1, it modifies CDM by putting more power on 
large scales, as is compatible with the CMBR anisotropy limits. It also removes the inconsis- 
tency between the age of the universe and that of the globular clusters for larger values of Hq. 
However, we face the problem that the upper limit of A from observations (A <, 10"^^ cm"^) 
is 120 orders of magnitude below the value for the vacuum energy density predicted by quan- 
tum field theory [88, 89]. (It is customary to associate a positive cosmological constant A with 
a vacuum density py = A/SttG.) According to gravitational lensing statistics, a universe with 
the a large cosmological constant should have more multiple-image systems than are actually 
observed. It was shown, for example, that HDF data may be incompatible with large A [75]. 

These (and other) problems generated a lot of interest in an open FRW model with linear 
evolution of the scale factor, a(t) oc t. This universe expands with a constant speed, hence the 
term coasting cosmology. Notable among such models is a recent idea of Allen [90], in which 
such a scaling results in an SU (2) cosmological instanton dominated universe. Yet another 
possibility derives from the Weyl gravity theory of Manheim and Kazanas [91]. Here again 
the FRW scale factor appoaches a linear evolution at late times. 

There are several motivations for investigating such models. Many of the problems of the 
standard model are naturally resolved in such a cosmology. Such a cosmology does not suffer 
from the horizon problem. Horizons occur in models with a(t) ~ for a<l. Linear evolution 
of a scale factor is supported in alternative gravity theories (eg. non-minimally coupled scalar- 
tensor theories), where it turns out to be independent of the matter equation of state. The 
scale factor in such theories does not constrain the matter density parameter, thereby curing 
the flatness problem. The age estimate in a coasting universe, deduced from the measurement 
of the Hubble parameter, is given by ^ = 1 /Hq. This is 50% greater than the age in standard 
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CDM cosmology, thus making it comfortably concordant with the ages of globular clusters. 
Finally, a linear coasting cosmology independent of the equation of state of matter, is a generic 
feature in a class of models that attempt to dynamically solve the A problem [92]. Such models 
have a scalar field non-minimally coupled to the large scale curvature of the universe. With 
the evolution of time, the non-minimal coupling diverges, the scale factor quickly approaches 
linearity and the non-minimally coupled field acquires a stress energy that cancels the vacuum 
energy in the theory. This model is consistent with nucleosynthesis in the early universe [93] 
and presents a good fit to the current SNE la data [94]. 

After the discovery of the first multipy imaged quasars, gravitational arcs and arclets, grav- 
itational lensing has rapidly become one of the most promising tools for cosmology. It is now 
well known that gravitational lensing is a useful probe of the geometry of the universe. With 
the rapid growth of the number of lensed systems, proposals were made to apply the statistical 
analysis to the samples of lenses in order to extract the cosmological parameters. The use of 
gravitational lensing statistics as a cosmological tool was first considered in detail by Turner et 
al. [95] and Fukugita et al [96]. More recently, Kochanek [97] and Falco et al [98] have laid the 
groundwork for using gravitational lensing statistics for a detailed analysis of the extragalac- 
tic surveys. It has been pointed out [96, 99] that the expected frequency of multiple imaging 
lensing events for high-redshift sources is sensitive to cosmology. In view of the successfuU 
results of the above mentioned works [93,94], it is tempting to use this test to constrain the 
power index of the scale factor of a general power law cosmology, a(t) oc In the next section 
we consider constraints on this index. The expected number of lens systems depends upon the 
index a through the angular diameter distances. By varying a, the number of lenses changes 
and on comparison with the observations we obtain constraints on a. 

One of the remarkable results of field theory is the existence of stable classical states of 
a field with non-vanishing energy. Such fields include topological solitons, domain walls, 
strings and monopoles. In addition to these field configurations, which are stabilized by their 
topological properties, non-topological solitons (NTS) produced by scalar fields have appeared 
in the literature and their relevance to cosmology has been assessed (for example, [100]). 
Variations on this theme include cosmic neutrino balls [101, 102], Q-balls [103] and soliton 
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stars [104]. NTS are rendered stable by the existence of a conserved Noether charge carried 
by the fields confined to a finite region of space. The theory essentially contains an additive 
quantum number N, carried either by a spin- 1/2 field ip (for example, fermion number), or a 
spin-0 complex field cp. In addition there is a scalar field a whose coupling gives cp or ip a 
mass. The soliton contains the interior in which a ^ (Tq, surrounded by a shell of width ^ 
over which a changes from ctq to 0. ii is the mass associated with a. a = and o" = ctq are two 
minima of the effective potential of the scalar field. The A'^-carrying field i/j, or (p, is confined 
to the interior where it is effectively massless at the local minimum of the potential, a = ctq. 
At the global minumum, a = 0, the field (ip or cp) has a non-vanishing mass, m. This leads to 
a stable configuration of massless particles trapped inside a region with a = (Tq separated from 
the true vacuum a = by a wall of thickness fti fT^. Confinement occurs for all particle states 
that are not on-shell in the exterior region. The NTS is stable as long as the kinetic energy 
of the particles inside the bag is less than Eftee — the minimum on-shell energy in the exterior 
region. This was explicitly demonstrated by Lee and Wick [105, 106]. 

The role of scalar fields in effective gravity models stems from the classical work of Brans 
and Dicke. In these approaches the gravitational action is induced by a coupling of the scalar 
curvature with a function of a scalar field. Lee and Wick's results can be carried over to 
a curved spacetime with a scalar curvature R non-minimally coupled to a through arbitrary 
function U (a) in the class of theories described by the effective action 



Here V(a) is effective potential of a scalar field, Lm is the matter field(s) Lagrangian which 
includes a Higgs coupling of cr to a fermion. Establishing and (Tout as the interior and 
exterior values of a, the effective gravitational constant would be given by 
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If U goes to infinity at some point, defined as a = without loss of generality, the theory gives 
rise to a solution with a spatial variation of the effective gravitational constant, - and 
=constant. This would be a generic feature of a Lee-Wick solution in which the scalar 
field is non-minimally coupled to the scalar curvature. It would give rise to a stable "ball of 

gravity". 

A gravity ball is characterised by an effective refractive index in its interior. This would 
cause bending of light incident on it. The G-ball thus acts as a spherical lens. In the last 
section we investigate gravitational lensing of an empty gravity ball situated at a cosmological 

distance. The lens — gravity ball — has interesting features which are not shared by other known 
lenses. In the last section we examine gravitational lesning properties of this special kind of 
NTS solutions referred to as G-balls [107]. 

4.2 Gravitational lensing statistics as a test for coasting cos- 
mology 

4.2.1 Linear coasting cosmology 

We consider a general power law cosmology with the scale factor given in terms of two pa- 
rameters 



for an open FRW metric. However, the light propagates through the inhomogeneous rather 
than the averaged smooth spacetime. The light ray feels the local metric which deviates from 
the smooth FRW metric — FRW metric excludes gravitational lensing in principle. We may 
assume, however, that the relation between the affine distance of the null geodesic and the 
redshift of the source object in a clumpy universe is the same as in the FRW cosmology. This 
may be regarded as a mathematical expression of the assumption that the universe is described 
by FRW geometry on large scales (for the metric see Eq. 2.14). 




(4.2) 
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Dimensionless Hubble parameter 



The expansion rate of the universe is described by a Hubble parameter, H = d/a. Using (4.2), 

m)=j. (4.3) 

The present expansion rate of the universe is defined by a Hubble constant, equal in our model 
to Hq = a /to (here and subsequently the subscript on a parameter refers to its present value). 
With the usual definition of the redshift z 

^^=1+. (4.4) 

a(z) 

we obtain to/t - (1 +z)^/°, and thus 

H(z) = Ho(l+zy^''. (4.5) 
The dimensionless Hubble parameter is 

/i(z) = ^ = (l+z)'/". (4.6) 

-no 

Present-day scale factor 

We define the present 'radius' of the universe (4.2) as 

ao = B^. (4.7) 

-no 

Angular diameter distance 

There are many ways to define a distance in cosmology — parallax distance, angular diameter 
distance, luminosity distance, proper motion distance, etc (see 2.2.4). It is the angular diameter 
distance that is relevant to the angular separation of images. 

Consider a source at r = ri which emitted light att -ti. The observed angular diameter of 
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the source, 6, is related to the proper diameter of the source, D, by 



a(h)ri 

The angular diameter distance Da is defined to be 

DA(z) = y = aai)ri = -^, (4.9) 
1+z 

where ri{z) is the coordinate distance to a redshift z- Angular diameter distance of a source at 
Zi measured by an (fictitious) observer at zi is given by 

D^{zuz,)=^. (4.10) 

Light travels on null geodesies, ds^ = 0. If a comoving observer is at coordinates (ro,Oo,(f)o) 
at time to, the geodesies intersecting tq = are lines of constant 9 and 0. Geodesic equation 
becomes 

dr^ 

Q^cdt^-a\t)--—^ (4.11) 

and the light signal emitted from coordinates {ri,9, (p) at time t will reach the observer at time 
to determined by 

f'^cdt' _ n dr 



Vl + r2 

We can convert the integrals over t into integrals over z by differentiating (4.4), 

dz ao . aod 

-J- = — 2^ = = -(i+z)H(z) , 

at a a a 



(4.12) 



to obtain 



Thus, equation (4.12) becomes: 

cdz dr 



H{z){l+z)a io \/T+r2 



(4.14) 
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Since ao/« = 1 + z> we use it to obtain 



c_ n apdz ^ c_ n dz(l+z) r 
aoJo H(z)(l+z)a aoJo H(z)(l+z) floio 



dz 



or, in terms of a dimensionless Hubble parameter h{z), Eq. 4.12 becomes 



dz n dr 



HoaoJo h(z) Jo Vl + r^ 



Comoving Coordinate Distance r For an open universe, ^ = -1, we obtain 



f - sinh ^ ri . 
io h(z) 



Thus, 



ri = sinh 



dz 



Da formulae If we take the value of ao from Eq. 4.7, we obtain 



ri = sinh 



1 n dz 



BJo Hz), 



Then, 



with the integral 



Da(z, a) = 



Be 



(l+z)Ho 



sinh 



1 r dz 



BJo h(z). 



Finally, the distance formula is 



Da(z, a) = 



Be 



(l+z)Ho 



sinh 



7i(l+z) a -1 « ^ 

a-l I J 



1 a 
B 



The distance between two objects at different redshifts is: 



DA(zi,Z2,a) = 



Be 



{l+Z2)Ho 



sinh 



1 a 

b'^ 
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(1+Z2)^-(1+Zl)" 



Limiting (a — >■ 1) formulae We know that h(z) = (1 with I, thus, h(z) - (l+z), 

then the integral 

/ -r— = ln(l+z), (4.20) 

and the Eq. 4.18 reduces to 



Be 1 
^A(z) = —-—smh. 
Ho l+z 



^ln(l+z) 



Be [(l+z)2/«-l] 
2/fo(l+z)^ 



The distance between two objects at different redshifts is, accordingly. 



Da(Zi,Z2) = 



Be [(l+Z2)^/"-(l+Zi)^/«] 
2^0 (l+z0i/«(l+z2)¥ 



(4.21) 



(4.22) 



Look-back time 



The look-back time. At = \t-to\, is the difference between the age of the universe when a 
particluar light ray was emitted and the age of the universe now: 



At(a;ao) = 



da 

'a 



(4.23) 



We can obtain the expression for the look-back time from Eq. 4.13 and Eq. 4.6. Since we 
define to as zero and time is increasing as we look back, we drop the negative sign in that 
expression. Thus, 

dt 1 

(4.24) 



dz //o(l+z)- 



and finally the expression we seek is: 



edt 



dZL Hoil+ZL)' 



(4.25) 



4.2.2 Basic equations of gravitational lensing statistics 



Given a lens model, the number counts of galaxies-lenses and their properties, we compute the 
probability pi that quasar / is lensed and the probability Pi(A9i) that quasar / is lensed and has 
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image separation A^,. We shall consider only early types of galaxies, elliptical and lenticulars, 
E/SO, neglecting the contribution of spirals as lenses. This is because the velocity dispersion 
V of spirals is small compared to E/SO galaxies, outweighing the larger number density. The 
numbers of E/SO and spirals are roughly 30% and 70%, respectively, whereas in our adopted 
sample the velocity dispersion of ellipticals is about 200 km s"^ and that of spirals is about 130 
km s"^ [108]. The calculations for the lensing statistics depend strongly on this quantity, since 
the lensing optical depth goes as the fourth power of the velocity dispersion (see Eq. 4.45). 

It is increasingly clear [97] that E/SO galaxies are effectively singular. The almost uniform 
absence of central images in the observed lenses and models of individual lenses imply core 
radius smaller than 100 pc and essentially suggest that the galaxy lenses are nearly singular 
[24, 109]. We will model the lensing galaxies as singular (= zero core radius) isothermal 
spheres (SIS). The detailed lensing properties of singular isothermal spheres are described in 
Section 2.3.5. Here we simply state the most important relations for the present calculations. 
A SIS lens with a velocity dispersion v at a redshift zl will produce two images of a quasar at 
a redshift zs separated by an angle 

A9 = 2aa^, (4.26) 

if angular position of the source is less than the critical angle P^i = c^d^Ls/^s- The deflection 
angle, a^, is given for all impact parameters as: 

Q;d = 47rQ'. (4.27) 

The angular separation of the two images is independent of the impact parameter as long as it is 
small enough to produce two images. The critical impact parameter is defined by aa = Di^P^r 
or, with the help of the previous definitions, 

«cr = 47r - — - — . (4.28) 

The quasars are treated as point sources of radiation. We assume that no evolution of galaxies 
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or quasars with cosmic time (quasar evolution is irrelevant, since the quasar redshift distribu- 
tion will be drawn from observed samples). In other words, the comoving number density of 
lenses is conserved 

ni. = no(l+ZL? , (4.29) 

where no is an average comoving density measured at the present epoch. Merging between 
galaxies and infall of surrounding mass onto galaxies are two posible processes that can change 
the comoving density of galaxies and/or their mass. Under the generic relation between the 
velocity dispersion and mass of early-type galaxies, Rix et al [1 10] and Mao & Kochanek [111] 
found that merging and/or evolution does not significantly change the statistics of lensing. 



Cross-section and optical depth of lensing. 

The cross-section a for "strong" lensing events is given by cr = nal^. Using (4.27) and the 
definitions from the previous subsection, we write the cross- section as 

The cross-section is largest approximately halfway in the distance between the source and the 
observer; and a vanishes at the two endpoints. To find the average effect on an image passing 
within Ucr of a mass ("scoring a hit"), we average over the cross-section to find the mean image 
separation as 

A0 = Sn(-Y— . (4.31) 
\cJ Ds 

The differential probability (or the optical depth) dr that a line of sight intersects a galaxy-lens 
at Zl in traversing the path of dzh from a population with number density is given by a ratio 
of differential light travel distance cdt to its mean free path between successive encounters 
with galaxies, l/ni^(z)a, 

cdt 

dr = niiz)a-—dzL ■ (4.32) 
dzi^ 
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For our case ni^(z) is given by (4.29) and the quantity cdt/dzh is obtained in the previous 
section (Eq. 4.25). Substituting for a and ni^(z), we get 



<ir = — -,z.. (4.33) 

Assuming that the brightness distribution of galaxies at any given redshift is described by a 
Schechter function, the comoving density of galaxies at redshift z and with luminosity between 
L and L+dL is 

*(L,z)<iL = „.(z)(^)"exp(-^)£L. (4.34) 



The parameter n*(z) is the average comoving density, L*(z) is the characteristic luminosity at 
which the luminosity function exhibits a rapid change in the slope in the logn, logL)-plane 
and a gives the slope of the luminosity function in the (logn, logL)-plane, when L<^L^. The 
comoving number density of galaxies, characteristic luminosity and mass of a galaxy at any 
redshift remain constant, therefore, n*(z) = n*(0) = constant and L*(z) = L*(0) = constant. "0" 
refers to present-day values. This is the most commonly used luminosity function for early 
type galaxies. Thus, Eq. 4.34 becomes 

$(L,z = 0)JL = n,(^^^ exp(^--^^^. (4.35) 

where n*, a and are the normalization factor, index of the faint-end-slope and the charac- 
teristic luminosity at the present epoch, respectively. These values are fixed in order to fit the 
current luminosities and densities of galaxies. From Eq. 4.33 and Eq. 4.35 we have 

We assume the velocity dispersion v is related to the luminosiy L by the empirical Faber- 
Jackson relation for E/SO galaxies: 

(v*) ' (4-37) 
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L\'"^ f L\dL 



Therefore, Eq. (4.36) becomes 



By integrating it over the luminosity function at redshift z, we obtain the differential probability 

167r^, , A^f. ^ \ f Di^Di^sV cdt 



dr=^(l + zO'ny,r\^a+- + lj\^-^^j ^dzi. (4.39) 

where F is the normal gamma function. If we define the 'dimensionless' image splitting as 
(f) = A9/8tt (v*/c)^, we can find the differential optical depth of lensing in traversing dzh with 
the angular separation between (f) and as ■^^d(j)dzL- Using luminosity- velocity relation 
(Eq. 4.37), we obtain the relation 



L f D 



Taking the differential. 



7/2 

(4.40) 



dL ^(DsX"dJ>^ ^^^^^ 



U 2 V£>LS 
Returning to the Eq. 4.33), we write dr/dzL as 

dr 167r3 .fDJh^Vcdt 

-i-=nL(z)^-v — - — — . (4.42) 

dzL \ Ds J dzL 

The differential optical depth of lensing in traversing dz^ with the angular separation between 

and (f)+d(j) is 



(4.43) 

4,°°fLY ( L\dLfLY^^ .16tt^ fDi^Di^sV cdt 







n^vl / - exp -- — - (1+Zl) 



c4 V Ds ; Jzl 



100 



Using the relations (Eq. 4.41), we obtain 



d(j)dzL = 



F*(1+Zl) 



3 Ho 



7/2 



cdt ( HqDi^Di^s 



LS 



exp 



D 



Is 



(4.44) 



—rdZL- 



Here we have introduced the useful dimensionless quantity F* which measures the effective- 
ness of matter in producing multiple images [95] 



167r3 .[^ , 4 



(4.45) 



Lens and source parameters. 



The variables relating galaxy numbers counts and the isothermal lens model are the number 
density of E/SO galaxies n*, the Schechter function slope a, the Faber- Jackson exponent 7 and 
the velocity dispersion of the dark matter for an galaxy v*. The Schechter function slope a 
controls the relative number of low- and high-mass galaxies. We use the following list of lens 
parameters for our calculations [108]: 



Survey 


a 7 v*{Km/s) (f)*(Mpc-^) F* 


LPEM 


+0.2 4.0 205 3.2x10-3 0.010 



We consider a sample of 867 (z > 1) high luminosity optical quasars which include 5 lensed 
quasars (1208+1011, H 1413+117, LBQS 1009 + 0252, PG 1115 + 080, 0142+100). This 
sample is taken from optical lens surveys, such as the HST Snapshot survey, the Crampton 
survey, the Yee survey, Surdej survey, the NOT Survey and the FKS survey [112]. The lens 
surveys and quasar catalogues usually use V magnitudes, so we transform my to a B-band 
magnitude using 5-V = 0.2 [113]. 

We make two corrections to the optical depth to get the lensing probability: magnification 
bias and selection function. 
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Magnification bias. Lensing increases the apparent brightness of a quasar causing over- 
representation of multiply imaged quasars in a flux-limited sample. This effect is called mag- 
nification bias. It is an enhancement of the probability that a quasar is lensed, that is why we 
have to include this correction in the probability function. The bias for a quasar at a redshift z 
with apparent magnitude m is given as [96, 97, 1 14] 



5(m,Au^,A2,z)= — J -] / — ^—^p(A)dA. (4.46) 

^"^ J Ja^ dm 



Here A is the total magnification, Ay^a the magnification at which the images have the mini- 
mum detectable flux ratio (for a SIS model, Aum =Aq = 2), m^ = m-i-2.51ogA is the lensing 
enhanced magnitude and ;?(> A) is the probability distribution for the two images having a to- 
tal magnification larger than A, p{A) = 8 /A^ for a SIS model. We can allow the upper limit on 
magnification A2 to be infinite, but in practice we set it to be A2 = 10"^. For the quasar apparent 
magnitude number counts we use Kochanek's broken power law [97] 

dN{m,z) ^ , ^ Q-a{m-mo{z)) ^. j Q-bim-m>(z))\ "1 (4.47) 

dm ^ ' ' 

where the bright-end slope index a and faint-end slope index b are constants and the break 
magnitude mo evolves with redshift as 

mo + (z+l) if z<l 
mo(z) ^^mo ifl<z<3 
mo-0.7(z-3) if z>3. 

Fitting this model to the quasar luminosity function data in [115] for z > 1, Kochanek finds 
that the 'best model' has a = 1.07, b = 0.27 and mo = 18.92 at B magnitude. Thus, magnitude 
corrected probability pi for the quasar / with apparent magnitude m, and redshift z, to get lensed 

is Pi = r(z,)5(m,-,Aiim,A2,z,). 



Selection function. Selection effects are caused by limitations on dynamic range, limita- 
tions on resolution and presence of confusing sources such as stars. In the SIS model the 
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selection function is modeled by a maximum magnitude difference Am(6) that can be detected 
for two images separated by A^. This is equivalent to a limit on the flux ratio between two 
images / = lo^-^^'^^^). The total magnification becomes Af = Ao(f+ 1)/(/- 1). So, the sur- 
vey can only detect lenses with magnifications greter than Af. This sets the lower limit on 
magnifications and in the bias function An^ gets replaced by Aj. To get selection function 
corrected probabilities, we divide our sample into two parts: the ground based surveys and the 
HST Snapshot survey and use the selection function for each survey as suggested by Kochanek 
[114]. 

The corrected lensing probability and image separation distribution function for a single 

source at redshift zs are given as [97] 

Pi(m,z)^Pi / \ , (4.48) 

J 5(m,z,Mo,M2) 

and 

,:,=,,(A«,Z!:M^, (4.49) 

p. 5(m,z,Mo,M2) 

where 

p,(Ae) = ^ r '^Ij.^ dzL. (4.50) 

Equation (4.49) defines the configuration probability. It is the probability that the lensed quasar 
i is lensed with the observed image separation. 

4.2.3 Testing the model against observations 

The above basic equations were used to perform the following test: 

The sum of the lensing probabilities p[ for the optical QSOs gives the expected number of 
lensed quasars, = 5^ Pi- Th^ summation is over the given quasar sample. We look for those 
values of the parameter for which the adopted optical sample has exactly five lensed quasars 
(that is, those values of the parameters for which hl = 5). 

We started with a two parameter fit. We allowed a to vary in the range (0.0 <a< 2.0) and 
B to vary in the range (0.5 < 5 < 10.0). We observe that for 5 > I, Da becomes independent 
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of it (as is also obvious from equation (4.18)). In the previous works constraining power law 
cosmology [93, 94], the value of 5 = 1 was found to be compatible with the observations. 
Incidentally, we can estimate the present scale factor of the universe as oq ~ c/Hq. Therefore, 
we used 5 = 1 in further analysis. 

Fig. 4.1 shows the predicted number of lensed quasars for the above specified range of a. 
We obtained = 5 for o; = 1.06. We further generated 10^ data sets using bootstrap method 
(see description in [116]) and found the best fit for a for each data set in order to obtain error 
bars on a. We finally obtained a = 1 .09 ± 0.3. 
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Figure 4. 1 : Predicted number of lensed quasars in the adopted optical quasar sample, with 
< 4", vs. power index a. 



The maximum likelihood analysis was also performed [116] to determine the value of a, 
for which the observed sample becomes the most probable observation. It was found that 
0.85 < a < 1.56 at Icr (68% confidence level), and 0.65 <a< 2.33 at la (95.4% confidence 
level). It is also interesting to notice that for a general power law cosmology the simplest 
constraint on a comes from the relation a = //q^o (Eq. 4.3). With updated value of //q = 70 ± 7 
km/sec/Mpc [80] and = 14 ± 2 Gyr [1 17], this constraint gives a = 0.98 ± 0.25. 
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4.2.4 Concluding remarks 

In this section we discussed the general power law cosmology and the ways to constraint its 
power index. We have derived angular diametere distances for this cosmology. We described 
the general formulation of the gravitational lensing statistics and used it to constraint the pa- 
rameters for that cosmology. These results together with other tests were used in [4] in the 
combined form and the conclusion made was that open linear coasting cosmology, a(t) = t, is 
consistent with the present observations. 



4.3 Gravity balls as gravitational lenses 
4.3.1 Gravity ball as a NTS solution 

Here we will discuss the general formalism of the problem for a spherically symmetric system 
consisting of fermion field ^, scalar field a and gravitational field g^,y. We follow the theory 
presented by Sethi and Lohiya [107]. Gravity balls are the NTS solutions of the field equations 
arising from (1) 

u(a)r - ^gi^^R] = -- [rj''^ + + r^^; + + w{aV^ - 2g'^'^c/(a)5] (4.5 1) 

dY dU 

g^''a,,,, + — -R— = (4.52) 

Here T^", T^^, and T^'^ are energy momentum tensors constructed from action for the 
scalar field, the fermion field, together with its Higgs coupling to a, and the rest of the matter 
fields, respectively. We consider an NTS with the scalar field held to a value do in the interior 
and making a fast transition to u = outside a thin shell. Thus, we have essentially three 
regions: (i) the interior of a soliton (r < Rq); (ii) a shell of thickness ~ fT^ and surface energy 
density s ^ /x(7o/6; and (iii) the exterior (r > Rq). The total energy of a NTS has contributions 
from: (1) the surface tension energy ~ sRq ; (2) the energy of the fermions Ef ^ N^l^ /Rq 
and; (3) the volume energy Ey ~ y((Jin)i?Q. For the degenerate case Vioi^) = 0, a NTS has total 
mass constrained by the stability against gravitational collapse to a value determined by the 
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surface tension s. The soliton mass, obtained by minimizing the total energy, is M = llnsRl. 
For s ~ (MeV)^ and N ^ 10^^, the size of the NTS is of the order of tens of kiloparces, while it 
is still away from the Schwarzschild bound. For configurations with R much greater than the 
Schwarzschild radius, the effects of gravity can be treated as a small perturbation. Thus, the 

form of the metric for the NTS satisfies a weak-field approximation. 



A. Interior: r < Rq + O {^fj, ^) . The metric inside is described as 



(4.53) 



The interior metric has one specific solution: v = -Cr^/6 and 



M = Mo + ■ 



c c 

— + — 
2 3 



(4.54) 



with C and C — constants depending on the fermionic energy inside the soUton. In the weak 
field approximation used here we consider the interior with 



ds^ e^"" - dr^ - [dO^ + sin^ d^^] 



(4.55) 



B. Exterior: r > Rq + O (^f^ In the exterior region we have essentially a Minkowskian 
metric 

ds^ = dt^ - dr^ - [dO^ + sin^ Od^p^] (4.56) 

To be consistent with the observations we will see that the uq has to be a small negative con- 
stant. The propagation of light inside the G-ball is equivalently described by using the Fermat 
principle with the effective refraction index Reff inside the ball given by ^eff = 1 - mq- This gives 
straight trajectories of light rays inside the ball. But since outside the ball Hgff = 1 , we see that 
G-ball behaves as a spherical lens (deflection of light occurs only at the boundaries). 
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4.3.2 Lens model for a G-ball and lensing properties 



We restrict our discussion to a cosmological model which is a variant of a Milne universe [107] 
with a scale factor a(t) = In this cosmology angular diameter distances are given by 

czi (2+zi) 



CZs (2 + Zs) 



2 ' 



2Ho(l+z,) 
with 

Zls = T— , 
1+Zl 

and/fo = lOO/i km s"^ Mpc"^ Subscripts '1', 's' and 'Is' stand for lens, source and lens-source, 
respectively. 

We consider the gravity ball to be of the size of a typical cluster of galaxies at the cosmo- 
logical distance and embedded in empty space region with no matter concentrations close to it. 
The deflector is transparent. We also assume a thin-lens approximation, since all the deflection 
occurs within Az < ±i?bau; the extent of the deflector is thus taken to be small compared with 
its distance from both the observer and the source. Typical distances from us to the cluster at 
z = 0.3 and to the source at z = 1 are ~ 1 Gpc and ~ 2 Gpc, respectively, while i?baii ~ 0.5 to 
1 Mpc. We consider deflection angles to be small. We also assume that the source, lens, and 
observer are stationary with respect to comoving coordinates. 

The magnification of images in the axisymmetric case is given by 

The tangential and radial critical curves folow from the singularities of the tangential and radial 
magnification 

...far „.f^V. (4.59) 
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Figure 4.2: Geometry of lensing by a gravity ball. i?baii — radius of the ball, i — angle of inci- 
dence, r — angle of refraction, — deflection angle, n — ratio of the refractive index inside the 
ball to the refractive index outside; n> I. 

Lens equation In the Fig. 4.2 we present the geometry of G-ball as a lens. From this figure 
and Snell's law the deflection angle is 



Since we consider all angles to be small, the following approximations are valid. Defining the 
quantity Oq = ^bau/^L, which is the radius of the region inside which the refraction occurs, we 
obtain the the expression for the deflection angle 




(4.60) 




(4.61) 



Thus, we obtain the lens equation for the gravity ball 




for < 6* < ; 



for ^ > • 



(4.62) 



where 9 =\d\ = ^d\-\-Q\ is the radial position of the image in the lens plane. 
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Multiple image diagram and conditions for multiple imaging. To illustrate the lensing 
properties described by equation (4.62) we show in Figure 4.3 the lensing curve for the gravity 
ball with parameters n = 1.0005, 9c = 5' and for Zsource = 1> Ziens = 0.3. The intersections of 
the lines P = const with the curve given by equation (4.62) give the solutions to the lensing 
equation. The source at P2 lies on a caustic point and Oj- is the radius of the radial critical curve. 
The source at (3i has three images, a point (3 = produces a ring — 'Einstein ring' — with the 
critical radius 9 = 9t and an image at ^ = 0. The function /3(9) is continuous except at 9c due to 
the edge of the ball. 




-2x10 ^ 2x10 ^ 

9 [rodians] 



Figure 4.3: Solution to the lensing equation by a gravity ball. The solid curve represents the 
lensing equation curve with n = 1.0005, 9c = 5', Zsource = 1. Ziens = 0.3, together with lines (3 = Pi 
(dashed lines) for various source positions The intersections of the lines /3 = const with 
the lensing equation curve give the number and positions of the lensed images. The source 
at P2 lies on a caustic point and 9^ is the radius of the radial critical curve; source at Pi has 
three images and point /? = produces a ring 9c = 9{, in addition, it has an image at ^ = 0. The 
function P{9) is continuos except at 9c due to the edge of the ball. 



Putting /3 = in (4.62) gives 



^-^fsin-f-sin-4)=0. (4.63) 
Ds V n9c, 
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Denoting a - 2Z)ls/^s> we obtain 

e = a( sin-i ^ - sin-i 4r I • (4-64) 
V nOcJ 

When the aUgnment of the source, lens and the observer is not perfect, we see features called 
arcs, which are the result of very strong distortion of a background sources. Arcs roughly 

trace the Einstein ring, so ^?arc ~ ^e- Since the radii of most known arcs do not exceed 30" 
small angle approximation is valid and we obtain 



One solution to this equation is trivial: 

^ = , (4.66) 

which corresponds to the image located at the centre of the lens. To find other solutions we 
write the expression in the curly brackets in (4.65) as 



and in order to find the conditions for the Einstein ring we rewrite it in the form 



Assuming for simplicity a=l (the ball (lens) half-way between the observer and the source) 
and using n from other estimates [107] to be from 1.001 to 1.0001 we obtain 
for n= 1.0001 

^E = 0.98^c, 

for n= 1.001 

= 0.73^c. 
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Here 9e is non-trivial solution of the equation (4.65). This leads to the appearence of the image 
in the form of a ring with the radial size of ~ 0.9 to 0.7 of the total size of the ball, depending 
on n. To find the condition for the appearence of multiple images we analyze (4.67). To obtain 
a physical solution we take 



1 



2 2 

an 



1 



92 



1 



a 



>0 



(4.69) 



The conditions for multiple imaging are 



^ball > 2Z)efr (1 — 
n 



(4.70) 



or expressed through n 



1 <n< 



2D, 



eff 



2Deff--/?ball 

where Z)eff = Z)ls^l/Ds is the effective distance. 



(4.71) 



Magnification and critical curves. From lens equation for a G-ball (4.62) and equation 
(4.58) we obtain expression for the total magnification of the images 



sin ( — -sm 



1—a 



1 



1 



, (4.72) 



where the first term represents /Xt and second term — //r. In Figure 4.4 we have plotted tangen- 
tial magnification fit and radial magnification fi^ against 9 for a G-ball with refraction index 
n = 1.0006. Singularities in these give the angular positions of tangential critical curves and 
radial critical curves, respectively. In the same Figure we have also plotted total magnification 
/X vs. 9. 



Limits on n from observations. We can rewrite equation (4.68) in terms of variables i?bau 
and Re to get an expression for i?bau in terms of n and a given Re. 



nDes 



(4.73) 
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n=1.0006, R,,„=1.5x10-^ (rad), Z3„„^^^=1. z,^^=0.3 
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Figure 4.4: The magnifications: tangential /ij denoted by green, radial /i^ denoted by red, and 
total is shown by blue line. Curves are plotted as a function of the image positions 6 for the 
parameters of the lens: n = 1.0006, 6c = 5', zi = 0.3, Zs = 1- The singularities of //r and fit give 
positions of the tangential and radial critical curves, respectively. 



Since i?baii > Re, we get a lower bound on the value of 7?baii if we know the radius of the Einstein 
ring. It is clear from the above equation that if the size of the G-ball is big enough then a large 
enough n can give us any desired radius for the Einstein ring. If we assume that all G-balls 
to be of the same size, then from observations we can infer a lower bound on the radius and 
the refractive index inside a ball by the following argument. The radius of any gravity ball has 
to be larger then the radius of the largest observable Einstein ring. Given this size of the ball 
the refractive index should be large enough to give a real Einstein ring for every other case. 
This situation is illustrated in the Figure 4.5, where observational data for several clusters with 
giant arcs (details are presented in the Table B.l, App. B) are used with the assumption that 
the raduis of the arc 6'arc ~ ^e- Together with the curve (4.73) for the cluster A370, we plotted 
the value of the radius i?arc of the A5 arc in that cluster (horizontal line), which is the largest 
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amongst the clusters considered. It is clear that the refractive index should be greater than the 
value where this line intersects the ^baii-" curve for this cluster. Assuming all radii to be the 
same, we can see that in order to have an arc, a gravity ball must have n more than n ^ 1 .00056. 
Thus, we obtained a lower limit on n. Of course, our present calculations are for the case of an 
empty gravity ball. Matter concentrations in the centre of the ball will increase the deflection 
angle. 




1.0002 1.0004 1.0006 1.0008 1.001 



n 

Figure 4.5: Rb^n vs. n for different clusters. Clusters are marked by different colours. Cor- 
responding curves have the same colour as the name of the cluster. Details are presented in 
Table B.l (App. B). Horizontal line is the value of radius of the arc A5 in the cluster A370. 



4.3.3 Simulations results 



To demonstrate how a G-ball at intermediate redshift gravitationally distorts background 
sources we performed computer simulations and present the results in Figures 4.6 and 4.7. 
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Figure 4.6: Illustration of six different images connnfigurations for an empty gravity ball. 
Zi = 0.3, Zs = 1, = 5'. Depending on the source position and n the lens produces images of 
(a) ring, (b) single arc, (c) arc and opposite image, (d) radial arc, (e) three images (with two 
opposite arcs), (f) two images on one side (with a straight arc). Small dot in the centers of the 
panels marks the center of the G-ball. Small star marks the position of the source. 
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Figure 4.7: Distortion field generated by a simulated gravity ball. The left panel shows the 
grid of randomly distributed background sources as it would be seen in the absence of the lens. 
The middle panel shows the same population once they are distorted by a foreground gravity 
ball with the parameters: Zs = 1, Zi = 0.3, 9c = 5' rad, n= 1.001. The right panel shows the same 
population distorted by a foreground (invisible) circular cluster (zs = 1.3, zi = 0.4, a = 1000 
kms~^y, the units for both axis are in arcsecs. 

We have simulated an empty gravity ball the size of a typical galaxy cluster at the reshift of 
Zi = 0.3, assuming for simplicity that all source galaxies are at the same redshift Zs = 1- The 
luminous area of each background galaxy is taken to be a circular disk of radius R with con- 
stant brightness. In all our numerical computations we used dimensionless Hubble parameter 
h - I. For the simulations of G-ball lensing we use the algorithm described in [73]. The 
image configurations by a simulated gravity ball are illustrated in Fig. 4.6 by changing the 
relative positions between the extended circular source and the lens. When the alignment of 
the source on the optical axis is perfect, the ring image — Einstein ring — appears (Fig. 4.6(a)). 
With a small displacement we see two opposite arcs approximately located at the Einstein 
radius (Fig. 4.6(e)). When the displacement becomes larger different features may appear: 
single tangential arc (Fig. 4.6(b)), radial arc (Fig. 4.6(d)), two opposite images (Fig. 4.6(c)) 
and straight arc (Fig. 4.6(f)). Figure 4.7 displays the distortion the gravity ball produces on the 
field of randomly distributed background sources. Figure 4.7(a) shows an unperturbed galaxy 
backgound projected randomly in a field of ~ 10' x 10' at a redshift of z = 1. Fig. 4.7(b) illus- 
trates the same field with an empty G-ball at the centre. To compare with the lensing effect 
a cluster modelled as a singular isothermal sphere produces on the background galaxies, we 
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presented in the Fig. 4.7(c) the simulation taken from [21]. We can see from the figures that an 
empty gravity ball can indeed act as a strong lens and produce the images of arcs and arclets. 

4.3.4 Concluding remarks 

In this section we have discussed a special kind of nontopological soliton, known as gravity 
ball, and examined the possibility for it to be a gravitational lens. We investigated its lensing 
properties, calculated the deflection angle of a point source, derived the lens equation and plot- 
ted the lens curve for the gravity ball with specific parameters. We have shown that depending 
on the parameters of the model there can be one, two or three images (two inside the Ein- 
stein radius and one outside). However, the G-ball gravitational lens has interesting features 
which are not shared by other known gravitational lenses. In the case of a source, lens and 
observer located on the optical axis together with the Einstein ring there is always an image in 
the centre. Besides, the lensing geometry is radically different — there is no effect outside the 
lens and thus, the Einstein radius is always less than radius of the ball. The refractive index 
of the gravity ball contributes to the decrease in the radius of the Einstein ring and increase 
of its thickness. For the extended sources we showed how a large gravity ball can induce sur- 
prisingly large distortions of images of distant galaxies. It is assumed that in order to produce 
large arcs the cluster must have its surface-mass density approximately supercritical, S > T^crtt- 
We found that the empty balls alone are able to produce arc-like features, including straight 
arcs, arclets and radial arcs. We have modelled the gravitational lensing effect of our gravity 
ball on a background source field and found that the ball produces distortion ('shear') of that 
field, consistent to some extent with the observations. We also obtained constraints on the size 
of the large gravity ball, which can be inferred from the existing observations of clusters with 
arcs. Empty G-balls can show in observations as arcs without evident mass distribution — a 
lens and/or through the distortion of the background random galaxy field. 

Observationally empty G-balls have not yet been found, but we conclude from our analysis 
that the existence of large G-balls cannot be ruled out by gravitational lensing effects. 
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Chapter 5 

Propagation of Light in Strong 
Gravitational Fields 

5.1 Introduction and main motivations 

As described in the Introduction to the thesis (2.1), gravitational lensing in strong gravitational 
fields leads to many astrophysically interesting effects. In this section we discuss yet another 
one. 

The idea that centres of the globular clusters host massive (up to 10"^ Mq) black holes (BH) 
goes back to 1975, when scientists were trying to detect them through X-ray emission. Since 
then other methods were proposed. In particular, the observations of the millisecond (ms) pul- 
sars in globular clusters can help, since the pulsar motion would be governed by the presence 
of massive black hole in the centre of the globular cluster ([125], [121] and references within). 
A large number of ms pulsars have been discovered in globular clusters — more than 35 ms 
and 4 long-period ones in 13 Galactic globular clusters have been reported. This is not sur- 
prising, since the dense stellar field in the globular clusters provides a favourable environment 
for frequent stellar encounters to form low-mass-X-ray-binaries (LMXB). These are believed 
to be the progenitors of the ms pulsars. By some estimates (for ex. [118]), the total number 
of pulsars in Galactic globular clusters, beaming at us, most likely lies between 300 and 2000. 
This provides an excellent opportunity to use them for revealing the existence of central black 
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holes. A possible mechanism for establishing the presence of massive black hole in the centre, 
which has not so far attracted any attention, is through the gravitational lensing of the pulsar's 
radiation by the central mass. It is only logical to suppose that the black hole in the centre of 
the globular cluster would be a rotating black hole. There exist a number of papers investi- 
gating the lensing properties of Kerr black holes. Bray [119] was the first to investigate the 
properties of the isolated Kerr black hole as a gravitational lens. He demonstrated that effect of 
spin of the BH comes in the second order terms of the expansion of null geodesies in the Ken- 
metric and considered the conditions for the multiple imaging. Simulteneously, Dymnikova 
[122] studied the time dilation, in the sense of Shapiro effect, and showed how the asymmetry 
in the time delay depends on the mutual orientation of a photon propagation direction and of 
the rotation axis. The photon is accelerated, if its motion is along the direction of the rotation, 
and it is delayed otherwise [124]. This effect can be traced back to the dragging of inertial 
frames. The asymmetry's is very enhanced in the regions corresponding to lensing and in- 
creases sharply for the photons traveling very close to the central body. If, for example, the 
motion is retrograde, the effect of the rotation is to impose an additional time delay acting in 
the same direction as the effect arising from the mass. For systems with the central rotating 
black hole the time lag can be of the order of minutes [123]. 

One more interesting effect was noted by Carini et al [120]. They showed that angular 
momentum of a Kerr black hole gives rise to a rotation of the plane of polarization of a lin- 
early polarized light. They noted that the rotation of the polarization plane induced by the 
angular momentum of the central gravitational source is analogous to Faraday effect of the 
electromagnetic wave propagating through magnetized plasma and thus may be referred to as 
the "gravitational Faraday rotation". 

Still more effects may also possibly exist — in particular, it may be interesting to study 
the differential frequency shift between the beams propagating in the equatorial plane (where 
this effect would be maximal) on the opposite sides of the black hole. Secondly, for a source 
of pulsed radiation (e.g. pulsars), it would be interesting to determine if a mechanism exists 
which can produce the difference in the time period between pulses in two images. In that case 
it could be possible that some (?) of the pulsars in certain globular clusters are gravitationally 
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lensed images of the same source. This is the main motivation with which the present work 
has started. 



5.2 Photons in a gravitational field 

It was Galileo's famous experiment that asserted that big balls fall at the same rate as small 
balls. By the more general statement, known as EP (Einstein equivalence principle) all par- 
ticles follow the same geodesic path, independent of the mass or internal structure. Photons 
move along null geodesies in curved spacetime, and all photons are deflected in a gravitational 
potential by the same angle independent of their energy or polarization. This is mostly used 
in studies of the gravitational lensing around black holes. However, some special situations 
are reported, where gravitational birefringence and dispersive effects are described. For ex- 
ample, in quantum electrodynamics, virtual electron loops induce curvature couplings in the 
effective action of low-frequency photons [126]. In such a case, photon propagation can be 
influenced by tidal effects coming from local spacetime curvature. Drummonds and Hathrell 
[127] found that the propagation of photons in spatially anisotropic background (Schwarzchild 
and gravitational-wave) was polarization-dependent. Similarly, gravitational birefringence 
was also found to occur in Reissner-Norstrom spacetimes [128]. The energy-dependent de- 
flection of light was shown in the context of string theory [130] and higher derivative gravity, 
which arises from QED radiative corrections [126]. There is some controversy in literature. 
For example, Masshoon firmly states that birefringence is possible only if the central body 
is a rotating body and it is a consequence of the coupling of the intrinsic spin of the photon 
with the rotation of the central body. That is, for a Schwarzschild black hole the amplitudes 
for the scattering of right circular polarized (RCP) and left circular polarized (LCP) waves are 
equal and, hence, the spherical symmetry of this field preserves the polarization of the incident 
radiation in the scattered waves. However, for a Kerr black hole RCP and LCP radiations are 
scattered differently and the deflection of the radiation by a rotating mass becomes polarization 
dependent [131]. 
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5.2.1 Spinning particles 



We would like to pay special attention here to the fact that photon is an elementary particle 
with spin. Spin angular momentum of matter, occurring in nature in units of h/2, is the 
physical notion which seems necessary for a successful extension of GR to microphysics. 
All elementary particles can be classified by means of irreducible unitary representations of 
the Poincare group and can be labeled by mass m and spin s. Mass is connected with the 
translational part of the Poincare group and spin with the rotational part. Mass and spin are 
elementary notions, each with an analogous standing not reducible to that of the other. 

The problem of the motion of a particle with internal angular momentum (spin) in external 
gravitational field consists of two parts: spin precession and spin influence on the trajectory 
of motion. Gyroscope precession in a centrally symmetric gravitational field was considered 
more than 70 years ago [132] and is now a well-described effect, since it is extremely impor- 
tant in the satellite industry. The situation is different with the second part of the problem, 
which refers to the spin influence on the trajectory. According to EP a spinning mass should 
fall at the same rate as non-spinning one. However, in a seminal paper, Papapetrou [133] 
showed that trajectories of spinning particles (in the limit of a vanishing mass) deviate from 
geodesies. His derivation was based upon taking a macroscopic mass, expanding it in moments 
about the center, and deriving the equations of motion for each moment. The basic result was 
that the deviation is proportional to the spin of the particle (vanishing for a zero spin) and the 
curvature of the space. It is presumed that the results are valid for fundamental particles, but 
there is no classical theory that yields the result [134]. The difficulty is that there has not been 
a satisfactory theory which derives these equations from a variational principle. There have 
been numerous attempts, though; a good review on equations for particles with spin is recently 
presented by Frydryszak [135]. Most recent claims, involving a generalized Clifford algebras, 
are presented by Pezzaglia [136]. The problem of the influence of spin on the trajectory of a 
particle in an external gravitational field is not merely of a theoretical interest. This problem 
is closely related to the gravitational waves field, where one has to calculate the gravitational 
waves radiation from the inspiralling black (Kerr) holes. In gravitational lensing this problem 
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can be of the importance. It could lead to the deviation of the photon trajectories from the 
geodesies and to the gravitational birefringence in an arbitrary gravitational field. In the next 
section, the modified geodesic equations of photon in an external gravitational field are inves- 
tigated. The results obtained in the Section 5.3 are used in the subsequent section (Section 5.4) 
in the context of recent claims [127, 128, 137] that photons acquire superluminal velocities in 
the Schwarzschild , FRW and gravitational wave background. We extend the results obtained 
by [138] to the Kerr background and find, in concordance to previous results for other metrics, 
that in Kerr metric this does not happen. In the last section we derive the Papapetrou equations 
for the photon directly. 

5.3 Modified geodesic equation for photons 
5.3.1 Theoretical background 

The action of the electromagnetic field in the curved spacetime is given by: 

S = J d'x^F^,F>'\ (5.1) 

where F^j, has its usual meaningF^j, = Afj_.y-Ai,.f^. Varying the action with respect to we 
obtain the source equations for a free electromagnetic field in curved space-time: 

^^^^, = 0. (5.2) 

These equations together with the property of the field: 

V[^F,.] = (5.3) 

constitute Maxwell equations in an empty 4-space. To obtain a wave equation for a photon, 
we will apply to the Eq. 5.3: 

(V^F,;,) + (V.F,^) + (VaF^,) = . (5.4) 
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Using the properties of a commutator of the covariant derivatives 



[V'",V„]A,-, = V'"(VAx)-V.(V„A,-0 , 
and from the field equations (5.3) we have 

□/VA + [V^ V.]FA^+[V^ Va]/v. = . (5.5) 

V ^ ' V ^ ' 

® ® 

The operator □ will be frequently referred to. In general coordinates it is to be taken as defined 
by 

It is also sometimes written in the form (see [139], p. 64) □ = ((...).o,)'", that is, we perform 
a covariant and a contravariant differentiation and contract them. In some texts it means also 
□ = VV^, where V is a covariant derivative. 
Following the definition of a commutator of covariant derivatives 

[V„ =R\c,X\-R\,,X\ , (5.7) 

we expand the numbered terms in the equation (5.5): 
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Using symmetry properties of a Riemann tensor 



nP _ _nP nP _ uP _ nP 

we obtain the wave equation for a photon in a curved space-time: 

□F., +R\y,,F% = . (5.8) 

In empty space this equation becomes 



n^^A+^Vn = o. (5.9) 

This equation was first derived by Eddington [139], though not from the action. This equation 
is quite unusual, therefore it deserves some more explanations. Let us consider the 4-potential 
A^. Since it is not gauge invariant, we will add the standard Lorentz gauge: 

(A% = 0. (5.10) 

With the boundary condition at the infinity, the value of becomes completely determined. 
From Maxwell equations in the presence of a source, we get 

In the Lorentz gauge, we have 

In empty space this becomes 04^ = 0, showing that is propagated with the fundamental 
velocity. We see that the velocity of propagation of electromagnetic potential and of electro- 
magnetic force is not the same. It should not be surprising, since A^ is not physically important, 
being fixed by arbitraty convention (Eq. 5.10). Bu the result (5.8) is very interesting, showing 
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that Ricci and Riemann curvature coupling terms give rise to the polarization dependent devi- 
ation of photon orbits from null geodesies. The derivation of the result (5.8) from the action 
principle was first given in [138]. 

In order to derive equation for the characteristics of photon propagation we make the ge- 
ometrical optics plane-wave approximation. In geometrical optics, one assumes that the elec- 
tromagnetic waves have a wavelength which is much smaller than the radius of curvature of 
the background geometry, or the scale of variation of the amplitude of the wave front. We 
also assume in our analysis that the influence of the electromagnetic field on the metric field 
can be neglected. Thus, "locally plane" waves can be represented by approximate solutions of 
Maxwell's equations of the form 

iV- = e"^"VM., (5.11) 

where f^i, is a lowly varying amplitude, and s, a rapidly oscillating phase. (It is understood 
that one takes the real part of the right hand side of Eq. 5.11). We should emphasize that 
although this is an approximation for the actual solution for the photon, it exactly determines 
the characteristics of propagation and hence the casual structure of the solution, at least in 
domains where the gradient of the s does not vanish [14]. The wave vector is given by the 
gradient of the phase 

k^^V^s. (5.12) 

Operating by V^^ on (5.11) and neglecting the derivative of a slow changing amplitude, we 
obtain 

^nFaf3 = iFali^n s = ik^Faf} ■ (5.13) 

Using (5.11), Eq. 5.8 becomes 

Considering a constant vector, after some algebra we obtain the wave equation (5.8) gives 

-k.k^U +KU +R'\J,p-Kfpx = . (5.14) 
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Since, by definition, wave vector is a gradient, V^'k^ - V fji"', and 

Light rays are defined as the integral curves of the vector field k"', that is, k°' = dxf^/ds. 
Therefore, 



2— — +r" —1=2 



^ ds ds 



(5.15) 



We can write the Eq. 5. 14 in the form: 



and obtain the dispersion relation for the wave vector: 



- {R^\i,fij.p+R\fpu-Rlfpx) ■ 



f 



■vX 



(5.16) 



Combining (5.15) and (5.16), we obtain the modified geodesic equation for photons: 



ds^ ds ds 2 



^''^xufp^p + R\Spv ~ R%fpx ) 



r 



vX 



(5.17) 



The nonzero right hand side of the modified geodesic equation leads to a polarization depen- 
dance in the gravitational redshift, deflection angle and Shapiro time delay of light. In order to 
investigate these equations we consider the Schwarzschild space-time. 



5.3.2 Schwarzschild metric 

Schwarzschild metric is represented by: 



dr^=( dt^-(\-—\ ' dr'-r^d9^-r^sin^ed(p\ (5.18) 



r J \ r J 
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Equations of motion (Eq. 5.17) in an empty space are: 



® 



Covariant derivative of a Riemann tensor: 



xRaPnk ■ (5.20) 

After expanding ®, the right hand side of the Eq. 5.19 is: 

- ri,R,,.k) +Rp,ux {dar^+K,f^+r^,,r')] . (5.21) 

We will write this equation in components. 



A. Let a = (p. 

Due to the axisymmetry of the metric. 



After expanding V^, the right hand side (Eq. 5.21) is: 

Using the non- vanishing components of the affine connection and Riemann tensor, we 
find: 

6m(2m-r)sin2^/7''^ 
RHS = j^. 
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B. a^t 



and we use the fact that in the absence of charges daf"'^ = 0. After expanding V, the 
RHS for the time component t is: 

RHS = . 

C. a = r 

After expanding V^, the RHS for the r component is: 



RHS = 



IP 



Using non-vanishing components of the affine connection and Riemann tensor, we find: 



RHS = 



2m 1 



.(^fty^ + (^ft^ 2 (2r-5m) ^ ^ ^^^^^ 2 (2r-5m)sin 9 ^ 



+ f drf (2m -r) + f^' dj"^' (2m - r) sin^ 9] 

D. a = 9 

After expanding V^, the RHS for 9 component is: 



RHS = 



{rydeRtnr+ {ff deRtete + {ry deRt^,^++rR,nrder+ 



+ f^Rtift^ def"^-f"f^ (J^ee^trtr + ^%Rtete) - {f"^) ^ '^t^Rt^t^p 



Using non-vanishing components of the affine connection and Riemann tensor, we find: 



RHS = 



2m 1 



,2 J2 



{f'P'f3(r-2m)cos9sm9+fr^^^—^+ 



+ def - + def"^ (2m - r) sin^ 9 
r 
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For the equatorial plane, ^ = 7r/2, the equations are simplified: 
RHS for r 



J.2 p 



RHS for (/? 



RHS for 9: 



RHS for t: 



RHS 



6m(2m-r)ff^ 



■t6 str 



RHS = 



6m(r-2m)/V' 



RHS = 



5.3.3 The final set of equations 

Combining together left hand sides and right hand sides and noticing that we considered the 
amplitude to be nearly constant, the equations of motion are: 

Equation for r: 



■ + ^-. — - I — I ^ — -(r-2m) — -(r-2m)sm^^ I — 1 = 



\ds I r{r—2m) \ds I 



ds I 



\ ds 



2m 1 



6_ ^^,,^2 ^ (2r-5m) ^^g,^2 ^ (2r-5m)sin 6 



(5.22) 



Equation for 9: 



d^9 ld9dr . „ J d^V _lm 1 



ds^ r ds ds \ds ) |/^' 



^^^^^r/^ + 3(r-2m)cos^sin^ (/^)' 



(5.23) 
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Equation for ip: 



(fip^2drdip^^^^^^dedip_ 6m(2m - r) sin^ 9 ff"^ 



ds^ r ds ds 



ds ds 



(5.24) 



Equation for?: 



d t 2m dt dr 
ds^ r{r-2m) ds ds 



= 



(5.25) 



And for ^ = 7r/2, they become 



[a] 
[b] 
[c] 
[d] 



JV_ _^2drdip _ 6m(2m - r) f"'f"f 



ds^ r ds ds 

d^t 2m dt dr ^ 

— r + = 0; 

ds'^ r(r- 2m) ds ds 

d^e _ 6m(r-2m)f'^f\ 
d?~ ? 7^' 



dh 



m 



ds^ r{r—2m) \ds 



dr 



■(r— 2m) 



2m 1 

-prj2 



d(p\'^ m{r—2m) f dt^ ^ 
ds J r^ \ds , 

2 (2r-5m) 



+(r) 



(5.26) 



Eqs. 5.26 give the following conserved quantity, the energy E: 



r-2m dt 
r ds 



= E 



(5.27) 



However, we notice from Eqs. 5.23 and 5.24 that the direction of the angular momentum is 

not preserved; photon that started at, say, 6 = 11/2, will not remain in the equatorial plane — the 
right hand side of the Eq. 5.26[c] gives the acceleration to a photon in the 9 direction. 



5.4 Velocity of photons in Kerr metric 

In this section we come back to the results of the Section 5.3 in view of the claims [127, 
137] that in higher derivative gravity which arises by QED radiative corrections, the photon 
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velocity in local inertial frame can exceed the velocity of light in the Minkowski background. 
This result was extended by Daniels & Shore to charged [128] and rotating black holes [129] 
with the same conclusions. However, Mohanty & Prasanna [138] showed that the claims of 
superluminal photon velocity are due to the neglect of the curvature coupling terms in the wave 
equation derived from the minimal F^i.F^" Lagrangian (see Section 5.3, Eq. 5.8) and that the 
velocity of photons is less than c in Schwarzschild and FRW metrics, as well as in a higher 
derivative gravity, (as long as the SEC (see Section 3.1) is satisfied). We extended these results 
to a Kerr gravitational background. 



5.4.1 e-matrix 

To obtain the photon velocity we will use the dispersion relation (Eq. 5.16). Due to the in- 
ternal Maxwell equations (5.3), only three components of F^,^ are independent. Choosing the 
components of the electric field vector = /o, as the independent components, we have from 
(5.3) and (5.13 

kofij+kifjo + kjfoi = 0. (5.28) 

Using (5.28) to substitute for fij in terms of electric field components fjo and /o,- in the wave 
equation (5.14), we obtain 



(5.29) 



Rewriting it through only one component, foj, gives 



k„k^+Rl+R\,'^ ] 5i+ { R{-RI^-1R'L-AR'^ ^ 



Oi 



foj = 



(5.30) 



This is the wave equation obeyed by the three components of the electric field vector. We can 
simplify this equation by defining the matrix e/: 



ko 



R"o+R'o^\S^+iR'-R 



■i ^"-0 



ko 



-2R%-4R'^ 



(5.31) 
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and thus, (5.30) becomes 

leSi + ei)foj = 0. (5.32) 

We can rewrite this equation as matrix equation and by diagonalising the ej, we find the equa- 
tions for normal modes: 

4foj = -k^foi , 

and 

{k^ + ei)foi = (/= 1,2,3), (5.33) 

where e, are the eigenvalues of the e^j matrix. By defining the projection operator Pjk - Sjk - 
HjUk, with n' = k'/\n\, these equations can be decomposed into a transverse part 

f^]^={5'j-n'nj)foi, (5.34) 

and a longitudal part 

f^f = n'njfoi. (5.35) 

By substituting back into (5.33), we obtain the wave equations for the transverse photons 
(polarization) 

[k^ + q] (5/-nX)/o, = 0, (5.36) 
and equation for the longitudal photons (polarization) 

[e + ei\n^nifoj^O. (5.37) 

5.4.2 Kerr metric 

First we have to find the components of the ej matrix. Kerr metric is represented in the form 

ds^^p^—(dtf--rsmHid<^ :^dt] -^{drf-p\def , (5.38) 

E \ y A 
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where 



A = r^ + a^-2Mr; 



(5.39) 



and parameters a and m have dimensions of length. Full form of e-matrix for the Kerr metric 
is presented in App. C. Here we will work in the linear in a approximation. Then, in the 
coordinate frame , r, ^, (^} the e-matrix is reduces to 



-^^cot^ 



COS g sing 



-^{r^ + 3ak^) - ^'^^y^'"^ cos 6 sine 



\ -'-^(rh-2k2cot0) ^(k2+h(r-2m)cote) 



6am , 



2m 
'73" 



(5.40) 



In the equatorial plane, ^ = 7r/2, the deviation from the geodesic is expected to be maximal. 
Thus, 

/ ^{r^ + 3aks) ^ 



Eigenvalues here are: 





llamki 
,6 



.2m f,,2 
r 



{r^ + 3ak3) 



6amki 
r5 



2m 

" r3 



(5.41) 



, 4m .2 T T \ 
Xi = — [r +3ak3) ; 

Ai = ^ (r^ + 6aA;3) ; 

A 



(5.42) 



Considering radial trajectories with n = (r,9(f) = (1,0,0), equation (5.36) yields the wave 
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equation for the transverse fields EiiEy. 



P-^ \ / ^2 
^2_2m I I 



= . (5.43) 



The dispersion relation yields 



.2 



2m 



^i = -r- (5.44) 



The velocity of propagation of the transverse photon is 

V=^= = . ' <!• (5.45) 

Thus, the photon velocity is subluminal. It becomes luminal (i.e. equal to c), when either m = 
(no black hole), or at large distances r ^ oo. 
For tangential trajectories: 

(i) n = (0, 1,0), and the equations for the transverse fields Ei and £'3 yield: 



k^ + ^ 



k^-^ 





(eA 







= . (5.46) 



The dispersion relation is {k^ + ^) {k^-^) and the root, corresponding to the propagating 
mode, is 



uj^-k^ = — . (5.47) 



The velocity of propagation is 



v^ = ^====<l. (5.48) 

The velocity is subluminal again. 

(ii) n = (0, 0, 1), and the equations for the transverse fields Ei and E2 yield: 



k'^ + ^{f^ + 3ak3) 

k^-^{r^ + 6ak3) 
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(eA 




[eJ 



= . (5.49) 



The equation for the propagating mode is 



(5.50) 



and the velocity of wave propagation is given by 



^3 + — (6aA:3 + r^) 



1 . 6am 



1+1:^(1+^) 



1/2 



(5.51) 



In this case, the condition for the velociy to be subuminal reduces to: 



18a m 



< 1 . 



For < m^, we know that the null horizon is at r± = m ± (m^ -a^Y^^. In this case, for r+, the 
condition goes to ~ ISa^m < (2m)', which is satisfied. For the extreme case, a = m, the nul 
Ihorizon is at r± = m, and we have IS/m^ < 1. 

We conclude that velociy of photons is always subluminal in the field of the Kerr black 
hole. 



5.5 Derivation of the Papapetrou equation for photons 

According to the EP of the general theory of relativity, the motion of structureless test parti- 
cles in a gravitational background field is determined only by the spacetime geometry: particle 
worldlines are the geodesies of the spacetime. Things become more complicated for test par- 
ticles which are not structureless but carry, for example, a non-vanishing charge or spin. In 
such cases, the worldline of the test particle, in general, is no longer a geodesic, but is modi- 
fied by electromagnetic and/or spin-gravity forces (see [140, 141] and references therein). It is 
far from obvious whether one can observe, in practice, the spin corrections to the equations of 
motion of elementary particles. However, the problem of influence of the spin on the trajectory 
of a particle in an external field is not only of pure theoretical interest. Spin-dependent correc- 
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tions certainly exist in differential cross sections of scattering processes. It was proposed long 
ago that it is posible to separate charged particles of different polarizations through the spin 
interaction with external fields in a storage ring of an accelerator [142]. Though this proposal 
is being discussed rather actively (see review [143]), it is not yet clear whether it is feasible 
technically. The EP can be put to test in an astrophysical setting. A recent proposal is based on 
the analysis of the differential time delay between the arrival of left and right-handed circularly 
polarized signals from the millisecond pulsar PSR 1937+214 [144]. The authors have reported 
the measured difference. However, far few papers exist on the theoretical foundations of pos- 
sible deviations of the photon motion from geodesies (see, for example [138, 144]). Here we 
report yet another approach to this long-standing problem, based on the Lagrangian approach. 
If an appropriate Lagrangian density is taken into account, then the photon equation of motion 
(modified geodesic equation) can be found from the Euler-Lagrange equations. 

5.5.1 Formalism 

The generalised concept exists that classical particles follow the path of least distance between 
the endpoints, even in curved spacetime . Thus, in (pseudo) Riemannian spaces geodesic equa- 
tion is found by variation of an action S identified with the parameter 5 of a curve, interpreted 
as its length. The same method has been used to find the geodesic equation for light where, 
however, one technical problem arises: photon's worldline is null. One way to avoid this dif- 
ficulty is to consider motion of a massive vector particle, then, if the equation obtained does 
not contain mass explicitly, simply put the mass to zero, with one further condition that the 
four velocity of the particle be a null vector. This method has been widely used for the scalar 
particles and is known to give the equation of geodesic regardless of mass. In this work we 
apply this approach, which makes it possible to use the length parameter in our considera- 
tions of the action principle. Usually in order to describe the behaviour of a field in a given 
gravitational background one solves the corresponding field equations for the given metric. If 
the goal is to describe waves, one can take the corresponding wave solution. In a spherically 
symmetric spacetime the solutions contain factors expressed in terms of spherical harmonics. 

135 



This method works well when the wavelength is comparable to the scales under consideration. 
It is not so in the case of light, propagating in the vicinity of a massive object and, to describe 
the propagation of light as an electromagnetic field, one would have to employ the spherical 
harmonics of a very high order. The corresponding solution would look too complicated and 
tell little about the behaviour of light. 

Many efforts have been spent in the last decades to work out a simple approach to this 
problem. The idea is find a satisfactory approximation to the wave as some curves, that could 
be called "rays", and, at the same time, would take into account the polarization of light. In this 
section we work out a simple approximation of this type. Our idea consists of the following: 
we consider a massive vector field obeying the Proca equation, describing the propagation of 
this field in some restricted domain of spacetime. The shape of the domain can be chosen as 
that of some world-tube transverse to the wave, with the cross-section comparable to or no 
more than two orders of magnitude greater than the wavelength. As this tube is timelike there 
must exist a timelike curve A in its interior, which specifies the time axis of a local coordinate 
system. If the tube is not too wide, this coordiniate system would cover the entire interior. If s 
is the proper time on the curve A, the curve can be chosen in such a way that the field equation 
reduces to 

^ + m% = 0, (5.52) 

the same way as it happens in standard Cartesian coordinates for Minkowskian spacetime. Cor- 
respondingly, the field Lagrangian is -A^+m^A^, where dot stands for the covariant derivative 
on s, if and only if the curve A is chosen properly. This Lagrangian should contain one more 
term, responsible for the shape of A, such that we get the geodesic equation, when we switch 
off the field. The form of this term is well known: 1 /2mx^, thus our final Lagrangian is 

2L = mx^ -A^+ rn^A^ , (5.53) 

and coupling between the field and the shape of the curve A is incorporated in the form of a 
covariant derivative A, which contains the product of connection, velocity x and the field. 
The derivation of the conservation laws is more convenient in orthonormal frames. In what 
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follows, el will denote the components of an orthonormal 1-fomi frame field, 



= e1(x)dx' , 



(5.54) 



and the components of its dual vector frame field. 



(5.55) 



Here frame indices are always a,b,c,...; coordinate indices are i, j, k,.... The metric tensor 
can be expressed as 

g = gijdx' (g) dxj = SabO" ® . (5.56) 

The connection 1-form for these frames may be introduced through the first structure equation 
(see App. D): 

J^° = a;/A^\ (5.57) 
and the connection coefficients 'jabc are that of the expansion of this 1-form in the local frames 

< = iJO" . (5.58) 



Thus, 



(5.59) 



or A" is a co variant derivative in orthonormal frame on a curve x'(s) with 7^^" being a spin 
connection. 

5.5.2 Equations of motion 

Each generalized coordinate has its conjugate generalized momentum: 



.h :b 9 



dA' 
ds 



dL 
dA" 



= mr]abx''-VdiA'A'%c , (5-60) 

(5.61) 



137 



The Euler-Lagrange equations are 



d_ _ dL 



d dL 
—E„ = 



ds " aA« 



Let us first consider equations for the Ea (5.63). RHS of this equation is: 



dL d 1 



aA^ dA" 2 L 



- ^ ^bcQ^A +m ^^A - 



d 
dA^ 



dA^ 
ds 



b-c Ad 



+m'A'riac = -A'7]bc%a^ + m^A'r]ac ■ 



And Euler-Lagrange equations for are therefore: 



ds 



Ea+A'r]bc%aX' -m^A'riac = 



With the help of Eq. 5.61, 



j-E,-E,jJx'^]-m^A'rjac = 



Expression in the brackets is a covariant derivative for Ea, thus, we obtain: 



Ea-m^A'r]ac = 



Again using Eq. 5.61, we finally obtain the Euler-Lagrange equations: 



A''r]ab+m^A%b = 0, 



which, as we can see, reduces to Proca equation for the four-vector field A^(x) [145] 



{D + m^)Ap^O 
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We describe spin of the particle in this model directly by a tensor of spin Sab- 

mSab = ^ (AaAb-AbAa) ■ (5.70) 
Thus, spin enters the expression for the generalized momentum (Eq. 5.61): 

Pa = Vabm^-%^Sj', (5.71) 
To obtain the equation of motion for spin, we write 

^ = \i [AbA^-AAb] - [AbA^-AAb] . (5.72) 



Due to the wave equation (5.69) all derivatives of A vanish and we are left with: 

dSdh 1 .„ r / . .-.id 



ds 



= \x^[AAi-AjA]^^[e^A) = 



[AAj-AjA^] {e^dlcbei + e\^Je\) = 

(tJ [AAb-AbA,] [AAa-AbaA,]) . (5.73) 

Here is the matrix introduced in the equation (5.55) and derivatives (e^) ^ are obtained from 
the first structure equation (see App. D for the derivation). Using our definition of a spin tensor 
(Eq. 5.70), we obtain: 

'^^"'=x^{jJS^b + lJS^) , (5.74) 



ds 
or 

^ = 0. (5.75) 
Ds 

Thus, spin is transported parallel to itself along the worldline. 

To derive Euler-Lagrange equations for the generalized momentum, we return to the coordinate 
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basis. In this case Lagrangian (Eq. 5.53) takes the form: 

lL = mgij£xi-gij (^ + V[/A^ {^ + Ti/A^ +m^gijA'A^ , (5.76) 

where T^i are the Christoffel symbols. This gives the generaUzed momentum: 

8L 

Pi=g^i= mgijxj - TT^A'A^ . (5.77) 

The Euler-Lagrange equations: 

^P,= ^. (5.78) 

as ax' 

LHS of this equation is: 

where derivatives 9, are defined as 9, = RHS is: 
31 \ ■ ■ 

— = -mdigjkxjx' + digmn {-A'-A- + m'A"^A-) -gmn^^A^A'^diT^ . (5.80) 
ox' I 

From the definition of the Christoffel symbols 

digmn = \ {gkn^ii + gkm^i) ■ (5.81) 



Using 



we obtain: 



x'xJdjgik = ^x'xj {djgik + digjk) , (5.82) 



mgtj — + -i'xv + - %,-,.) - m^r^^-A^ - VIA' + 

r:;^'"A„-r-A,— -xM'A,a,r^+g,„/A'A"a,rr, = . (5.83) 

as 
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First two terms of this equation are the covariant derivative of mf, and using the Proca equation 
(5.69) after some algebra we get: 



With the definition of the curvature tensor and spin tensor (5.70), the equation takes the evi- 
dently covariant form: 

8ij^=R'jii^%'- (5.85) 

This equation coincides with the Papapetrou equation in the linear-in-spin approximation 
[133]. It must be pointed out that in case of a zero spin this equation becomes a geodesic. 
If we reparametrise the curve with some new parameter A in such a way that 

dx' dx^ 

we can rewrite Eq. 5.85 as 

8ij^=R'ja^S^ ■ (5.87) 

This equation is valid for the massless particles with spin as well, since mass does not enter 
the equation explicitely. 



5.6 Conclusions 

In this chapter we considered the propagation of photons in the strong gravitational field. 
Strong field can have additional effects on the light propagation which do not exist (or are of 
vanishingly small order) in the weak-field approximation. Gravitational lensing of light in the 
field of black holes (and, may be, neutron stars) is strongly influenced by these effects and 
the basic theoretical formulation underlying it is necessary. However, contrary to previous 
claims, we confirmed that the photon's speed does not become superluminal even in the field 
of the Kerr black hole. If, in fact, photons do not follow geodesies in the fields of Kerr and 
Schwarzschild black holes, as can be followed from the Papapetrou equation, or they do ex- 
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perience gravitational birefringence, all the calculations for the gravitational lensing in these 
metrics done so far have to be revised. The main aim of our future work, following the narra- 
tion of this chapter is to solve the equations of motion of photons in the fields of Schwarzschild 
and Kerr black holes and obtain the trajectories for RCP and LCP. We expect to prove that, 
even in the first order approximation to the geodesic, the polarizations separate. We also plan 
to apply the results of Sec. 5.4 to the phenomenon of gravitational Faraday rotation in of pulsar 
radiation. 
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Chapter 6 



Summary and Final Comments 

Over the last few decades the status of gravitational lensing has been raised from that of a 
theoretical speculation to that of a precision measurement. Indeed, gravitational lensing is 
now widely used to establish concordance of a given cosmological model. In this thesis we 
have applied gravitational lensing in the above spirit. 

Gravitational lensing provides an exciting new probe for detecting exotic objects in the 
universe that may be described by matter fields which are not yet detected. Proposals have 
been made to look for signatures of cosmic strings [25], boson stars [146] and neutralino stars 
[147] from their gravitational lensing effects. There is no compelling evidence that any of the 
observed GL systems are due to these objects, however, it is essential to develop new lens 
models with objects which are not forbidden on the theoretical grounds. In chapter three we 
have explored the consequences of the existence of matter violating the energy conditions. 
We develop the formalism of gravitational lensing by a point negative mass lens. Negative 
mass microlensing simulations are carried on, showing the resulting shapes of the images, the 
intensity profiles, the time gain function, the radial and tangential magnifications, and other 
features. We provide an in-depth study of the theoretical peculiarities that arise in negative 
mass microlensing, both for a point source and for extended source situations for various 
source intensity profiles and discuss the differences between them. We have also numerically 
analyzed the influence of a finite size of the source on the negative-mass gravitational lensing 
event. 
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Since the primary motivation for our theoretical work is observational, in order to set the 
experimental context more explicitly, we take wormholes as a useful and appealing theoretical 
scenario for the appearance of negative masses. But it applies, mutatis mutandis, to any pos- 
sible object possessing a negative mass. Thus, we provide a useful comparison arena where 
to test the possible existence of wormholes or other kinds of negative mass objects. Our work 
gives a method of analyzing observational predictions quantitatively. The next step would be to 
test these predictions using archived, current, and forthcoming observational microlensing ex- 
periments. The only search done up to now included the BATSE database of 7-ray bursts [62]. 
There is still much unexplored territory in the gravitational microlensing archives. Follow- 
ing our suggestions, the MOA (Microlensing Observations in Astrophysics) group is currently 
adapting their systems to look for the effects of negative mass microlensing in the first direct 
online search for exotic objects in the universe [148]. 

It has been recently suggested that there can exist large classical violations of the energy 
conditions [76]. Disregarding the fundamental mechanism by which the EC are violated (e.g. 
fundamental scalar fields, modified gravitational theories, etc.), if large-scale localized vio- 
lations of null energy condition exist in our universe, we show that it would be possible to 
detect them through cosmological macrolensing. Contrary to the usual case, where arc-like 
structures are expected, finger-like "runaway" filaments and a central void appear. Compari- 
son of the effects of negative masses on background fields with macrolensing effects produced 
by equal amounts of positive mass, located at the same redshift, brings out notorious differ- 
ences. These results make the cosmological macrolensing produced by matter violating the 
null energy condition observationally distinguishable from the standard situation. Whether 
large-scale violations of the EC, resulting in spacetime regions with average negative energy 
density indeed exist in the universe can therefore be decided through observations. 

Alternative cosmologies have been more and more actively discussed in the literature fol- 
lowing serious theoretical and observational difficulties with standard cold dark matter FRW 
cosmology [78]. An open FRW model with linear evolution of the scale factor, a(t) oc t, has 
generated a lot of interest. Many of the problems of a standard model are naturally resolved 
in such a cosmology. This model is consistent with nucleosynthesis in the early universe [93] 
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and presents a good fit to the current SNE la data [94]. It is now well known that gravitational 
lensing is a useful probe of the geometry of the universe. It has been pointed out [96, 99] that 
the expected frequency of multiple imaging lensing events for high-redshift sources is sensi- 
tive to a cosmology. In chapter four we used this test to constrain the power index of the scale 
factor of a general power law cosmology, a{t) oc t". Expected number of lens systems depends 
upon the index a through the angular diameter distances. By varying a, the number of lenses 
changes and on comparison with the observations gives us constraints on a. We describe the 
general formulation of the gravitational lensing statistics and use it to constraint the parameters 
for that cosmology. These results lead to the conclusion that open linear coasting cosmology, 
a(t) = t, is consistent with present observations. 

Linear evolution of the scale factor is a generic feature in a class of models that attempts 
to dynamically solve the cosmological constant problem [92]. Such models have a scalar field 
non-minimally coupled to the large scale curvature of the universe. With the evolution of time, 
the non-minimal coupling diverges, the scale factor quickly approaches linearity and the non- 
minimally coupled field acquires a stress energy that cancels the vacuum energy in the theory. 
Such an evolution of the scale factor is also possible in alternative effective gravity and higher 
gravity theories. It is also possible in the "toy" model [107] that combines the Lee-Wick con- 
struction of non-topological soliton (NTS) solutions in a variant of an effective gravity model 
proposed by Zee [149]. The conditions of the theory are sufficient for the existence of large 
NTS's with the size equal or larger than typical halos of galaxies. The interior and exterior 
of such a domains would be regions with effective gravitational constant Gg^ = constant and 
= 0, respectively. Such a domain is characterized by an effective refraction index in its 
interior, which would cause bending of light incident on it. In chapter four we examine the 
possibility of a special kind of NTS, gravity ball, to be a gravitational lens. We investigate its 
lensing properties, calculate the deflection angle of a point source, derive the lens equation and 
plot the lens curve for the gravity ball with specific parameters. We have shown that depending 
on the parameters of the model there can be one, two or three images (two inside the Einstein 
radius and one outside). For the extended sources we showed how a large gravity ball can 
induce surprisingly large distortions of images of distant galaxies. It is assumed that in order 
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to produce large arcs a cluster must have its surface-mass density approximately supercritical, 
S > T,crif We found that the empty balls alone are able to produce arc-like features, including 
straight arcs, arclets and radial arcs. We have modeled the gravitational lensing effect of our 
gravity ball on a background source field and found that the ball produces distortion ('shear') 
of that field, consistent with observations. We also obtain constraints on the size of a large 
gravity ball from the existing observations of clusters with arcs. Empty G-balls can show in 
observations as arcs without evident mass distribution — a lens and/or through the distortion of 
a background random galaxy field. Observationally empty G-balls have not yet been found, 
but we conclude from our analysis that the existence of this special kind of NTS, large G-balls, 
cannot be ruled out by gravitational lensing effects. 

Most cosmological studies of gravitational lensing are performed in the weak field approxi- 
mation. However, weak field approximation becomes invalid in the vicinity of compact objects 
like black holes and pulsars. Propagation of light in such backgrounds must be described as 
occurring in strong-field regime. Chapter five is dedicated to the studies of the propagation of 
light in the fields of Schwarzschild and Kerr black holes. Strong fields can exhibit additional 
effects on the light propagation which do not exist (or are vanishingly small) in the weak-field 
approximation. Gravitational lensing of light in the field of black holes (and, may be, neutron 
stars) is strongly influenced by these effects and it is necessary to develop the basic theoretical 
formulation underlying it. Contrary to previous claims [127, 137], we confirm that the pho- 
ton's speed does not become superluminal even in the field of the Kerr black hole. We derive 
the modified geodesic equation for photons in Schwarzschild background and present a way to 
solve it. Alternatively, we derive the Papapetrou equation for photons in the first order approx- 
imation. If, in fact, photons do not follow geodesies in the fields of Kerr and Schwarzschild 
black holes, as can be followed from the Papapetrou equation, or they do experience gravita- 
tional birefringence, gravitational lensing calculations in these metrics done so far have to be 
revised. This is a part of our ongoing research work. 
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Appendix A 

Wormhole solution to the Einstein 
equations 

Building of the wormhole solution of the Einstein equations was done in quite a non-traditional 
way, that is why we will describe it in details. The traditional way of obtaining solutions to 
Einstein equations is to start by picking one's favourite Lagrangian for the matter fields that 
could support the wormhole's spacetime. One then calculates the relevant stress-energy tensor 
and solves the Einstein field equations. Finally, one checks for the presence of curvature singu- 
larities in the solution so obtained — the goal is to avoid any singularity whatsoever (like event 
horizon or naked singularity), since the main idea is to obtain a traversable stable wormhole. 

Till the seminal paper of Morris and Thorne [54], all attempts at following the above de- 
scribed prescription failed. The resulting wormholes were pathological in one or more man- 
ners. Morris and Thome adopted a more engineering oriented approach. First, assume the 
existence of a suitably well-behaved geometry. Then calculate the Riemann tensor associ- 
ated with this geometry and use the Einstein field equations to deduce what the distribution of 
stress-energy must be. Finally, decide whether the deduced distribution is physically reason- 
able. 
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Desired properties of traversable wormholes 



It is interesting to note that wormholes as objects of study in mathematical relativity predate 
black holes: Within one year after Einstein 's final formulation of his field equations, the 
Viennese physicist Ludwig Flamm recognized that the Schwarzschild solution of Einsteiiis 
field equations represents a wormhole [150]. Possible roles of the Schwarzschild and other 
wormholes in physics were speculated upon in the 1920s by Herman Weyl [151], in the 1930s 
by Einstein and Rosen [152], and in the 1950s-1960s by John Wheeler [153]. However, all the 
above described configurations were unsuitable for the use of travel. Here we will state the 
properties we would like the wormholes to possess in order to be traversable — a Lorentzian 
wormhole that is (quasi)-permanent and macroscopic so as to allow for traverse by human 
beings. We will state them verbally without the mathematical details. 

(i) The metric should be both spherical symmetric and static. 

(ii) The solution must everywhere obey the Einstein field equations, since we assume the 

correctness of GR theory. 

(iii) To be a wormhole, the solution must have a throat that connects two asymptotically flat 
regions of spacetime. 

(iv) There should be no horizon, and tidal gravitational forces experienced by a traveler must 
be bearably small. 

(v) A traveler must be able to cross through the wormholes in a finite and reasonably small 
proper time as measured not only by a traveler, but also by observers outside the worm- 
hole. 

(vi) The matter and fields that generate the wormhole 's spacetime curvature must have a 
physically reasonable stress-energy tensor. 

(vii) The solution should be perturbatively stable (especially as a traveler passes through). 

(viii) It should be possible to assemble a wormhole. For example, the assembly should both 

require much less than the mass of the universe and much less time than the age of the 
universe. This is closely related to the property (vi), in the sense that we shall make 
the wormhole construction material as compatible as possible to our present prejudices 
about the forms of matter allowed by physical laws. 

148 




rcos0 



Figure A.l: Left: Embedding diagramme for a wormhole that connects two different uni- 
verses. Right: Embedding diagramme for a wormhole that connects two distant regions of our 
universe. Each diagramme depicts the geometry of an equatorial slice {0-t[/2) through space 
at a specific moment of time {t = constant). 

Thus, adopting Schwarzschild {t,r,9,(p) coordinates, the metric can be put into the form 

ds^^-e^^^^''^dt^ + —4^—— + r^\dd^ + sm^ed(t?] . (A.l) 
\-b±{r)/r '- -' 

• Radial coordinates: r, such that Ivrr = circumference; / — proper radial distance from the 
throat. They are related by /(r) = ± j^^dr/{\ -b±(r')/r'y^/^. Two coordinate patches 
([ro,+oo]), each covering one universe, join at tq, the throat of the wormhole. 

• Solution is determined by two freely specifiable functions of r: shape function, b(r), 
defined by dl / dr = ±(l -b/r)~^/^ and redshift function, 0(r), defined by gn = -e^'^. The 
mass of the wormhole, as seen from spatial infinity, is given by b± = 2GM±, the two 
mouths of the same wormhole can in general have different masses; besides, there is no a 
priory requirement that 0+(oo) = 0-(oo), which implies that time can run at different rates 
in the two universes. Constraints on b(r) and 0(r) to produce a traversable wormhole: 

(a) Spatial geometry must have wormhole shape (see Fig. A.l). No horizons or sin- 
gularities, (j) is everywhere finite; t measures proper time in asymptotically flat 

regions <^=^ ^ 0, as / — > ±00. 

(b) Radial velocity of a traveler as she passes radius r, as measured by a static observer 
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there, in temis of distance traveled dl, radius traveled dr, coordinate time lapse 
dt and proper time lapse as seen by the traveler, dT, is v = dl/e'^dt = ^dr/{\ - 
hjrfl^e'^dt. Defining 7 = [1 -(v/c)^]"^/^ we find 

dr 

'^ = ^ il-b/ryndr - ^^-'^ 
The material that generates the wormhole's spacetime curvature: 

(a) Stress-energy tensor as measured by static observers in the orthonormal frame: 

Tfi = pc^ = (density of mass-energy); Tff = -r = -(radial tension); 

Tgg = T^^ - p- (lateral pressure) . (A.3) 

(b) Einstein field equations: 

b' _ b/r-2(r-b)(j)\ r 



T = 



p = -[(pc^-r)(t>'-r']T. (A.4) 



87rGc-2r2' SttGc-V^ ' ^ 2 

(c) The constraints placed on the shape function b(r) give rise, via the Einstein field 
equations, to constraints on p, r and p, that generate the spacetime curvature. The 
most severe constraints occur in the wormhole's throat. The requirement that the 
wormhole be connectible to asymptotically flat spacetime leads to the flaring-out 
condition 

^"^i > at or near the throat. 

The constraint tq > poc^ means that in the throat the tension of the wormhole's ma- 
terial must be so large as to exceed the total density of mass-energy poc^. Material 
with this property, r > pc^ > 0, is called "exotic". An observer moving through the 
throat sees an energy density (projection of the stress-energy tensor (Eq. A.3) on 
traveler's time basis vector eg, = 76^ =f 7(v/c)er) ) given by 

Tq'o' = 7^(Poc^ - ^0) + To • (A.5) 



If such an observer moves sufficiently fast, with a radial velocity close to the speed 
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of light, 7 1, he will see a negative density of mass-energy. A rough explanation 
of th negative mass-energy density seen by some observers is that bundles of light 
rays (null geodesies) that enter a wormhole at one mouth and emerge from the 
other must have cross-section at areas that initially decrease and then increase. The 
conversion from decreasing to increasing can only be produced by gravitational 
repulsion of matter through which the rays pass, a repulsion that requires negative 
energy density. 

The above material was compiled using Refs. [12, 54]. 
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Appendix B 

Lensing clusters sample 



Table 1 : Lensing clusters sample 



cluster 


Zcl 


arc" 


■Zarc 


r (arcsecs)'' 


ref." 


A2280 


0.326 


large arc 


RiO.7 
or 1.0? 


14" 


1 


MS 0440+0204 


0.190 


1 GLA + 












3 bright arcs 


0.53? 


22" 


2 


A370 


0.374 


arclet A5 


1.306? 


63.7" 


3 


PKS 0745-191 


0.1028 


bright arc 


0.433 


18.2" 


4 


A2390 


0.231 


GLA; straight 












image 


0.913 


37" 


5 


A963 


0.206 


GLA+counter 




image separation= 








arc 


0.77 


30"; ^E= 15" 


6 


CI 2244-02 


0.329 


GLA 


2.237 


14" 


7 


CI 0024-1-1654 


0.391 


five arcs 


1.39? 


35" 


8 


CI 0500-24 


0.316 


GLA 


0.913 


~26" 


9 


S295 


0.299 


GLA 


0.93 


~20" 


10 



" Lensing phenomenon. GLA — giant luminous arc. 
* Distance from the centre of cluster where arc(s) are detected. 
Only one reference is provided for lensing feature. References are listed in [154]. 
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Appendix C 



Full matrix in t^ip^r^ 6 components 

SMr (-3^2 + 2r2 - 3a^ cos 29) 
{a^ + 2r^ + a^cos2eY 

8Mr (-3fl2 + 2r2 - 3a^ cos 29) {5a^ + 6aki +4r^-a^ cos 29) 
{a^ + 2r^+a^cos29)^ 

1 6Mr (-3a^ + 2r^- 3a^ cos 26') (2<32 + 3aki + r^- cos 26') 
(a2 + 2r2 + a2cos2^)'* 

%ahM{a^ + r{-2M+r)) {a^-6r^+a^cos29)sm29 _ 
(a2 + 2r2+a2cos2e)'* 

8aM (^6k2r{-3a' + 2r^-3a^cos29)- ''^-'"'''''''-''^^^^^^ 

{a^ + 2r^+a^cos29)'^ 

Sak2M{a^ + r(-2M+r)) {a^-6r^ + a^cos29) sin 2^ 
(a2 + 2r2 + a2cos26')'^ 

16aM{3k3r(3a^-2r^ + 3a^cos29) + k2 {-2a^ + (2M-r)r+a^cos29) {a^-6r^+a^cos29)cot9) _ 

(fl2 + 2r2 + fl2cos2^)'* 

8aM (6r2 - - cos 26*) {-3a^ - 5a% + a{4M -3r)r + 4ki (2M - r)r + a{3a^ + aki + 3r^- AMr) cos 26*) cot 9 

(a2 - 2Mr + r^) {a^ + 2r^ + cos 2^) ^ 

8aM (6r2 -a^-a^ cos 26*) (-3fl^ - 4a2A;i +fl(2M- 3r)r + 2A:i(2M- r)r+a{3a^ + 2aA:i + 3r2 - 2Mr) cos 26*) cot 9 

(a2 + 2r2 + fl2cos2^)^ 
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Appendix D 
Structure equations 

D.l Connections 

Consider a Riemann space M endowed with a coordinate system {x'} and let 1 -forms 

= e'^(x)dx^ (D.l) 

constitute an orthonormal frame in every point of the space M, such that its metric tensor may 
be expressed through the Kronecker delta: 

g = gijdx' (g) dx^ = Sabi^" (8) ly''. 

Any given connection is characterized by two tensors, torsion and curvature. The torsion can 
be expressed as a vector-valued 2-form T" (we usually assume zero torsion) and the curvature 
as a (l,l)-tensor-valued 2-form Q^. The defining relations for these two tensors are known as 
Maurer-Cartan structure equations: 

1st Structure Equation: T'' = du'' + A (D.2) 
and 

End Structure Equation: Q.^ = dco^ + co^ A u;/ , (D.3) 
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where is called the 1-form of connection. Due to the orthonormality of the frame, the 
connection 1-form satisfies the condition 



and the connection coefficients ^yabc are that of the expansion of this 1-form in the local frames 

< = 7>^ (D.5) 

The Riemann curvature is represented by the 2-form Q^, defined by the second structure equa- 
tion, and coefficients of its expansion in the local frame of 2-forms l^u" Au''} constitute the 
components of the Riemann curvature tensor 



D.2 Derivatives of the components of the frame 

Suppose, we have an orthonormal frame of 1 -forms 

= e^idx' , (D.7) 

and dual to it vector frame 

dx' = e'y. (D.8) 
From Eqs. D.2 and D.5, we write 1st structure equations as (torsion is zero) 

du" + -f Ju' Au^ = 0. (D.9) 

Using (D.7), we write 
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where the derivatives of the tetrad components can be written by analogue with a gradient of a 
function, df = fidx' = {£f) dx' = (e^ of)v' = f^,v', as 

de't={e,oef^v^={e^) y . (D.IO) 

Thus, with the use of (D.8), 

du" = {el) ^ A dx' = {e'{) ^ e\u' A . 
Then, we can write (D.9) as 

{e1)/,v'^u' = -^,^u'^vK (D.ll) 
The 1st structure equation can then be written as 

{eD^, = -%Se\. (D.12) 
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